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CHAPTER  1.  INTRODUCTION 


1.  General  -  The  System  Description.  The  System  Description  (SD)  for 
the  European  Theater  Network  Analysis  Model  (ETNAM)  has  been  prepared 
for  the  Defense  Communications  Agency  (DCA)  under  contract  DCA  100-70- 
C-0039,  Modification  $G01,  dated  January  22,  1970*  Submission  of  this 
document  is  in  fuxfillment  of  paragraph  3*2  of  the  contract. 

a.  The  System  Description  consists  of  four  related  volumes  bound 
separately.  Volume  I  is  the  User's  Manual  and  is  intended  primarily 
for  use  by  the  action  officers.  Volume  II  is  the  Operator's  Manual 
and  was  written  by  the  National  Military  Command  System  Support  Center 
(NMCSSC)  and  is  intended  for  the  use  of  computer  operators  and  others 
interested  in  the  operating  characteristics  of  the  model.  Volume  III, 
this  document,  is  the  Analytical  Manual  and  is  intended  for  the  use  of 
analysts  and  programmers.  Volume  IV  is  the  ETNAM  Data  Base,  consisting 
of  two  parts,  and  is  intended  for  users  and  analysts  as  a  source  of 
network  and  mobility  resource  data. 

b.  While  it  is  intended  that  these  manuals  be  capable  of  use  inde¬ 
pendently,  it  should  be  noted  that  the  Operational  Capability  Descrip¬ 
tion  (OCD)  and  the  Capability  Design  Specification  (CDG)  for  ETNAM 
contain  additional  information  of  value  to  both  the  action  officers  and 
analysts.  The  OCD  contains  a  detailed  statement  of  the  types  of  appli¬ 
cations  for  which  the  ETNAM  model  has  been  designed.  The  CDS  (particu¬ 
larly  Appendices  A  and  b)  contains  a  detailed  and  rigorous  mathematical 
statement  of  the  algorithms  used  in  the  ETNAM  model.  Appendices  A  and 
B  of  the  CDS  are  repeated  in  this  manual  as  Enclosures  1  and  2. 

c.  The  System  Description  has  been  prepared  in  accordance  with  the 
DCA  Instruction  210-175-3  "Documentation  Standards"  for  this  type 
document.  The  order  of  the  major  headings  and  the  information  selected 
for  inclusion  have  been  based  upon  this  instruction.  The  User's  Manual 
(UM)  follows  the  suggested  outline  exactly.  The  suggested  outline  for 
this  Analytical  Manual  (AM)  has  been  modified  slightly  to  permit  the 
presentation  of  all  the  information  about  each  of  the  four  modules  of 
the  system  in  separate  chapters. 

d.  This  Analytical  Manual  is  intended  to  provide  a  basic,  nan- 
mathematical  introduction  to  the  ETNAM  model  and  the  ETNAM  system  of 
computer  programs.  It  provides  a  basic  guide  to  the  preparation  of 
input  data  and  the  operation  of  the  computer  program  from  an  analytical 
point  of  view.  It  discusses  the  methodology  used,  the  adequacy  of 
input  data  for  European  Theater  strategic  mobility  problems,  and  alter¬ 
native  solutions  available.  Each  of  the  four  modules  of  the  computer 
program  is  discussed  in  some  detail.  The  analyst  should  be  able  to 
use  this  information  to  obtain  a  basic  understanding  of  the  model  and 
of  the  required  input  data  and  its  format. 
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e.  It  should  be  noted  at  the  outset  that  the  System  Description 
will  not  be  adequate  in  itself  to  prepare  the  user  for  the  operation 
of  the  ETNAM  model.  The  formulation  of  strategic  mobility  problems 
for  ETNAM  requires  some  basic  familiarity  with  network  terms  and 
methods.  Understanding  of  the  solution  algorithm  of  the  model  requires 
an  understanding  of  some  sophisticated  mathematical  techniques. 

f.  The  basic  formulation  of  the  problem  to  be  solved  and  the  inter¬ 
pretation  of  the  results  should  be  carefully  controlled  Dy  the  user 
himself.  In  the  applications  of  the  ETNAM  model  completed  to  date, 
those  with  full  user  participation  have  been  significantly  more 
successful  than  those  conducted  entirely  by  the  model  developers. 

g.  It  has  also  been  discovered  that  for  many  users  the  quickest 
and  most  successful  way  to  learn  the  uses  of  the  ETNAM  model  is  to  study 
the  formulation  and  results  of  previous  applications.  A  number  of 
these  test  cases  are  now  available  at  the  Research  Analysis  Corporation 
and  in  the  Systems  Analysis  Division  of  the  Office  of  the  Deputy 
Director  for  Logistics  Strategic  Mobility,  J-4,  JCS. 

2.  Problem- Process/Event  Simulation.  The  definition  of  the  problem 
to  be  solved  by  ETNAM  was  developed  after  a  study  of  users '  require¬ 
ments.  A  feasibility  study  was  conducted  to  determine  the  best  method 
available  for  use  in  solving  this  problem  and  a  "breadboard"  version 
of  the  present  ETNAM  model  was  developed  and  tested.  These  develop¬ 
ments  and  others  leading  to  the  present  ETNAM  model  exe  described  in 
detail  In  the  Feasibility  Study,  the  Operational  Capability  Description 
(OCD),  the  Capability  Design  Specification  (CDS)  and  the  Verification 
Test  Plan  (VTP).  Those  analysts  interested  in  more  details  of  the 
development  of  the  system  should  refer  to  these  documents. 

a.  Operating  Objectives.  In  general  terms  the  ETNAM  model  was 
designed  to  provide  a  tool  for  the  detailed  analysis  of  the  capabilities 
of  a  transportation  system  and  to  determine  the  optimum  allocation  of 
resources  in  that  system.  The  model  has  four  objectives  of  operation. 

(1)  In  the  maximum- flow  (max-flow)  objective,  the  model  will 
determine  the  maximum  flow  of  men  and  materials  that  can  be  sustained 
by  a  transportation  network  given  the  origin-destination  pairs,  the 
inventory  of  mobility  resources  available,  the  capacities  of  the  road, 
rail,  air,  waterway,  pipeline  and  other  types  of  links  in  the  trans¬ 
portation  network,  and  the  capacities  of  specified  nodes  of  intersection 
of  three  or  more  links. 

(2)  In  the  minimum-cost  (min-cost)  objective,  the  model  will 
determine  the  set  of  routes  and  the  allocation  of  resources  to  those 
routes  that  will  minimize  cost,  given  the  delivery  requirements,  the 
maximum  number  of  each  type  of  mobility  resource  that  may  be  assigned, 
the  relative  costs  of  each  of  these  resources,  and  the  capacities  of 
each  of  the  links  and  specified  nodes  of  the  transportation  network. 
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(3)  In  the  minimum-time  (min-time)  objective,  the  model  will 
determine  the  set  of  routes  and  the  allocation  of  resources  to  those 
routes  that  will  minimize  flow  unit-hours,  given  the  delivery  require¬ 
ments,  the  maximum  number  of  mobility  resources  that  may  be  assigned, 
the  relative  speeds  of  each  of  the  modes,  and  the  capacities  of  each 
of  the  links  and  specified  nodes  of  the  transportation  network.  The 
flow  units  may  be  all  in  tons  if  all  the  commodity  units  have  been  con¬ 
verted  to  equivalent  tons,  or  they  may  be  in  the  commodity's  natural 
unit.  This  option  will  be  explained  later  in  this  section. 

(4)  Finally,  there  is  a  mixed  minimization  of  cost  and  time 
objective  in  which  th?  mode],  will  determine  the  set  of  routes  and  the 
allocation  of  resources  to  those  routes  which  will  minimize  the  combi¬ 
nation  of  both  cost  and  time  according  to  weights  established  by  the 
user. 


(5)  In  summary,  the  model  may  be  used  in  the  max-flow,  rain¬ 
coat  ,  min-time  or  mixed  minimization  of  cost  and  +ime  operating  objectives. 

b.  Max-Flow.  The  max-flow  objective  may  be  called  a  capabilities 
objective.  That  is,  when  operated  in  this  objective  the  model  deter¬ 
mines  the  maximum  capability  of  the  transportation  system,  given  the 
limitations  of  the  available  resources  and  the  capacity  of  the  network. 

It  should  be  emphasized,  however,  that  this  maximum  flow  is  strictly 
dependent  upon  the  origin- destination  pairs  provided  as  input  data. 

As  a  simple  example,  consider  a  problem  involving  only  two  origin- 
destination  pairs:  one  pair  fifty  airline  miles  apart  and  the  second 
one  hundred  airline  miles  apart.  If  resources  are  limited,  the  model 
will  assign  all  available  resources  to  the  first  pair  and  possibly  none 
to  the  second  pair  because  this  allocation  will  provide  the  maximum 
flow.  Thus  the  maximum- flow  operational  objective  of  the  model  will  not 
be  useful  unless  the  meaning  of  flew  maximization  is  understood.  In 
general,  the  max-flow  objective  is  not  used  except  when  a  general  order- 
of -magnitude  estimate  of  a  network  capability  is  desired. 

c.  Min- Cost.  The  minimum-cost  objective  is  most  often  used  because 
it  provides  an  estimate  of  requirements  for  meeting  a  specific  delivery 
plan.  In  this  operating  objective  the  model  determines  the  minimum- 
cost  set  of  routes  and  mobility  resources  necessary  to  fulfill  delivery 
rt  arements.  It  may  be  used  to  determine  the  feasibility  of  a  movement 
plan  or  to  determine  the  mobility  resource  requirements  to  meet  delivery 
goals. 


d.  Min -Ti:ae.  The  minimum-time  operating  objective  determines  the 
mobility  resource  requirements  and  the  routes  that  will  move  men  and 
materials  through  the  network  to  their  destination  in  the  minimum  time. 
It  is  important  to  understand  that  the  model  minimizes  flow  unit-hours. 
That  is,  it  minimizes  the  sum  of  the  products  of  all  deliveries  and  the 
elapsed  time  from  origin  to  destination.  Users  may  wish  to  examine 
both  the  min-cost  and  the  min-time  solutions  to  the  same  basic  problem. 
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In  general,  a  comparison  of  a  min-time  solution  with  a  min-cost  solution 
will  show  that  men  and  materials  to  be  moved  have  been  shifted  from 
slower  (but  lower  cost)  modes  such  as  waterway  to  faster  (but  higher 
cost)  modes  such  as  air. 

e.  Mixed  Minimization  of  Cost  and  Time.  The  mixed  minimization  of 
cost  and  time  objective  is  not  a  separate  operating  objective  as  far  as 
ETNAM  is  concerned,  but  rather  a  problem  formulation  option  that  the 
analyst  may  use  under  the  min-cost  operating  objective.  Briefly,  the 
method  assigns  time  as  a  resource  to  be  minimized  along  with  the  other 
applicable  costs.  For  each  commodity  -  transportation  mode  pair  a  time 
"productivity"  is  given.  This  is  the  amount  of  time  it  takes  to  move 
one  unit  of  the  specified  commodity  one  mile  on  the  specified  transpor¬ 
tation  mode.  A  price  is  attached  to  the  resource  time.  In  a  min-cost 
run,  time  is  minimized  to  the  extent  that  the  total  cost  is  minimized 
which  depends  on  the  prices  of  the  other  resources.  It  is  not  necessary 
to  follow  this  procedure  for  all  commodities.  It  is  frequently  the  case 
that  some  commodities  should  be  given  time-priorities  while  others 
should  be  shipped  via  the  least  cost  routes.  For  example,  PAX  and  Unit 
Equipment  might  have  elapsed  time  included  as  one  of  the  resources  used 
(and  priced)  while  other  commodities  are  handled  in  the  usual  min-time 
way.  This  approach  would  have  the  effect  of  placing  PAX  and  Unit  Equip¬ 
ment  shipments  on  the  most  expeditious  routes  while  the  other  commodities 
(e.g.,  POL)  would  be  sent  via  the  least  cost  routes  without  the  time 
minimization  consideration. 

f .  Problem  Formulation.  A  problem  for  ETNAM  may  be  formulated  in 
either  of  two  ways.  One  way  is  to  convert  all  flows  and  requirements 
for  commodity  movement  into  some  equivalent  unit  such  as  short  tons. 

This  allows  the  network  capacities  for  the  links  and  nodes  to  be  given 
in  short  tons  per  day.  Another  way  is  to  express  all  commodities  and 
commodity  requirements  in  terms  of  their  natural  units  (e.g.,  numbers 
of  passengers,  barrels  of  gasoline,  tons  of  coal,  etc.).  This  conven¬ 
iently  allows  the  expression  of  network  capacities  in  terms  of  vehicle3- 
per- time-period,  which  is  the  way  in  which  most  traffic  statistics  are 
collected.  The  "equivalent  unit"  formulation  was  the  basis  for  the 
original  ETNAM  development  and  is  described  in  more  detail  in  Enclosures 
1  and  2.  The  "natural  unit"  formulation  has  subsequently  been  incor¬ 
porated  in  the  ETNAM  system  and  is  fully  described  in  Enclosure  3» 

g.  Outputs .  The  outputs  of  ETNAM  provide  considerable  detail  about 
the  transportation  system  under  study.  The  solution  itself  is  presented 
in  several  forms.  The  routes  chosen  and  the  resources  assigned  to  them 
are  provided  by  commodity.  The  movements  by  commodity  over  each  of  the 
links  of  the  network  are  presented.  If  the  network  capacities  are 
expressed  in  some  equivalent  unit,  like  short  tons,  a  unit-mile  or  ton- 
mile  flow  by  mode  is  given.  The  total  flow  through  origin  and  through 
destination  nodes  is  shown.  Total  flow  through  up  to  fifty  intermediate 
nodes  may  be  accumulated  at  the  user's  option.  Summaries  of  commodity 
flow  and  of  resource  utilization  are  also  presented.  Analytical  data 
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is  presented  that  indicates  the  decrease  in  cost  (or  the  increase  in 
flew)  that  could  be  obtained  if  one  or  more  units  of  each  mobility 
resource  were  available.  The  decrease  in  cost  (or  the  increase  in 
flow)  that  would  result  from  one  more  ton  of  capacity  on  each  of  the 
links  or  constrained  nodes  in  the  network  is  also  presented.  In 
linear  programming  these  analytical  data  are  referred  to  by  several 
names  including  shadow  prices,  marginal  costs,  and  dual  solution. 

3-  Discussion  of  Methodology.  The  ETNAM  System  consists  of  an  Input 
Module,  a  RACAT  Module,  an  Output  Module,  and  a  POP  (Post  Optimal  Pack¬ 
age)  Module.  Figure  1  shows  them  schematically.  These  module*  re 
discussed  in  detail  in  Chapters  2,  3,  4  and  5  respectively.  The  3THAM 
system  is  based  upon  a  generalized  algorithm  called  the  Resource  Allo¬ 
cation  and  Chain  Analysis  Technique  (RACAT)  that  was  developed  speci¬ 
fically  for  this  purpose.  The  Input  and  Output  Modules  are  primarily 
data  processing  modules  designed  solely  to  edit  and  format  the  inputs 
and  outputs  for  the  RACAT  Module.  The  POP  Module  gives  upper  and  lower 
limits  to  which  the  right-hand- sides  (i.e.,  link  capacities,  resource 
inventories,  node  capacities  and  requirements)  may  be  changed,  one  at 
a  time,  without  requiring  a  change  in  the  current  solution  or  basis. 

If  this  information  is  desii id  then  POP  is  employed  after  a  solution 
from  RACAT  is  obtained.  POP  is  fully  described  in  Enclosure  4-. 

a.  Input  Module.  All  inputs  to  the  ETNAM  system  are  prepared  man¬ 
ually.  That  is,  there  is  no  automated  system  for  the  preparation  of 

the  input  data  involving  the  selection  of  existing  data  from  the  National 
Military  Command  and  Control  System  (NMCCS)  data  banks.  A  typical  run 
of  the  ETNAM  system  requires  less  than  one  thousand  punched  cards.  Of 
these,  approximately  two  thirds  are  network  data  cards  that  are  based 
upon  theater  data  and  hence  apply  to  all  problems  for  that  theater. 

The  Input  Module  is  thus  a  relatively  simple  computer  program  that 
edits  and  formats  the  input  data  and  prepares  reports  that  can  be 
checked  by  the  user.  Complete  descriptions  of  these  reports  are  con¬ 
tained  in  the  User's  Manual. 

b.  Output  Module.  This  module  takes  the  output  of  the  RACAT 
Module  and  edits  and  format  it  for  presentation  to  the  user.  The  raw 
output  of  the  RACAT  Module  is  the  mathematical  solution  to  the  problem 
and  is  not  easily  understood  or  interpreted.  Some  special  problems  must 
be  solved  in  organizing  and  editing  this  data  fer  presentation,  but 
standard  data  processing  techniques  are  used.  In  general,  the  Output 
Module  reports  have  been  found  to  be  useful  and  complete  in  the  appli¬ 
cations  processed  to  date.  These  reports  are  also  described  in  the 
User's  Manual. 

c.  RACAT  Module.  This  module  contains  the  mathematical  algorithm 
and  is  the  core  of  the  system.  It  is  an  extension  to  the  multi-commodity 
network  flow  algorithm  suggested  by  Ford  and  Fulkerson  in  1958  and 
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FIGURE  1.  ETNAM  SYSTEM  FLOW  SCHEMATIC 
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extended  to  the  minimum- cost  formulation  by  Tomlin  in  19 66  .  This 
algorithm  finds  the  optimum  solution  to  network  flow  problems  where 
there  are  distinct  origin-destination  pairs  to  be  connected  by  routes 
in  such  a  way  that  all  requirements  are  satisfied  and  none  of  the 
individual  link  and  node  capacities  nor  the  resource  inventories  are 
exceeded.  There  are  other  models  that  will  find  solutions  to  the 
problem,  but  these  are  not  optimal  solutions  and  there  is  no  way  to 
determine  how  far  these  solutions  are  from  the  optimum. 

(1)  A  simple  example  of  a  multi- commodity  problem  is  the  pro¬ 
blem  of  the  telephone  company  on  Mother's  Day.  Thousands  of  sons  are 
calling  thousands  of  mothers  and  the  telephone  company  must  connect 
each  son  to  the  corresponding  mother  in  such  a  way  that  none  of  their 
telephone  lines  is  overloaded.  There  are  excellent  single- commodity 
models  that  will  find  the  least  cost  solution  to  single- commodity 
problems.  If  the  single-commodity  model  were  used  to  solve  the  tele¬ 
phone  company's  problem  outlined  above,  it  wouJ d  in  effect  connect 
each  son  to  the  nearest  mother  --  an  unsatisfactory  solution  for  most 
sons  and  mothers. 

(2)  A  number  of  multi- commodity  computer  models  have  been 
written  using  the  basic  Ford  and  Fulkerson  method.  The  Logistics 
Research  Project  of  the  George  Washington  University,  the  Civil  En¬ 
gineering  Department  of  the  Johns  Hopkins  University,  the  University  of 
California  at  Berkley,  and  others  have  such  models.  The  RACAT  algorithm 
is  different  from  these  in  that  it  includes  node  capacity  and  resource 
constraints  as  well  as  link  capacity  constraints. 

(3)  Mobility  resources  are  likely  to  be  scarce  in  most  military 
theaters.  In  the  European  Theater  the  transportation  network  itself  is 
one  of  the  most  dense  and  most  sophisticated  in  the  world.  The  capa¬ 
cities  of  the  roads,  railroads,  waterways,  pipelines,  and  airways  -- 
while  they  may  constrain  individual  movements  —  are  large,  and  there 
are  sufficient  alternatives  that  almost  any  delivery  requirements  can 

be  met.  The  problem  is  to  simultaneously  select  the  routes  and  allo¬ 
cate  resources  to  those  routes  in  an  optimal  fashion.  The  RACAT  algor¬ 
ithm  is  an  optimizing  algorithm  that  simultaneously  selects  routes  and 
allocates  resources  while  maximizing  flow,  minimizing  cost,  or  minimizing 
time. 


(4)  The  RACAT  algorithm  (and  the  original  Ford-Fulkerson  algor¬ 
ithm)  is  a  combination  of  the  simplex  algorithm  of  linear  programming 
and  the  shortest  chain  algorithm  of  graph  theory.  In  highly  simplified 
terms  the  simplex  algorithm  can  be  said  to  consist  of  4  steps  as  follows : 


L.  R.‘  Ford  and  D.  R.  Fulkerson,  "A  Suggested  Computation  for  Maximal 
Multi-Commodity  Network  Flows",  Mgmt .  Sc . ,  October,  1958;  J.  A.  Tomlin, 
"Minimum-Cost  Multicommodity  Network  Flows",  ORSA,  December,  1966. 
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(a)  Find  a  feasible  solution  consisting  of  a  set  of  routes 
and  the  corresponding  resources  that  will  meet  the  requirements. 

(b)  Select  a  route  not  now  used  that  will  improve  that 

solution. 

(c)  Update  the  solution.  This  may  mean  dropping  a  route 
in  the  previous  solution. 

(d)  Continue  repeating  steps  (b)  and  (c)  until  there  is 
no  available  route  that  will  improve  the  solution.  When  there  is  no 
such  route,  an  optimum  solution  has  been  found. 

(5)  In  the  RACAT  version  of  the  simplex  algorithm,  the  shortest 
chain  algorithm  is  used  to  find  the  improving  routes.  That  is,  the 
shortest  chain  algorithm  is  substituted  for  step  (b). 

(6)  In  order  to  fully  understand  the  solution  procedure  and 
the  computer  code  of  the  RACAT  Module,  it  is  necessary  to  have  a  back¬ 
ground  knowledge  of  the  simplex  algorithm  and  certain  aspects  of  graph 
theory.  While  this  background  knowledge  is  not  necessary  to  use  the 
model  or  understand  its  results,  such  knowledge  will  be  necessary  to 
anyone  who  plans  to  maintain  or  modify  the  computer  program.  A  biblio¬ 
graphy  of  books  and  articles  on  the  subject  is  included  as  Enclosure  6. 
Enclosures  1,  2  and  3  are  technical  papers  on  the  RACAT  algorithm  and 
contain  a  rigorous  mathematical  treatment  of  the  algorithm.  Analysts 
who  wish  to  understand  the  procedure  in  all  its  detail  should  refer  to 
these  technical  papers.  Sample  data  listings  appear  in  Enclosure  5. 

d.  POP  Module.  This  module  is  a  post  optimal  procedure  giving  the 
analyst  data  sensitivity  information.  The  sensitivity  analysis  that  is 
done  is  right-hand-side-ranging.  The  right-hand-sides  referred  to  are 
the  constraint  limits  of  the  problem,  thrt  is,  network  capacity  and 
resource  inventory  constraints  and  movement  requirements.  During  a  POP 
computer  run  any  number  of  right-hand-sides  may  be  analyzed  one  at  a 
time.  That  is,  while  the  "range"  of  a  given  right -hand- side  is  found, 
the  others  are  held  fixed. 

(1)  The  range  of  a  right-hand-side  is  the  set  of  values  the 
right-hand-side  may  assume,  keeping  the  other  right-hand-side  values 
fixed,  without  requiring  a  change  in  the  solution.  This  information  may 
give  insight  to  the  analyst  regarding  binding  or  potentially  binding 
constraints  for  that  solution.  It  must  be  remembered  that  this  is  a 
local  sensitivity  analysis,  that  is,  a  change  in  the  problem  could  give 
quite  different  sensitivity  results. 

(2)  Although  POP  makes  use  of  the  properties  of  the  simplex  pro¬ 
cedure,  it  is  net  necessary  to  understand  this  procedure  co  be  able  to 
read  the  output  reports.  However,  on >  should  be  aquainted  with  the 
mathematical  techniques  used  to  understand  the  computer  code.  Technical 
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information  may  be  found  in  the  linear  programming  references  in  the 
bibliography  in  Enclosure  6.  Enclosure  4  gives  a  rigorous  mathematical 
treatment  of  the  solution  procedure  employed  by  POP. 

4.  Adequacy  of  Available  Data. 

a.  In  general,  three  types  of  data  are  needed  to  formulate  a  pro¬ 
blem  for  the  ETNAM  system.  These  are  network  data,  resource  data,  and 
control  parameters. 

(1)  The  network  data  tend  to  be  fixed  for  all  problems  involv¬ 
ing  a  given  theater.  New  highway  and  railway  links  are  constructed  and 
old  ones  expanded,  but  these  changes  are  infrequent.  Thus,  the  same 
set  of  network  punched  cards  can  normally  be  used  to  run  a  whole  series 
of  runs  in  a  given  theater. 

(2)  Basic  resource  data  change  very  little.  New  vehicles  and 
trar,  sportation  techniques  are  developed,  but  these  changes  are  infrequent. 
Nevertheless  resource  productivities  and  resource  prices  (ten-year 
systems  costs)  are  normally  developed  carefully  for  each  run  or  series 

of  runs.  The  reason  is  that  resource  productivities  are  a  function  of 
the  operating  procedures  that  are  planned  for  the  vehicles  as  well  as 
their  physical  characteristics.  The  productivity  of  a  carrier  will  be 
quite  different  if  a  twenty-hour  operating  day  is  planned  than  if  an 
eighteen-hour  operating  day  is  planned.  The  relative  prices  to  be  used 
will  also  vary  from  study  to  study  depending  upon  the  ground  rules 
agreed  upon.  It  is  anticipated,  therefore,  that  resource  data  will  be 
updated  for  each  study. 

(3)  The  control  parameters  change  from  run  to  run  even  within 
the  same  study  and  must,  therefore,  be  prepared  each  time.  The  control 
parameters  specify  the  operating  objective  of  the  individual  run,  the 
number  of  iterations  between  inversion  and  other  parameters  having  to 
do  with  the  control  of  the  program  during  the  solution  of  the  problem. 
These  specifications  are  reasonably  straight-forward  and  the  development 
of  this  data  is  not  difficult. 

b.  There  are  three  major  sources  of  data  of  the  first  two  types 
that  have  been  regularly  used  to  date.  The  Defense  Intelligence  Agency 
(DIA)  has  available  extensive  network  data  for  many  theaters  throughout 
the  world.  Some  users  prefer  to  make  their  own  estimates  of  capacity 
based  on  DIA  base  data.  Additional  information  on  the  European  Theater 
is  available  in  the  DCA  report  (prepared  by  RAC)  on  the  European  Theater 
Transportation  System  (ETTS).  Volume  IV  of  the  SD,  the  ETNAM  Data  Base 
(also  prepared  by  RAC;,  contains  related  information  on  Europe,  South¬ 
east  Asia  and  Northeast  Asia. 

c.  The  Research  Analysis  Corporation  has  on  file  a  European  network 
in  proper  punched  card  form.  These  network  data  are  based  on  DIA  data 
modified  as  a  result  of  the  ETTS  report.  It  has  been  used  for  a  number 
of  applications. 
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d.  In  general  terms,  the  available  data  for  the  operation  of  the 
ETNAM  model  has  been  found  to  be  adequate  for  the  European  Theater. 

This  does  not  mean  that  there  is  not  room  for  improvement.  The  ETNAM 
model  is  very  sensitive  to  resource  productivities  and  to  resource 
prices.  Very  small  changes  in  either  of  these  categories  of  data  may 
cause  significant  changes  in  the  solution.  For  tnis  reason,  systems 
costs  and  productivities  should  be  carefully  reviewed  by  the  user. 

e,  Based  on  experience  in  the  application  of  the  ETNAM  model  to 
date,  it  is  strongly  recommended  that  the  user  review  and  study 
resource  prices  and  productivities  prior  to  the  first  run  of  the  model. 
Such  a  review  insures  that  the  user  is  satisfied  with  these  data  and, 

in  addition,  seems  to  assist  the  user  in  understanding  the  basic  approach 
of  the  ETNAM  model.  While  there  is  no  doubt  that  there  are  adequate 
data  available  for  the  operation  of  the  system,  only  the  user  can  deter¬ 
mine  if  these  data  are  satisfactory  for  his  purposes. 

5.  Alternative  Solutions.  There  are  two  senses  in  which  alternative 
solutions  to  problems  structured  for  the  ETNAM  model  may  be  considered. 
There  are  alternative  solutions  in  a  mathematical  sense  and  there  are 
alternative  solutions  in  an  operational  sense. 

a.  The  solution  generated  by  the  ETNAM  model  is  an  optimum  solution 
to  the  network  flow  problem  presented  to  it.  Mathematically  there  are 
only  two  possible  cases:  one,  there  is  one  and  only  one  optimum  solu¬ 
tion;  and  two,  there  is  an  infinite  number  of  optimum  solutions.  In 
operational  terms,  the  problem  may  have  only  one  set  of  routes  and  one 
set  of  flows  on  those  routes  that  is  optimum.  That  is,  all  other 
possible  sets  cost  more  (or  take  more  time,  or  achieve  less  flow). 

b.  Maximum  flow  problems  frequently  have  alternative  solutions  in 
the  mathematical  sense.  That  is,  there  may  be  another  route  or  routes 
that  could  be  substituted  and  the  same  (maximum)  flow  would  be  achieved. 
If  this  is  the  case  then  there  is  actually  an  infinite  number  of  solu¬ 
tions.  If  we  have  two  solutions  we  can  have  an  infinite  number  of 
linear  combinations  of  the  two.  Operationally  this  is  the  same  as  say¬ 
ing  that  if  we  can  add  another  route  we  can  have  an  infinite  number  of 
variations  by  simply  changing  the  levels  of  flow  on  the  routes  in  our 
solution. 

c.  While  there  will  be  cases  where  the  user  is  interested  In  alter¬ 
native  solutions  in  the  mathematical  sense,  he  will  more  often  be  inter¬ 
ested  in  operational  alternatives.  That  is,  the  military  analyst  may 
prefer  an  alternate  solution  that  is  not  a  mathematical  optimum  but  has 
characteristics  that  have  advantages  not  directly  representable  in  the 
mathematical  formula -ion  of  the  problem. 

d.  A  simple  example  may  serve  to  illustrate  this  problem.  The 
mathematical  optimum  solution  to  a  problem  may  call  for  the  heavy  use 
of  a  link  in  the  network  that  is  considered  to  be  vulnerable  to  enemy 
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action  by  the  military  analyst.  Thus  the  military  analyst  may  wish  to 
reduce  the  flow  over  that  link  even  though  he  knows  that  such  a  reduc¬ 
tion  will  increase  the  cost  of  the  overall  movement. 

e.  Operational  alternative  solutions  are  readily  obtainable  using 
the  ETNAM  model.  They  do,  however,  involve  re-runs  of  the  model  with 
the  new  restrictions  incorporated  in  the  system.  In  the  example  cited 
above,  the  capacity  of  the  heavily  used  link  would  be  reduced  and  the 
problem  re-run.  ETNAM  has  a  re-start  feature  which  permits  a  resumption 
of  the  algorithm  at  a  particular  point  in  the  solution  process.  Thus, 
if  an  alternative  solution  is  desired  it  is  not  necessary  to  "start 
from  scratch." 


O 
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CHAPTER  2.  INRJT  MODULE 


1.  Purpose  and  Scope.  The  primary  purpose  of  the  Input  Module  is  to 
read  data  and  to  perform  fundamental  editing  functions  in  structuring 
the  data  for  the  RACAT  Module  algorithms.  This  chapter  describes  the 
inputs  required  and  outputs  generated  by  the  Input  Module  and  discusses 
the  methods  used  in  the  editing  processes. 

2.  Inputs  and  Outputs. 

a.  Input  Data  Requirements.  Two  groups  of  input  data  are  required 
to  operate  the  ETNAM  system.  The  first  set  must  be  in  fixed  sequence 
while  the  other  may  be  entered  in  any  order.  Figure  2  is  a  schematic 
description  of  the  input  deck  arrangement  for  the  Input  Module.  This 
section  describes  the  general  formats  for  each  type  of  input  data. 
Enclosure  5  has  some  sample  data  listings. 

(l)  Sequenced  Data.  The  Input  Module  has  been  designed  to 
operate  primarily  upon  input  data  which  are  not  required  to  be  in  any 
fixed  sequence  for  entry.  The  exceptions  to  this  scheme  which  must  be 
in  sequence  are : 

Title  Information 

Run  Parameters 

Output  Titles  (in  any  order  within  this  group) 

Network  Data  (in  any  order  within  this  group) 

Node  Constraint  Data  (if  desired) 

Each  of  these  is  discussed  below. 

(a)  Title  Information.  The  first  data  card  is  used  as  a 
title  entry.  All  80  columns  are  available  and  any  alphanumeric  entry 
is  permitted.  This  entry  becomes  a  header  for  each  page  of  the  reports 
in  the  Input,  Output  and  POP  Modules. 

(b)  Run  Parameters.  Entries  for  this  card,  which  must 
follow  the  Title  card,  are  as  follows : 

Columns  Information 

1-8  Run  objective  (MAX  FLOW,  MIN  COST, 

or  MIN  TIME) 

15  Number  of  substitution  groups 

allowed  (s3) 

19-20  Number  of  transportation  modes  (s20) 

24-25  Number  of  commodities  (s20) 

29-30  Number  of  resources  (*50) 

3^-35  Number  of  node  constraints  (£50) 

38-U0  Type  of  link  capacities 

0  =  some  common  unit  of  flow  (e.g., 
short  tons) 


2-1 


m 


INyMWWWWIWI  W1  f¥U4 - 


'*> '  •*i',r ! "  'S*i*  '■  f; 


1  =  vehicles -per- time-period  (V/TP) 

43-45  Number  of  iterations  before  inversion 
(see  RACAT  MODME,  Chapter  4) 

47-50  Maximum  number  of  iterations  permitted 
(see  RACAT  MQDUIE,  Chapter  4) 

(c)  Output  Titles.  These  cards  must  follow  the  run  para¬ 
meter  card  but  may” be  entered  in  any  order.  The  4  characters  in  columns 
1-4  identify  the  type  of  titles.  Since  the  titles  are  assigned  to 
modes,  resources,  and  commodities  by  number,  they  should  be  entered 
in  order  on  the  card.  The  types  of  entry  and  dimensions  are  a3  follows: 


Columns  1-4  Type  of  Titles 


Dimension 


NAMO  Names  of  modes  £20 
NACO  Names  of  commodities  £20 
HARE  Names  of  resources  £50 


and  the  format  for  all  types  is: 


Columns  Information 


1-4 

10-13 

16-19 

22-25 

28-31 

34-37  > 

40-43 

46-49 

52-55 

58-61 

64-67 


Card  type 


10  titles  per  card  of  4  characters 
each  (left  adjusted) 


(d)  Network  Data. 


1.  The  first  card  of  this  Bet  must  be  a  title  card 
with  "ARCS"  in  columns  1-4.  The  network  data  cards  follow  and  a  blank 
card  at  the  end  signifies  that  all  network  data  cards  have  been  entered. 
It  will  be  noted  here  that  throughout  the  SD,  the  word  arc  is  used  inter 
changeably  with  the  word  link  except  where  the  specific  word  "ARC"  or 
"ARCS"  is  vised  as  input  data.  The  format  for  data  cards  is  as  follows: 


Columns  Information 


1-6 

7-12 

13-18 

19-24 

25-32 


FROM  node,  alphanumeric,  left  adjusted 
TO  node,  alphanumeric,  left  adjusted 
Mode  of  arc,  integer 
Condition  of  arc,  real  number 
Distance  of  arc,  real  number 
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33-40  Capacity  of  arc,  real  number 
50  Direction  of  arc  (0  =  undirected; 

1  =  directed) 

51-70  Commodity  availability  (l  column  for 
each  commodity.  0  =  will  not  handle; 
1  =  will  handle) 


2.  The  network  data  cards  may  be  entered  in  any 
order  since  they  are  sorted  in  the  Input  Module.  It  will  be  noted 
again  that  the  analyst  has  the  option  of  stating  the  capacities  of 
the  links  in  terms  of  tons-per-time-period  or  in  terms  of  vehicles- 
per-time-period.  Whichever  choice  is  made  must  be  used  consistently 
throughout  the  network. 


(e)  Node  Constraint  Data.  If  the  capacity  (either  tons 
or  V/TP)  of  certain  nodes  is  to  be  constrained,  the  pertinent  data  must 


appear  after  the  network  listing, 
arrays : 

NODCQN  (I) 
NOEMOD  (I) 

RHS(NARC  +  NKES 

+  I) 

Where  NARC 
NRES 


These  data  are  stored  in  the  following 


=  Name  of  the  constrained  node  I 
=  Number  of  the  mode  of  con¬ 
strained  node  I 

=  Throughput  limit  of  constrained 
node  I 

=  The  number  of  arcs 
=  The  number  of  resources 


The  throughput  limits  become  right -hand- side  entries  for  the  constrained 
nodes.  An  error  message  is  generated  if  more  than  50  nodes  are  entered. 
The  nodes  which  are  origins  or  destinations  are  not  allowed  to  be  con¬ 
strained.  The  format  for  these  data  cards  is : 


Columns  Information 

I- 6  Name  of  node  (left  adjusted) 

7-10  Number  of  mode  (right  adjusted) 

II- 20  Throughput  limit  (node  capacity 

constraint ) 

(2)  Unsequenced  Data.  The  remaining  input  data  can  be  entered 
in  any  sequence  since  each  card  specifies  in  columns  1-4  the  name  of 
the  type  of  data  to  be  entered.  Methods  for  entering  each  type  of  data 
are  discussed  below.  Each  card  type  is  specified  by  a  4  character 
mnemonic.  The  general  format  for  these  data  cards  except  origin/ 
destination  pairs  and  tolls  is : 


Columns 


Information 


6-10  Integer  value 

11-15  Integer  value 

16-20  Integer  value 

21-30  Real  value 

31-35  Real  value 

36-40  Integer  value 

41-45  Integer  value 

46-50  Integer  value 

54-60  Real  value 

61-65  Real  value 

The  values  in  card  columns  31-35  and  61-65  are  used  only  for  resource 
data  and  is  discussed  in  that  section. 

(a)  PROD  or  PAID  -  Resource  Productivities. 

1.  Either  PROD  or  PALD  (but  not  both)  may  be  used 
for  any  trio  of  commodity,  mode,  resource  productivity.  These  data 
are  stored  in  the  arrays  MASRES  and  RESMAS  as  follows: 

MASHES  (I,  l)  =  Number  of  commodity 
MASRES  1 I,  2)  -  Number  of  mode 
MASRES  (I.  3)  =  Number  of  resource  type 
RESMAS  (I)  =  Productivity  for  this  type 

where  I  is  the  index  for  the  non- zero  productivities. 

2.  There  are  two  methods  of  inputing  the  producti¬ 
vities  which  are  entered  in  their  natural  units  (e.g.,  fraction  of  a 
bus  per  person  per  mile  or  fraction  of  a  truck  per  ton  per  mile): 

a.  When  the  productivity  is  supplied  by  the 
user  the  PROD  card  will  be  used. 

b.  When  the  user  supplies  the  payload  of  the 
resource  and  the  total  mileage  traveled  per  day  (including  a  return 
trip  empty)  for  a  typical  unit  of  that  resource,  the  model  will  cal¬ 
culate  the  productivity.  This  will  require  use  of  the  PALD  card. 

When  the  distance  is  greater  than  1  the  model  will  compute  the  inverse 
of  the  product  of  the  payload  and  one-half  the  distance.  If  the  dis¬ 
tance  is  1  (i.e.,  transfer  units)  the  productivity  is  the  inverse  of 
the  payload. 

3.  The  productivities  are  entered  sequentially  into 
the  matrix  in  any  order.  If  no  value  is  entered  for  a  given  trio  of 
commodity,  mode  and  resource,  that  productivity  is  assumed  to  be  zero. 

4.  A  choice  of  resources  may  be  given  to  a  commodity 
mode  pair.  This  is"  known  as  substitution  of  resources.  ETNAM  can 
accommodate  up  to  three  substitution  groups.  If  a  commodity-mode  pair 


> 

>  First  set 

>  Second  set  (if  desired) 
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is  to  have  a  choice  of  up  to  3  sets  of  resources,  each  set  is  put  into 
a  different  substitution  group.  But  if  another  commodity-mode  pair  has 
no  competitive  resources,  then  each  substitution  group  has  the  identi¬ 
cal  resource  sets  for  that  commodity-mode  pair  and  must  be  repeated  for 
each  substitution  group  in  the  input  data. 

£.  For  purposes  of  inputing  the  data,  the  substitu¬ 
tion  group  into  which  a  resource  productivity  belongs  needs  to  be 
identified.  At  most  50  resources  may  be  used.  These  resources  are 
numbered  sequentially  from  1  to  at  most  50  in  substitution  group  1. 

The  same  resources  are  numbered  from  51-100  and  101-150  in  substitution 
groups  2  and  3  respectively. 

6.  Error  messages  are  produced  if  any  subscript  is 
blank  or  zero  or  if  more  than  700  entires  are  attempted.  The  format 
for  these  data  cards  is: 


Columns 

Information 

1-4 

PROD 

PALD 

6-10 

Number  of  commodity 

Number  of  commodity 

11-15 

Number  of  mode 

Number  of  mode 

16-20 

Number  of  resource 

Number  of  resource 

21-30 

Productivity  value 

Payload 

31-35 

Blank 

Total  mileage 
(includes  empty 
return  trip) 

(a  second  set  if  desired  according  to  standard 

format ) . 

(b)  INVE  -  Resource  Inventories.  These  data  are  stored 
in  the  array  RHS.  These  data  become  right-hand- side  entires  for  the 
resources.  A  subscript  is  required  to  identify  the  number  of  resource 
type.  An  error  message  is  generated  if  the  subscript  is  blank  or  zero 
or  if  more  than  50  resource  types  are  entered.  The  format  for  the  INVE 
card  type  i s : 

Columns  Information 

1-4  INVE 

16-20  Number  of  resource 

21-30  Total  in  inventory 

(a  secOiid  set  if  desired  -  columns  46-50  and  51-60). 


1.  These  data  are  stored  in  NORMS  (I,  4),  where 
I  -  number  of  o/D  pair.  These  data  are  origin/destination  pairings  by 


commodity  number  and  a  delivery  requirement  for  that  pair  if  appropriate 
The  delivery  requirement  is  entered  in  the  natural  unit  of  the  commodity 
Normally  a  delivery  requirement  is  imposed  on  an  o/D  pair.  The  alter¬ 
native  to  that  is  to  leave  some  or  all  o/D  pairs  for  a  given  commodity 
"open"  (i.e.,  do  not  specify  a  requirement)  but  specify  a  requirement 
for  that  commodity  for  the  "open"  o/D  pairs  in  a  DELI  card  (described 
below).  The  RACAT  algorithm  will  assign  appropriate  flows  to  the 
"open"  O/D  pairs,  the  total  of  which  is  equal  to  the  requirement 
specified  in  the  appropriate  DELI  card. 

2.  If  an  O/D  pair  has  a  commodity  delivery  require¬ 
ment  of  zero  but  18"  not  to  be  an  "open"  requirement,  that  is,  it  is  to 
have  a  requirement  in  a  subsequent  run,  then  the  variable  ID  is  set  to 
1.  This  effectively  reserves  a  requirement  row  in  the  mathematical 
structure  of  the  problem  for  subsequent  use.  If  the  pair  is  to  be 
"open",  ID  is  left  blank  and  a  DELI  card  must  be  present. 

j}.  An  error  message  is  generated  if  more  than  150 
o/D  pairs  are  entered  or  if  any  commodity  number  is  >20. 

4.  The  first  card  of  this  set  must  be  a  title  card 
with  "ODPR"  in  columns  1-4.  The  format  of  the  data  cards  following 
it  is : 

Columns  Information 

1-4  ODPR 

7-10  Number  of  commodity,  integer 

14-19  Origin  node,  left  adjusted 

20-25  Destination  node,  left  adjusted 

-6-35  Delivery  requirement  for  this  pair  if 
any,  a  real  number 

36  ID  =  0,  if  delivery  requirement  is  present 

=  0,  if  delivery  requirement  is  zero 

but  the  O/D  pair  is  "open". 

=  1,  if  delivery  requirement  is  zero 
but  the  O/D  pair  is  not  "open". 

(a  second  set  if  desired  in  columns  37-40,  44-49, 

50-55,  56-65,  66). 

The  last  card  must  have  "ODPR"  in  columns  1-4  and  "999"  in  columns 
8-10.  All  integers  in  the  above  format  are  right  adjusted  in  their 
respective  fields. 

(d)  DELI  -  Delivery  Requirements.  These  data  are  stored 
in  the  array  DELREQ  (i),  where  I  =  number  of  commodity.  Any  O/D  pairs 
left  "open"  need  a  total  delivery  requirement  by  commodity.  Closed 
O/D  pairs  do  not  have  corresponding  DELI  cards.  An  error  message  is 
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produced  if  a  commodity  number  is  blank,  zero,  or  >20.  The  format  for  the 
DELI  card  type  is: 

Columns  Information 

1-4  DELI 

16-20  Number  of  commodity 

21-30  Total  delivery  requirement  for  this 

commodity. 

(a  second  set  if  desired  -  columns  46-50  and  51-60). 

(e)  PHIS  -  Prices  on  Resources.  These  data  are  stored  in 
the  array  PHI  (i),  where  I  =  number  of  resource  type. 

1.  Prices  on  resources  are  required  to  operate  the 
model  in  a  min-cost  objective.  Resource  prices  must  not  be  used  for 
the  max-flow  or  min- time  operating  objectives.  An  error  message  is 
generated  if  the  resource  type  number  is  blank,  zero,  or  greater  than 

50. 


'  2.  Prices  are  entered  in  their  actual  dollar  cost 
form.  The  Input  Module  normalizes  these  prices  by  searching  for  the 
minimum  non- zero  price  and  then  converting  all  the  actual  prices  to 
a  relative  scale.  The  format  for  the  PHIS  card  type  is: 

Columns  Information 

1-4  PHIS 

16-20  Number  of  resource 

21- 30  Price 

(a  second  set  if  desired  -  columns  46-50  and  51-60). 

(f )  TOLL  -  Tolls  on  Certain  Arcs. 

1.  These  data  are  stored  in  the  array  TOLLS.  Tolls 
may  be  placed  on  selected  arcs  by  specifying  three  items:  mode;  FROM 
node; and  TO  node.  The  program  will  loca  ,e  the  desired  arc  row  and 
impose  the  specified  toll  on  the  particular  arc.  An  error  message  is 
produced  if  the  specified  arc  cannot  be  located  in  the  arc  matrix. 

2.  The  format  for  the  TOLL  card  is  similar  to  the 
0DPR  cards.  The  first,  card  of  this  set  must  be  a  title  card  with  "TOIL" 
in  column  1-4.  The  format  of  the  data  cards  following  it  is: 

Columns  Information 


1-4 

7-10 


TOLL 

Number  of  arc  mode,  integer-right  adjusted 


Columns 


Information 


14-19  Name  of  FROM  node,  left  adjusted 

20- 25  Name  of  TO  ncde ,  left  adjusted 

26-35  TOLL  to  be  Imposed,  a  real  number 

(a  second  set  if  desired  in  columns  37-40,  44-49, 

50-55,  56-65). 

The  last  card  must  have  "TOLL"  in  columns  1-4  and  "999"  in  columns  8-10. 

(g)  RATE  -  Rates  of  Speed.  These  data  are  stored  in  the 

array  RATES  (i),  where  I  =  number  of  mode.  Rates  of  speed  by  mode  may 

be  entered  in  order  to  calculate  delay  times  on  the  arcs.  If  the  model 
is  to  be  operated  in  a  minimum  time  objective,  these  data  are  required. 
An  error  message  is  generated  if  the  number  of  mode  is  blank,  zero,  or 
greater  than  20.  The  format  for  the  RATE  card  type  is : 

Columns  Information 

1-4  RATE 

16-20  Number  of  mode,  integer-right  adjusted 

21- 30  Speed,  a  real  number 

(a  second  set  if  desired  -  columns  46-50  and  51-60). 

(h)  ARCC  -  Arc  Data  Update.  This  feature  is  discussed 
under  the  section  "Methods"  on  performing  network  data  updates.  The 
format  for  the  ARCC  card  type  is : 


(i) 


Columns 

Information 

1-4 

ARCC 

11-15 

Data  field 

16-20 

Constant  to  be  added 

21-30 

Multiplier 

36-40 

Limit  value 

41-45 

Data  field 

46-50 

Limit  value 

CONF  - 

Conversion  Factors. 

1.  These  data  are  stored  in  the  array  CONF  (i), 
where  I  =  commodity  number.  They  must  not  be  used  when  node  and  link 
capacities  are  expressed  in  vehicles-per-tirae-period. 


2.  When  the  link  and  node  capacities  are  expressed 
in  some  equivalent  unit  like  short  tons-per-time- period,  then  conversion 
factors  are  specified  for  any  commodity  not  expressed  in  the  same  units 
as  the  capacities.  Normally,  the  capacities  are  specified  in  short 
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tons-per-day;  thus  any  commodity  such  as  passengers  or  POL  will  need  to 
be  converted  to  short  tons  to  be  moved  over  the  network.  Both  the 
delivery  requirements  and  the  productivities  are  converted  by  the 
Input  Module. 

The  format  for  the  CQNF  card  type  is : 

Columns  Information 

1-4  CQNF 

16-20  Number  of  commodity 

21-30  Conversion  factor 

(a  second  set  if  desired  in  columns  46-50,  51-6o). 

(j)  TRIP  -  Payload  factor. 

1.  If  the  link  capacities  are  stated  in  vehicles- 
per- time -period  these  data  cards  must  be  present.  The  data  are  stored 
in  the  array  TRIP  (I,  J,  K),  where: 

I  number  of  substitution  group 

J  =  number  of  commodity 

K  =  number  of  mode 

2.  The  factors  are  the  normal  payloads  of  a  vehicle 

(resource)  by  commodity  by  mode  for  each  substitution  group.  If  sub¬ 
stitution  groups  are  used,  the  factors  for  substitution  groups  2  and 

3  need  not  be  entered  unless  they  differ  from  substitution  group  1. 

If  a  certain  commodity  is  not  allowed  on  a  mode,  i.e.,  passengers  on  a 
pipeline,  no  factor  is  entered,  but  for  all  allowable  ccmmodity/mode 
combinations  there  must  be  a  factor. 


An  error  message  is  generated  if  any  commodity  or 
mode  number  is  blank,  zero,  or  320.  The  format  for  the  TRIP  card  type 
is : 

Columns  Information 

I- 4  TRIP 

6-10  Number  of  substitution  group 

II- 15  Number  of  commodity 

16-20  Number  of  mode 

21- 30  Trip  payload  factor 

(a  second  set  if  desired  according  to  standard  format). 

(k)  END--  Signifies  End  cf  Input  Data.  This  must  be  the 
final  card  in  the  input  data  file.  (Columns  1-3  =  END).  An  error  message 

f 

V' 
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is  produced  if  a  data  card  does  not  contain  one  of  the  above  names  for 
the  type  of  input  in  columns  1-4.  Figure  2  is  a  schematic  description 
of  the  input  deck  arrangement  for  the  input  module. 

li.  Output  Reports.  The  input  data  reports  of  the  Input  Module  are 
generated  by  the  WRITES  routine.  The  formats  used  in  these  reports  are 
described  as  follows :* 

(1)  Title  and  Run  Parameters  Report. 

(2)  Rates  of  Speed  by  Mode  Report.  Identifies  the  modes  and 
their  rates  of  speed.  Up  to  20  modes  are  permitted  and  will  be  reported. 

(3)  Commodity  Identification  Report.  It  identifies  each 
commodity  and  gives  its  conversion  factor  if  any.  Up  to  20  commodities 
are  permitted. 

(4)  Network  Data  Listing.  This  report  is  produced  after  the 
arcs  are  sorted  on  the  FROM  node.  The  "SYM  and  Commodities"  is  a 
packed-word  output  in  octal  format  on  the  CDC  6400  and  in  hexa-decimal 
format  for  the  IM  System/ 360. 

(5)  Table  of  Unlq  Nodes . 

(6)  Node  Constraint  Report.  Up  to  50  node  constraints  are 
permitted.  The  report  will  appear  only  if  node  constraints  are  used. 

(7)  Resource  Inventories.  Up  to  50  resource  types  are  per¬ 
mitted  and  will  be  reported. 

(8)  Resource  Productivities  Report.  These  reports  are  gen¬ 
erated  for  each  substitution  group  by  commodity  by  mode. 

(9)  Trip  Payload  Report.  This  report  is  generated  for  each 
substitution  group  only  if  vehicles-per-time-period  link  constraints 
are  used. 

(10 )  Origin-Destination  Data  Report  **or  a  Minimization  Run. 

A  maximization  run  would  not  show  delivery  requirements. 

c.  Output  File.  The  output  file  of  the  Input  Module  becomes  the 
input  file  for  the  RACAT  Module.  It  is  a  binary  file  and  is  in  the 
following  sequence: 

(l)  Title  and  Run  Parameters. 


♦See  Enclosure  1  of  SD,  Volume  I,  User's  tonual,  for  examples  of 
computer-produced  reports  based  on  a  test  problem. 
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(2)  Network  Data. 

(3)  Node  Constraints  (if  any). 

(4)  Resource  Inventories. 

(5)  Resource  Normalized  Prices  and  Smallest  Non-Zero  Price 
(Min- cost  operating  objective  only). 

(6)  Origin/Destination  Pair  Indexes  for  List  Processing. 

(7)  Origin/Destination  Pair  Data. 

(8)  Delivery  Requirements  (minimi zation  objective  only). 

(9)  Substitution  Indexes  for  List  Processing. 

(10 )  Bit  Map  for  Resource  Productivities  Substitution  Changes 

(if  substitution  is  used).  ~ 

(11)  Resource  Productivities. 

(12)  TRIP  Payload  Factors  (if  link  capacities  are  in  v/TP). 

(13)  Coi.version  Factors  for  Commodities. 

(14)  Output  Titles  (for  Output  Module  only). 

(15)  Table  of  Unique  Nodes  (for  Output  Module  only). 

3-  Methods .  The  logic  utilized  by  the  Input  Module  consists  basically 
of  eight  phases  of  operations  (or  subroutines): 

Read  input  data  (INEDIT)  . 

Sort  and  edit  data  (SORTAR)* 

Set-up  origin/ destination  pairs  (ODRCWS). 

Write  input  reports  and  file  for  RACAT  Module  (WRITES). 

Convert  requirements  and  productivities  (CONVER). 

Set  up  table  of  unique  node  names  (NODIIDC). 

Perform  network  update  (MODIFY). 

Perform  and  report  on  error  checks  (ERROR). 

Each  phase  is  discussed  below  and  is  shown  in  the  Input  Module  Functional 
Flow  Chart,  Figure  3. 

a.  Read  Input  Data  (INEDIT).  The  functions  of  the  INEDIT  routine 
are  to  enter  the  input  data  specified  in  the  previous  section  into  the 
appropriate  arrays  and  to  monitor  the  editing  subroutines. 

(l)  Basic  network  data  are  entered  and  stored  in  the  following 

ways  : 
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MATARC  (I,  I) 
MATARC  (i,  2) 
MATARC  (I,  3) 
COND  (T) 

COST  (I) 

RHS  (I) 

II 

NCD  (20) 


NC0M4D  (i) 


=  Arc  No. 

=  FROM  node 
=  TO  node 
=  Mode  of  arc 
=  Arc  condition 
=  Arc  cost  (distance) 

=  Arc  capacity 

=  Direction  (0  =  undirected;  1  =  directed) 

=  Commodity  availability  (0  =  will  not 
handle  commodity  k;  1  =  will  handle 
commodity  k). 

--  The  array  in  which  the  direction  (bit  21) 
and  the  commodity  evailability  (bits  1-20, 
right  to  left  in  the  CDC  6400  and  right 
to  left  in  the  IEM/360)  bits  are  stored 
for  arc  I. 


('  )  A  packing  technique  is  used  for  the  arc  direction  and 
commodity  availability  matrix  in  order  to  minimize  core  utilization. 

The  bits  pertaining  to  arc  I  are  packed  in  NCOMMD  (i).  The  firct  20 
bits  of  e  word  are  used  for  the  20  commodities;  i.e.,  commodity  1  =  bit 
1,  etc.  If  the  particular  arc  can  accommodate  commodity  k  then  bit 
number  k  of  the  word  is  turned  on  (set  to  l);  if  it  cannot,  it  is  left 
off  (zero).  An  additional  bit  is  used  to  specify  arc  direction.  The 
21st  bit  of  a  word  is  for  symmetry.  If  the  arc  is  directed,  the  bit 
is  turned  on  (set  to  l);  if  the  arc  is  undirected,  it  is  left  off  (zero). 
See  Chapter  3  for  more  details. 

(3)  Delay  times  on  the  arcs  are  computed  by  the  equation 


TRTIME  (I)  = 


COND  (I)  x  COST  (I] 
RATES  (j) 


where , 

I  =  Arc  row  number 
J  =  Mode  of  the  arc 

(4)  Run  parameters  are  set  in  the  INEDIT  subroutine  according 
to  the  following  dimension'’ ; 

NARC  =  Number  of  arc  rows  (£700 ) 

NRES  =  Number  of  resource  constraints  (£50 ) 

NODR  =  Number  of  movement  requirement  rows  (*150 ) 

NRCF  =  Total  number  of  rows  in  the  basis  (£900) 

NCONST  =  Number  of  node  constraints  (£50 ) 


Note  that  NROW  =  NARC  +  NRES  +  NODR  +  NCCNST  could  equal  950  if  all  were 
at  their  maximum.  This  is  not  allowed.  NROW  is  restricted  to  be  £900. 
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b.  Sort  and  Edit  Data  (SORTAR).  Some  elementary  editing  operations 
are  performed  on  the  input  data  upon  completion  of  reading  the  network 
data  and  all  other  input  data. 

(1)  Network  Data  Sorting.  An  internal  sort  routine  1b  used 
to  sort  the  entire  network  data  on  the  FROM  node.  The  sorting  process 
is  used  to  increase  the  efficiency  of  the  RA.CAT  Module.  The  sort  rou¬ 
tine  used  is  a  "Shell  Sort"  (developed  hy  C.  Shell)  modified  to  acco¬ 
mmodate  the  special  data.  This  routine  is  extremely  fast  and  effectively 
has  as  its  maximum  number  of  comparisons : 

(N  Logg  N] 

where  N  =  Number  of  elements  of  the  array  to  be  sorted. 

Run  indices  are  computed  based  on  the  input  data  as  follows : 

MP  =  NARC  +  NRES 

NODRCW  =  MP  +  NCONST 

NROW  =  NOURCW  +  NODR 

(2)  Resource  Productivities  Sorting.  The  modified  Shell  Sort 
routine  is  used  again  to  sort  the  resource  productivities.  The  sort 
hierarchy  is  as  follows : 

Substitution  group 
Commodity  number 
Mode  number 

Type  of  resource. 

Thus  the  MASRES  matrix  is  sorted  as  is  the  RESMAS  vector  containing  the 
resource  productivities  to  correspond  with  the  proper  elements .  Indices 
are  established  to  define  the  first  and  last  entry  of  each  substitution 
group.  The  resource  type  numbers  in  substitution  groups  two  and  three 
are  reduced  to  the  range  1-50  by  subtracting  50  and  100  respectively. 

This  yields  the  actual  index  of  the  resource. 

(3)  Resource  Cost  Comparisons.  The  cost  of  the  resources 
(productivity  multiplied  by  the  price)  in  each  substitution  group  is 
compared.  The  least  cost  resource  is  placed  in  group  one,  second  in 
group  two,  most  costly  in  group  three.  Other  data  which  are  given  by 
substitution  group  are  rearranged  accordingly.  This  is  done  to  aid 
the  RACAT  Module  in  obtaining  a  lower  priced  feasible  solution  at  the 
end  of  Phase  I  than  it  would  otherwise  during  a  min- cost  run  with  sub¬ 
stitution.  of  resources.  The  resource  prices  have  no  influence  on  the 
Phase  I  objective  function.  Consequently,  the  first  available  resource 
set  is  used.  With  the  resources  rearranged  as  described  here,  the 
cheapest  resource  set  will  be  used  first. 
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c.  Set  Up  Origin/Destination  Pairs  (ODROWS). 

(1)  The  function  of  subroutine  ODROWS  is  to  manipulate  the 
origin/destination  (o/D)  pair  data  and  to  establish  the  appropriate  rows 
in  the  basis  for  a  minimization  run.  First  the  o/D  pair  data  are  sorted 
by  commodity.  The  storage  assignments  are  as  follows: 

NORDES  (i,  l)  =  Commodity  number 
NORDES  (I,  2)  =  Origin  node 
NORDES  (I,  3;  =  Destination  node 

NORDES  (I,  4)  =  Row  number  for  basis  (for  minimization 
run) 

(2)  The  "fixed"  delivery  requirement  for  each  pair,  if  speci¬ 
fied,  is  entered  into  DEIREQ  (20  +  I),  where  I  =  number  of  o/D  pair. 

The  first  20  words  of  DELREQ  Eire  reserved  for  any  "open"  delivery 
requirements  for  the  £20  commodities.  Thus  an  o/D  pair  may  be  "open", 
where  a  delivery  requirement  ■'s  not  specified,  or  "fixed"  where  a 
delivery  requirement  is  specified  in  the  input  data. 

(3)  The  limits  of  the  O/D  pairs  are  then  obtained  and  stored 
in  (NCD  (I),  I  =1,  20 )  for  the  £20  commodities.  These  indices  are 
used  to  process  the  list  in  the  RACAT  Module. 

(4)  For  a  minimization  run,  the  O/D  pair  data  are  sorted  within 
each  commodity  by  "open"  and  "fixed."  "Open"  o/D  pairs  are  assigned  a 
rov,  number  in  the  basis  corresponding  to  the  particular  commodity.  The 
right-hand-side  delivery  requirement  for  this  row  is  the  amount  to  be 
delivered  for  the  "open"  O/D  pairs  of  that  commodity.  The  "fixed"  O/D 
pairs  are  assigned  individual  row  numbers  in  the  basis  since  they  have 
specific  right-hand-side  delivery  requirements.  The  right-hand-side 
vector  for  the  o/D  pair  rows  is  then  built  on  the  bas  ..3  of  hese  row 
assignments . 

d.  Write  Input  Reports  and  File  for  RACAT  Module  (WRITES).  This 
routine  is  used  to  generate  the  reports  and  file  discussed  in  the 
previous  section  under  "Output  Reports"  and  "Output  Files".  The  methods 
used  are  straight-forward  report  routines. 

e.  Convert  Requirements  and  Productivities  (CQNVER). 

(l)  This  routine  performs  the  conversion  of  delivery  require¬ 
ments  and  productivities  for  all  commodities  not  normally  specified  in 
the  same  units  as  the  link  capacities.  The  conversion  is  not  made  if 
the  capacities  are  expressed  in  vehicles-per-time-period  (V/TP),  since 
in  that  formulation  the  commodities  are  accepted  in  their  natural  units. 
Commodities  such  as  passengers  and  POL  are  better  specified  in  their 
natural  units  as  people  and  barrels  respectively.  When  not  in  the 
v/TF  formulation,  it  is  still  convenient  to  input  the  commodities  and 
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their  corresponding  productivities  and  movement  requirements  in  terms 
of  their  natural  units.  Thus  a  conversion  factor  is  used  by  the  pro¬ 
gram  to  equate  all  the  commodities  to  a  single  equivalent  unit  - 
usually  short  tons. 

(2)  The  delivery  requirements  for  a  commodity  requiring  con¬ 
version  are  divided  by  the  specified  factor.  For  example,  70  passengers 
to  be  moved  with  a  conversion  factor  of  7  would  be  set  equal  to  10 
equivalent  short  tons  by  CCNVER. 

(3)  Resource  productivities  applying  to  any  commodity  requir¬ 
ing  conversion  are  multiplied  by  the  specified  factor.  For  example, 

a  troop  carrier  productivity  of  .001  (40  men,  500  miles  per  day  includ¬ 
ing  a  return  trip  empty)  would  be  converted  to  .007  vehicles  per  ton 
mile. 


f .  Perform  Network  Update  (MODIFY ) . 

(l)  As  mentioned  earlier,  the  ARCC  data  card  provides  the  cap¬ 
ability  of  updating  certain  network  data.  There  are  8  "fields"  of  arc 
data  that  are  defined  as  follows: 

Field  Description 

1  FROM  node 

2  TO  node 

3  Mode  of  arc 

b  Condition  of  arc 

5  Cost  (distance)  of  arc 

6  Toll  of  arc 

7  Capacity  of  arc 

8  Direction  of  arc 

The  desired  update  of  most  of  the  above  fields  is  accomplished  by  the 
use  of  a  linear  combination.  The  FORTRAN  expression  is 

[F*]  =  Q  •  (F  ]  +  X 

where 

[Fj]  =  The  data  field  to  be  modified 

[Fj]  =  Any  data  field  (i  =  j  is  acceptable) 

Q  =  Multiplier 
X  =  Value  to  be  added 

Thus  the  new  value  of  field  i  is  set  equal  to  the  old  value  of  field  j 
times  a  multiplier  plus  seme  value. 
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(2)  Provision  is  also  made  for  specifying  limits  between  which 
the  updating  is  to  be  accomplished.  This  is  performed  by  the  inequality: 

Y  £  [F  ]  s  Z 

where 

[Fi ]  =  The  data  field  to  which  the  limits  apply 

Y  =  The  lower  limit  for  modification 
Z  =  The  upper  limit  for  modification 

Thus  one  may  update  field  Fj  for  the  entire  network  or  within  specified 
limits  as  indicated  by  the  inequality.  The  [F. ]  in  the  inequality  may 
refer  to  smother  field  if  desired. 

(3)  For  example,  assume  it  is  desired  to  impose  a  condition  of 
1.0  to  all  arcs  of  the  network  which  are  mode  6.  We  would  use  both 
equations  as 


F4  =  0.  ‘  F4  +  1.0 


and 


6  s  F3  s  6 

which  would  sot  data  field  4  (condition  of  arc)  equal  to  zero  times 
field  4  plus  1.0  for  all  arcs  which  meet  the  criterion  that  data  field 
3  (mode  of  arc)  have  a  value  26  and  s6. 

(4)  The  above  expressions  are  applicable  to  any  combination  of 
the  third  through  eighth  network  data  fields.  Fields  1  and  2  cannot  be 
used  because  they  are  alphanumeric  fields  (FROM  and  TO  nodes).  If  the 
updating  is  performed  on  the  integer  field,  field  3  (mode),  the  new 
value  is  rounded. 

g.  Set  Up  Table  of  Unique  Node  Names  (RODIDX).  The  NODIBX  routine 
finds  the  unique  nodes  (there  must  be  at  least  four  in  the  network  for 
the  binary  search  technique  to  work  properly)  and  places  them  in  the 
array  LA.HTA3  (j).  The  index  number  of  the  array  is  then  placed  in 
MATARC  (I,  l)  and  MATARC  (i,  2),  the  FROM  and  TO  nodes,  and  written  on 
the  file,  TAPE1,  for  use  in  the  RACAT  Module.  An  error  message  is  gen¬ 
erated  if  the  array  is  >400.  The  nodes  that  are  constrained  and  the 
0/D  pair  nodes  are  assigned  their  respective  index  numbers.  While  both 
the  Input  and  Output  Reports  print  alphanumeric  node  names,  RACAT  uses 
integer  numbers. 

h.  Perform  and  Report  on  Error  Checks  (ERROR). 

(l)  An  input  data  error  routine  is  used  to  report  the  major 
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input  errors  and,  if  necessary,  terminate  the  run.  The  data  processing 
continues  until  all  data  are  completed  in  order  to  report  all  the  dis¬ 
covered  errors  in  one  pass.  The  specific  error  messages  relating  to 
each  type  of  input  data  are  mentioned  in  the  section  in  this  chapter 
on  "Read  Input  Data". 

(2)  In  addition  to  the  errors  by  card  type,  the  routine  checks 
for  possible  run  parameter  limits,  proper  amounts  of  output  titles,  and 
excessive  row  entries  and  will  generate  error  messages  if  appropriate. 
The  error  messages  are  reported  before  the  normal  input  data  reports 
are  produced.  Upon  completion  of  the  input  data  reports,  the  error 
routine  decides  whether  any  discovered  error  is  a  major  one  and,  if 

so,  terminates  the  run. 

(3)  Several  self-explanatory  error  messages  are  reported  in 
other  routines. 
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CHAPTER  3.  RACAT  MODULE 

1.  Purpose  and  Scope.  The  purpose  of  the  RACAT  Module  is  to  solve  the 
multi-commodity  network  flow  problem  with  mobility  resource  allocation 
as  formulated  by  the  Input  Module  in  the  data  file.  The  technique 
permits  the  maximization  of  flow  or  the  minimization  of  cost  (or  time) 
subject  to  the  constraints  of  the  capacity  of  the  network  and  the 
availability  of  resources.  Due  to  the  mathematical  sophistication  of 
the  algorithms,  this  chapter  will  be  confined  to  a  discussion  of  the 
general  methods  used  to  solve  the  problems.  Specific  information  with 
respect  to  the  algorithms  will  be  found  in  Enclosures  1,  2  and  3- 

2.  Inputs  and  Outputs.  Inasmuch  as  the  RACAT  Module  is  a  set  of 
algorit’ ms  for  solving  specific  problems,  the  inputs  and  outputs  are 
well-defined  and  structured  so  as  to  facilitate  efficiency  toward 
arriving  at  solutions.  This  section  describes  the  input  data  require¬ 
ments  and  output  files  generated  by  the  RACAT  Module. 

a.  Input  Data  Requirements.  There  are  two  types  of  input  data 
required  to  operate  the  RACAT  Module:  first,  the  data  file  (TAPfil) 
written  by  the  Input  Module;  and  secondly,  one  card  specifying  certain 
parameters  and  the  mode  of  operation  to  be  used  —  NORMAL  or  RESTART. 
However,  the  restart  procedure  requires  more  data  and  will  be  discussed 
later. 


(1)  Data  File.  The  contents  and  order  of  the  data  file  are 
specified  in  Chapter  2,  INRJT  MODULE.  The  file  name  is  TAPE1  and  is 
written  in  binary  logical  records.  As  discussed  in  Chapter  2,  the  data 
are  structured  specifically  for  RACAT. 

(2)  RACAT  Data  Card.  One  punched  card  is  required  to  operate 
the  RACAT  Module  and  is  used  to  specify  the  operating  mode  and  optional 
control  parameters.  Two  modes  of  operation  are  possible  -  NORMAL  and 
RESTART.  These  are  discussed  in  detail  under  the  section  of  this 
chapter  entitled  "Methods."  The  format  of  the  RACAT  data  card  and  the 
definitions  of  the  parameters  are  as  follows: 


Columns 

Information 

1-7 

NORMAL  or  RESTART 

11-15 

NITER 

16-20 

INVC 

21-25 

IXPRNT 

26-30 

NREPRT 

31-35 

ITSAVE 

36- 40 

IREST 

41-45 

IRTARC 

46-50 

IRTRES 

51-55 

ISTMI 

NITER 

Maximum  number  of 

iterations  until  job  is  terminated. 
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INVC  Frequency  of  inversion. 

IXPRNT  Frequency  of  debugging  print. 

NREFRT  Iteration  number  when  solution  report  is  to  be  obtained 

(before  optimum)  and  job  terminated. 

ITSAVE  Frequency  of  saving  basis  and  bit  map  on  TAPE9* 

IREST  Equal  to  0,  1,  2  or  3 

0  =  Normal  restart  from  file  and  mark  slacks. 

1  =  Restart  from  file  and  do  not  mark  slacks. 

2  =  Restart  from  cards  and  mark  slacks. 

3  =  Restart  from  cards  and  do  not  mark  slacks. 

IRTARC  Equal  0  or  1 

0  =  Do  not  print  arcs  on  iteration  log. 

1  =  Print  arcs  on  iteration  log. 

IRTRES  Frequency  of  printing  resources  on  the  iteration  log. 

ISIMP  Frequency  of  printing  simplex  multipliers. 

The  parameters  on  the  RACAT  data  card  are  optional;  if  not  desired, 
the  appropriate  data  field  is  left  blank  or  set  to  0.  NITER  and  INVC 
may  have  already  been  set  in  the  Input  Module.  Those  values  may  be 
retained  by  leaving  the  appropriate  data  fields  blank  or  setting  them 
to  0.  They  may  be  rese:t  by  specifying  their  new  values  on  the  RACAT 
data  card. 

b.  Outputs.  Four  output  files  are  generated  by  the  RACAT  Module. 

(1)  The  data  file  (TAPEl)  is  used  to  output  certain  informa¬ 
tion  obtained  by  the  optimal  solution.  The  data  listed  below  are 
added  to  the  data  file  (TAPEl)  when  an  optimal  solution  is  reached. 
These  are  explained  in  detail  under  "Methods"  in  this  chapter. 

(a)  The  vector  identification  numbers  present  in  the 
basis  at  optimization. 

(b)  The  solution  vector. 

(c)  The  simplex  multipliers  (shadow  prices). 

(d)  The  value  of  the  objective  function. 

(e)  The  number  of  iterations  required  to  reach  the 

solution. 

(2)  In  addition,  a  vector  file  (TAPE2)  is  created  during  the 
algorithm  which  specifies  the  vectors  as  they  are  generated  to  enter 
the  basis.  The  creation  and  maintenance  of  the  vector  file  is 
discussed  under  "Methods"  in  this  chapter.  TAPE2,  a  binary  file,  is 
structured  as  follows: 

(a)  Index  number  of  this  vector 

(b)  Origin/ destination  number  of  this  vector. 
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(c)  Row  number  of  delivery  requirement  intersect. 

(d)  Chain  delay  time. 

(e)  Number  of  links  in  the  chain. 

(f)  Link  row  numbers  in  the  chain. 

(g)  The  vehicle  payload  factors  (if  link  capacities  are 
expressed  in  V/TP). 

(h)  Resource  requirements  per  unit  flow  on  this  chain. 

(i)  Constrained  nodes  in  the  chain,  if  any. 

(3)  A  save-basis  file  (TAPE9)  is  created  for  restart  purposes. 
This  file  is  explained  under  the  Restart  section  of  'Methods."  TAPE9, 
a  binary  file,  is  structured  as  follows: 

(a)  The  number  of  structural  vectors  (or  routes)  in  the 

basis. 

(b)  The  identifying  number  of  each  structural  vector  in 

the  basis. 

(c)  Tne  bit  map  identifying  the  slack  vectors  in  the 

basis. 

(4)  A  POP  file  (TAPE14)  is  created  for  use  by  the  Post  Optimal 
Package  Module.  This  file  iG  a  binary  file.  It  is  discussed  further 
in  Chapter  5.  During  the  RACAT  computer  run,  TAPE14  is  used  as  the  eta 
file,  where  the  product  form  of  the  inverse  is  stored  if  it  cannot  be 
contained  in  core.  The  structure  of  the  file  TAPE14  at  a  normal  ending 
of  a  RACAT  run  is  as  follows: 

(a)  Hie  product  form  of  the  inverse  of  the  basis  (count 
parameters,  bit  map  and  packed  eta  vectors  with  the  zeros  removed). 

(b)  Type  of  record  indicator,  run  title,  necessary 
parameters  (NARC,  NRES,  NCONST,  NROW,  ITER,  ISORTI,  NCOM)  and 
conversion  factors. 

(c)  Right-hand- side  values. 

(d)  Vector  identification  numbers  present  in  the  optimum 

solution. 

(e)  The  solution  vector. 

(f)  0/D  data. 
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3.  Methods.  The  RACAT  Module  set  of  algorithms  for  solving  the  multi- 
commodity  network  flow  problem  with  resource  allocation  has  been 
designed  to  operate  in  either  a  maximization  of  flow  or  a  minimization 
of  cost  (or  time)  objectives.  To  accomplish  the  optimization,  a  series 
of  18  functional  subroutines  is  used.  Each  of  these  subroutines  is 
discussed  in  this  section.  These  discussions  are  general  in  nature  and 
describe  the  algorithms  in  a  functional  framework.*  Three  bit  manipu¬ 
lating  subroutines  axe  also  used  and  are  briefly  described. 

a.  Subroutines.  These  functional  subroutines  are  discussed  in 
this  section  by  describing  the  algorithms  for  solving  the  two  operating 
objectives  of  the  model  -  maximization  of  flow  and  minimization  of 
costs  --  and  defined  in  terms  of  their  primary  function  as  follows: 

(1)  RACAT.  Reads  and  loads  the  input  data  file. 

(2)  SETUP.  Sets  up  the  algorithms  for  the  appropriate  opera¬ 
ting  objective. 

(3)  PEAS  El .  Obtains  a  feasible  solution  to  a  minimization 
problem  (Riase  i). 

(4)  SMPLEX.  Generalized  product  form  of  the  inverse  linear 
programming  algorithm. 

(5)  ALHIAS.  Computes  simplex  multipliers. 

(6)  MAXIO.  Controls  shortest  chain  algorithm  for  a  maximiza¬ 
tion  problem. 

(7)  MINIO.  Controls  shortest  chain  algorithm  for  Riase  II  of 
a  minimization  problem. 

(8)  SUBST.  Computes  resource  vector  each  iteration  for 
shortest  chain  costs. 


0 


(9) 

subproblem. 

SCHAIN. 

Shortest  chain  algorithm  for  column  generation 

(10) 

SORTAJ. 

Internal  sort  for  inversion  routine. 

(11) 

REINV. 

Operates  inversion  routine. 

(12) 

RESTART 

.  Operates  restart  operation  to  introduce  vectors 

into  the  bas 

is . 

* 

See  Enclosures  1,  2  and  3  for  the  mathematical  derivation  of  the 
RACAT  Algorithm. 

L 
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(13)  WRITE.  Genera  J  output  at  optimum. 

(14)  ERROR.  Terminates  RACAT  due  tc  some  major  error  during 
algorithm. 

(15)  NODBC.  Apportions  the  simplex  multipliers  of  the 
constrained  nodes  to  the  appropriate  link. 

(16)  S BASIS.  Saves  the  basis. 

(17)  SORTIE.  Supplies  the  vehicle  payload  factors  for 
constrained  nodes  and  the  shortest  chain. 

(18)  ETAFLE.  Manages  the  basis  inverse  in  product  form  and 
packed  format.  Die  basis  inverse  is  modified  from  iteration  to 
iteration.  The  modification  is  in  the  form  of  a  multiplication  by  an 
elementary  column  matrix.  The  elementary  column  matrix  has  one  column 
which  differs  from  the  identity  matrix  of  the  same  dimensions.  This 
column  is  known  as  the  eta  vector  corresponding  to  that  elementary 
column  matrix.  Only  the  index  giving  the  position  of  the  eta  vector 
in  the  elementary  column  matrix  and  the  eta  vector  itself  need  be 
stored  every  iteration.  Usually  the  eta  vector  is  very  sparse,  that 
is,  it  has  many  zero  entries.  ETAFLE  stores  only  the  non- zero  numbers 
along  with  a  bit  map  indicating  their  positions  in  the  eta  vector.  It 
also  retrieves  these  numbers  when  needed  and  makes  the  necessary  compu¬ 
tations  of  only  the  non- zero  numbers  as  appropriate  each  iteration  to 
conserve  computer  time.  ETAFLE  is  discussed  in  detail  at  this  point 
because  its  function  is  similar  for  all  operating  objectives.  It  is 
called  by  the  subroutine  ALPHA  in  computing  the  simplex  multipliers, 

by  SMPLEX  in  removing  and  entering  a  variable  into  the  basis  and 
generating  the  corresponding  eta  vector.  REINV  and  RESTART  also  use 
it  to  generate  a  new  eta  file. 

(19)  IONBIT,  IOFFBIT,  KBIT.  These  are  bit  manipulating  sub¬ 
routines.  The  arguments  for  all  three  are  the  same,  (NWORD,  NLMBIT). 
IONBIT  sets  the  bit  in  NWORD  given  by  NLMBIT  to  1.  IOFFBIT  sets  it  to 
0.  KBIT  tests  that  bit  whether  it  is  on  or  off. 

b.  Maximization  Objective.  Die  simplest  operating  objective  is 
that  of  flow  maximization.  Figure  4  is  a  simplified  flow  chart  of  the 
RACAT  Max- Flow  Algorithm.  The  problem  as  formulated  in  words  is  to 
determine  the  maximum  flow  over  all  possible  chains  from  origins  to 
their  respective  destinations  subject  to  the  capacity  constraints  of 
the  arcs  and  the  availability  of  resources.  The  RACAT  Module  is 
designed  to  solve  the  maximization  problem  by  an  iterative  sequence  of 
routines  to  accomplish  the  algorithm  as  specified  in  Figure  4.  Diis 
sequence  is  discussed  below.  Figure  5  is  a  general  flow  chart  of  the 
RACAT  Module  max-flow  process  showing  the  functional  subroutines 
sequence  for  a  normal  run.  The  functions  of  each  subroutine  are  as 
follows : 
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Form  an  initial  basis 

=  +  p  +  t 


Calculate  B' 


Calculate  the  simplex  multipliers 
c  g  B  =  ( i  <  &  m  (  t] » ••• »  ffp  / 

P\ i  •••<  Pt) 


Are  any  of  the  simplex 
multiplier?  negative? 


Calculate  dj*,  the 
cost  of  each  arc  i  per 
commodity  It  in  terms  of 
the  simplex  multipliers 
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appropriate 
slack  vector 


Solution 
is  optimum 


FIGURE  4.  SIMPLIFIED  FLOW  CHART  OF  RACAT 
MAX-FLOW  ALGORITHM 
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(1)  RACAT.  The  RACAT  subroutine  is  used  to  load  the  input 
data  file  into  the  proper  storage  locations.  It  reads  TAPE1,  the  data 
file  created  by  the  Input  Module  discussed  in  Chapter  2.  No  editing  is 
performed  since  the  Input  Module  structures  the  input  file  specifically 
for  RACAT.  It  sets  up  various  tables  of  the  network  data  to  be  used 

by  the  shortest  chain  algorithm. 

(2)  SETUP.  This  subroutine  performs  initialization  of  the 
appropriate  variables,  vectors,  and  matrices  and  then  sets  up  the 
operating  objective.  The  algorithm  starts  with  an  identity  matrix 

as  a  beginning  basis  and  a  solution  vector  equal  to  the  right-hand- side 
values.  The  number  of  rows  for  a  maximization  is  the  following  sum. 

NRCW  --  NARC  +  NRES  +  NCONST 

where , 

NARC  =  Number  of  arc  constraints 
NRES  =  Number  of  resource  constraints 
NCONST  =  Number  of  node  constraints. 

(3)  MAXIO.  This  subroutine  operates  as  the  controller  for  a 
maximization  run. 

(a)  Its  functions  each  iteration  (see  Figure  6)  are 
delineated  as: 


1.  Initialize  variables. 

2.  Access  an  origin/destination  pair. 

3.  Call  SUBST  routine  for  appropriate  resource  vector. 

4.  Assign  pseudo-costs  to  arcs  of  the  network. 

jj.  Call  SCHAIN  routine  to  generate  a  vector. 

6.  Call  SORTIE  routine  if  V/TF  constraints  formulation 
is  used. 

7.  Call  NODEC  routine  if  node  constraints  are  present. 

8.  Evaluate  vector  to  enter  basis  against  simplex 
criterion. 

Q.  Repeat  steps  2-6  for  all  0/D  pairs  to  obtain  the 
"best"  vector  to  enter  the  basis. 


10.  Write  this  vector  on  the  vector  file. 


•<f 
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11.  Call  SMPLEX  routine  to  update  basis. 

12.  If  no  vector  meets  the  criterion,  terminate 
algorithm. 

(b)  Each  vector  file  record  is  written  as  follows: 

1.  NAJ  =  Number  of  vector.  A  sequential  numbering 
scheme  beginning  with  NRCW  +  1. 

2.  NP  =  Number  of  the  origin/destination  pair 
associated  with  this  vector. 

3.  NR  =  No  useful  information  for  the  "Maximization 
Objective." 

4.  CTME  =  The  delay  time  for  this  chain. 

5-  NAR  =  Number  of  arcs  contained  in  the  chain. 

6.  (AJ(l),  I  =  1,  NAR)  =  Arc  row  numbers  in  the  chain. 

7*  (AJ2(l),  1=1,  NAR)  =  Vehicle  payload  factor  for 

the  arc  whose  row  number  is  in  the  chain. 

8.  ( COLUMN ( I ) ,  I  =  NARC  +  1,  NODRCW)  =  Resource 

coefficients  and  constrained  nodes  for  this  chain. 

(4)  SUBST.  This  subroutine  performs  the  necessary  calculations 
each  iteration  to  compute  the  appropriate  resource  vector  for  each 
transportation  mode  for  a  commodity  number.  If  ISUB  =  1  (i.e.,  no 
substitution  alternatives  permitted)  the  computation  involves  merely 
processing  the  master  list  of  resource  productivities  (MASRES,  RESMAS) 
and  building  the  appropriate  vector.  If  ISUB  =  2  or  ISUB  =  3,  each 
feasible  set  of  resource  productivities  is  evaluated  against  the  current 
simplex  multipliers  to  obtain  the  best  alternatives  from  which  to 
generate  pseudocosts  for  the  arcs.  Figure  7  shows  a  general  flow 
chart  of  the  SUBST  routine  used  in  the  maximization  c- srating  objective. 

(5)  SCHAIN.  This  subroutine  consists  of  a  shortest  chain 
algorithm  described  by  Dreyfus-'-  as  the  computationally  most  efficient 
procedure  and  credited  to  E.  W.  Dijkstra.  The  problem  is  that  of 
finding  a  minimal  cost  chain  from  one  node  to  another  in  a  network  in 
which  each  arc  (i,  j)  has  an  associated  c o-t  (or  length)  d(i,  j)  s  0. 


^See  Dreyfus,  S.  E.,  An  Appraisal  of  Some  Shortest-path  Algorithms, 
Rand  Corporation,  Memo.  RM-5433-1-PR,  September  19^8. 


3-10 


FIGURE  7.  GENERAL  FLOW  CHART  OF  SUBST  ROUTINE-MAXIMIZATION  OBJECTIVE 


Thus  the  restriction  that  all  costs  d(i,  j)  £  0  precludes  the  possibility 
of  negative  cycles  and  enables  the  algorithm  to  find  the  minimal  cost 
chain. 


(a)  Figure  8  illustrates  the  shortest  chain  algorithm 
used  in  the  SCHAIN  routine.  The  algorithm  has  been  modified  to  choose 
the  chain  with  the  minimum  distance  in  case  of  ties.  This  method  has 
in  most  cases  proved  to  be  more  efficient  with  respect  to  speed  in 
obtaining  an  optimal  solution  to  the  larger  problem.  An  alternate 
method  is  to  choose  the  chain  with  maximum  capacity  in  case  of  a  tie. 
This  method  is  especially  suited  for  a  max-flow  problem.  Thus  the 
problem  structure  determines  which  of  the  two  methods  is  best.  The 
origin  and  destination  nodes  are  specified  and  the  pseudocosts  applied 
to  the  arcs  are  computed  in  the  MAXIO  routine. 

(b)  The  SCHAIN  routine  checks  each  arc  before  it  is 
labeled  to  insure  the  arc  can  accommodate  the  commodity  number  with 
which  the  origin/destination  pair  is  associated.  This  procedure  is 
accomplished  by  inspection  of  the  appropriate  bit  as  "packed"  into  a 
word  which  is  discussed  in  Chapter  2,  INRJT  MODULE.  If  the  appropriate 
bit  is  "on,"  the  arc  can  accommodate  the  particular  commodity  and  is  a 
candidate  for  the  shortest  chain;  if  the  appropriate  bit  is  "off,"  the 
arc  cannot  accommodate  the  commodity  and  thus  is  not  labeled.  Checking 
for  directed  arcs  is  handled  similarly. 

(c)  An  error  message  is  produced  and  the  algorithm  is 
terminated  if  no  chain  can  be  found  between  any  origin  and  destination. 
This  indicates  a  deficiency  in  the  network  data  and  results  in 
termination  of  the  RAC AT  Module. 

(6)  SMPLEX.  The  SMPLEX  subroutine  is  a  linear  programming 
code  based  on  the  product  form  of  the  inverse.  The  algorithm  updates 
the  basis  inverse  (product  form)  by  the  use  of  an  eta  vector  file.^- 
Figure  9  illustrates  the  general  logic  of  the  SMPLEX  routine.  Much 
of  the  eta  file  is  kept  in  core,  if  it  gets  too  large  for  in-core 
operation,  it  is  written  on  the  binary  file  TAPE14.  Each  iteration 
the  eta  file  is  read  to  calculate  the  proper  row  to  leave  the  basis 
and  a  new  eta  vector  is  computed  and  written  on  the  file.  In  the 
maximization  mode  the  "cost"  of  a  unit  flow  is  unit;  thus  the  coeffi¬ 
cients  are  updated  in  this  fashion.  'The  objective  function  is  then 
computed  by 

z  =  E  c.  x 

j  J  J 


xSee  G.  Hadley,  Linear  Programming,  Addis  on- Wes ley  Publishing  Co. , 
Inc.,  Reading  Mass.,  19&2,  for  a  complete  discussion  of  the  revised 
simplex  method. 
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Start  SCHAIN 


FIGURE  9.  GENERAL  FLOW  CHART  OF  SMPLEX  ROUTINE 


where, 


c  =  "cost"  of  vector  j  (unity  if  vector  j  is  in  the  basis; 

J  zero  if  vector  j  is  not) 

Xj  =  flew  on  vector  j 

thus  the  value  of  z  is  the  total  flow;  at  optimum  it  is  the  maximum 
flow  attainable. 

(7)  ALPHAS.  This  subroutine  is  used  to  generate  the  simplex 
multipliers  each  iteration  after  SMPLEX  performs  the  basis  update. 

Figure  10  illustrates  the  ALPHAS  routine  for  the  maximization  objective. 
Note  that  if  any  simplex  multiplier  is  found  to  be  negative,  a  slack 
vector  is  created  and  immediately  entered  into  the  basis  by  the  SMPLEX 
routine. 


(8)  WRITE.  The  WRITE  routine  is  used  at  optimum  to  produce 
certain  optimal  data  which  are  required  for  the  output  reports.  As 
discussed  earlier  in  this  chapter  under  "Inputs  and  Outputs,"  these 
data  are  added  to  the  input  data  file  (TAFEl)  when  an  optimal  solution 
is  reached.  It  also  calls  the  routine  ETAFLE  to  write  ths  file  used 
by  the  POP  Module.  Definitions  and  content  of  the  data  written  on 
TAPE1  follow. 


(a)  Vector  I.D.  Numbers  in  Basis.  NCHARC  is  the  array 
of  vector  identification  numbers  present  in  the  basis  at  optimization 
with  their  corresponding  row  numbers.  As  discussed  under  MAXIO,  the 
vector  numbers  (NAJ)  are  assigned  in  sequential  fashion  beginning  with 
NRCW  +  1  as  they  are  entered  into  the  basis  during  the  algorithm. 

Since  some  vectors  may  leave  the  basis,  only  those  present  at  opti¬ 
mization  are  relevant  to  the  solution.  The  I.D,  numbers  permit  the 
Output  Module  to  identify  these  from  the  vector  file  (TAFE2)*. 

(b)  The  Solution  Vector.  The  solution  vector,  SOLVEC, 
is  defined  as  x^  =  B-1b 

where , 

B  ^  =  Basis  inverse 
b  =  Right- hand- side  vector. 

Since  the  SMPLEX  routine  uses  the  product  form  of  the  inverse,  the 
SOLVEC  values  are  updated  each  iteration  by  an  evaluation  of  the  eta 
vectors.  The  SOLVEC  values  represent  the  variables  in  the  basis;  i.e., 


♦See  the  section  of  this  chapter  on  "Restart  Mod''"  for  a  discussion 
of  how  these  I.D.  numbers  also  are  used  for  restarting. 
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FIGURE  10.  GENERAL  FLOW  CHART  OF  ALPHAS  ROUTINE 
MAXIMIZATION  OBJECTIVE 


the  values  of  all  active  vectors  are  obtained  from  SOLVEC.  The  utili¬ 
zation  of  these  values  is  discussed  in  Chapter  4,  OUTPUT  MODULE. 

(c)  The  Simplex  Multipliers  (shadow  prices).  Simplex 
multipliers  are  computed  each  iteration  by  subroutine  ALPHAS  by  reading 
the  eta  file.  At  optimum,  these  multipliers  represent  the  shadow 
prices  (or  movement  constraints)  of  each  constraining  row  of  the  basis. 

(d)  The  Value  of  the  Objective  Function.  As  discussed 
under  subroutine  SMPLEX  in  this  section,  the  objective  function  is 
computed  each  iteration.  At  optimum  it  represents  the  maximum  flow 
attainable  in  the  maximization  objective  and  it  is  reported  in  the 
Output  Module. 

(e)  The  Number  of  Iterations  Required  to  Reach  the  Solution. 
The  iteration  count  is  sequential  beginning  with  unity  and  is  reported 

in  the  Output  Module  for  information  purposes.  An  iteration  is  defined 
as  an  update  of  the  basis  (the  SMPLEX  routine)  required  either  in 
introducing  a  real  vector  or  entering  a  clack  vector  into  the  basis. 

c.  Minimization  Objective.  Figure  11  is  a  simplified  flow  chart 
of  RACAT  min-cost  algorithm  (Phase  I  and  II )  similar  in  form  to  Figure  h 
of  the  max- flow  process.  The  problem  as  formulated  in  words  is  to 
find  the  set  of  routes  and  the  allocation  of  resources  which  will 
satisfy  fixed  movement  requirements  at  minimum  total  cost  (or  time). 
Delivery  requirements  are  set  for  each  commodity  and  the  program  deter¬ 
mines  the  allocation  which  will  meet  these  requirements  at  minimum 
cost.  The  min-cost  solution  is  not  independent  of  time.  Since  the 
quantity  of  a  given  resource  used  is  :i  function  of  travel  time,  time 
is  minimized  for  a  given  mode  of  transportation.  However,  in  the 
min-time  operating  objective,  routes  are  selected  and  resources  assigned 
so  as  to  meet  the  delivery  requirements  with  minimum  total  flow 
unit-hours  (e.g.,  ton-hours).  The  minimization  objective  of  the  RACAT 
Module  solves  either  of  these  problems  with  the  same  algorithm;  thus 
this  discussion  is  applicable  to  both.  Figure  12  is  a  general  flow 
chart  of  the  RACAT  Module  min-cost  process  showing  the  subroutine 
sequence  for  a  normal  run.  The  functions  of  each  routine  are  discussed 
below. 


(1)  RACAT.  The  RACAT  subroutine  performs  the  same  functions 
in  the  minimization  problem  as  it  does  in  a  max- flow  problem. 

(2)  SETUP.  As  in  the  maximization  objective,  SETUP  performs 
initialization  of  the  appropriate  variables,  vectors,  end  matrices  and 
sets  up  for  the  operating  objective.  A  minimization  run  also  starts 
with  an  identity  matrix  as  a  beginning  basis  and  a  solution  vector 
equal  to  the  initial  right  hand  side  values.  However,  the  number  of 
rows  for  a  minimization  is  set  as 

NRCW  =  NARC  +  NRES  +  N CONST  +  NODP 
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FIGURE  12.  (continued) 
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NARC  =  Number  of  ere  constraints 


NRES  -  Number  of  resource  constraints 
NC0NS7  =  Number  of  nod-  constraints 

NODP  =  Number  of  origin-destination  movement  requirements 

(3)  PHASE1.  The  PHASE1  routine  operates  as  the  controller  for 
a  minimization  run  in  the  Phase  I  operation,  i.e.,  in  an  attempt  to 
attain  a  feasible  solution  to  the  problem.* 

(a)  As  such  it  performs  similar  functions  as  the  MAXIO 
routine  in  the  maximization  objective  arid  as  MINIO,  the  Phase  II 
controller  to  be  discussed  later.  Figure  13  is  a  general  flow-chart 
of  PHASE1.  The  function  of  PHASE1  at  each  iteration  is  to: 

1.  Initialize  variables. 

2.  Access  an  origin/destination  r,air  (in  sequence). 
Call  SUBST  routine  for  appropriate  resource  vector. 

4.  Assign  pseudocosts  to  arcs  of  the  network. 

5-  Call  SORTIE  routine  if  V/TP  constraints  formulation 
is  used. 

6.  Call  NODEC  routine  if  node  constraints  are  present. 

7.  Call  SCHAIN  routine  to  generate  a  vector. 

8.  Evaluate  vector  to  enter  basis  against  simplex 
criterion. 

If  the  vector  does  not  improve  the  objective 
function,  go  to  step  2. 

10.  Write  this  vector  on  the  vector  file. 

11.  Call  SKPLEX  routine  to  update  basis. 

12.  If  no  vector  meets  the  criterion,  terminate 
RACAT  run. 


*See  Enclosure  1,  2,  and  3  for  a  mathematical  discussion  of  the 
Phase  I  operation. 
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FIGURE  13.  GENERAL  FLOW  CHAFT  OF  PHASEI  ROUTINE 


13.  If  the  value  cf  the  objective  function  13  zero 
(i.e.,  the  sum  of  the  artificial  variables  is  0) 
terminate  FHASE1  with  a  feasible  solution. 
Otherwise  go  to  step  2. 

(b)  Each  vector  file  record  in  the  minimization  mode  is 
written  as  follows : 

NAJ  =  Number  of  vector.  A  sequential 

numbering  scheme  beginning  with 
NROW  +  1. 

NP  =  Number  of  the  origin/destination 

pair. 

NR  -  Number  of  the  row  in  the  basis  in 

which  the  0/D  intersect  occurs. 

This  row  number  identifies  the  g/D 
row  to  which  this  vector  applies. 

CTME  =  The  delay  time  for  this  chain. 

NAR  =  Number  of  arcs  contained  in  the 

chain. 

(Aj(l) ,  1=1,  NAR)  =  Arc  row  numbers  in  the  chain. 

(AJ2(l),  1=1,  NAR)  =  Vehicles  payload  factor  for  the 

arc  whose  row  number  is  in  the 
chain. 

(COUJMN(l),  I  = 

NARC  +  1,  NODRCW)  =  Resource  coefficients  and  con¬ 
strained  nodes  for  this  chain. 

(c)  As  indicated  in  the  FHASE1  flow  chart  (Figure  13)  a 
feasible  solution  is  attained  when  the  value  of  the  objective  function 
reaches  zero.  This  means  that  the  delivery  requirements  have  teen  met 
in  a  feasible  fashion  and  Phase  I  is  completed.  If  a  feasible  vector 
cannot  be  generated  at  any  iteration  to  enter  the  basis,  then  an 
infeasible  problem  has  been  encountered  and  the  RACAT  Module  terminates. 

(4)  SUBST.  Operation  of  the  SUBST  routine  for  the  minimization 
objective  is  similar  to  the  description  given  earlier  under  the  maximi¬ 
zation  objective.  The  criteria  used  for  selection  of  the  "best"  sub¬ 
stitution  alternative  is  the  only  difference  and  is  explained  in 
Enclosure  2. 

(5)  S CHAIN.  See  the  S CHAIN  discussion  under  "Jteximization 
Objective"  which  is  also  applicable  here. 
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(6)  SMPLEX.  The  discussion  in  ‘Maximisation  Objective"  of  the 
generalized  SMPLEX  linear  programming  code  based  on  the  product  form  of 
the  inverse  is  also  applicable  to  the  minimization  run  with  a  few 
exceptions.  Since  the  purpose  of  Phase  I  is  to  either  remove  the 
artificial  variab2.es  (requirements)  from  the  basis  or  to  have  them  in 
the  basis  at  a  zero  level,  costs  (c.)  of  all  structural  and  slack  variables 
are  set  to  zero  and  costs  of  unity  assigned  to  each  artificial.  In 
Phase  II  the  actual  cost  (c  )  is  assigned  each  legitimate  vector  and 
is  computed  each  iteration. ^  Thus  the  final  value  of  the  objective 
function  may  be  computed  by: 


as  before.  The  value  of  z  in  a  minimization  run  is  the  total  cost  of 
satisfying  the  delivery  requirements. 

(7)  ALPHAS.  The  ALPHAS  routine  performs  the  same  operations 

in  computing  the  simplex  multipliers  as  in  the  "Maximization  Objective." 
However,  if  any  simplex  multiplier  on  a  link  row  or  on  a  resource  row 
is  found  to  be  positive,  the  corresponding  slack  vector  is  created  and 
immediately  entered  into  the  basis  via  the  SMPLEX  routine. 

(8)  MINIO.  A  transition  operation  is  required  in  proceeding 
from  a  Phase  I  feasible  solution  to  a  Phase  II  processing.  This  involves 
a  simple  resetting  of  cost  coefficients  and  control  parameters.  MINIO 
then  assumes  control  for  the  Phase  II  solution  process.  It  performs 

the  same  general  functions  as  the  PHASE1  subroutine  in  that  it  accesses 
0/D  pairs  and  generates  a  vector  to  enter  the  basis.  Each  vector  to 
enter  is  written  on  the  vector  file  Just  as  is  done  in  Phase  I.  Phase  II 
is  terminated  when  the  MINIO  routine  is  unable  to  find  a  vector  qualified 
to  enter  the  basis.  That  is,  there  is  no  new  route  that  will  improve 
the  solution.  At  this  point  an  optimal  solution  for  the  minimization 
problem  has  been  reached.  Figure  14  illustrates  tls  MINIO  routine. 

(9)  ERROR.  The  ERROR  routine  is  called  when  either  of  ten  (10) 
major  errors  is  discovered.  The  type  of  error  is  reported  and  the 
RACAT  Module  is  terminated  with  a  report  of  the  active  vectors  present 
in  the  current  basis.  These  errors  are  in  general  only  discernible 
within  RACAT  and  indicate  a  major  problem  in  the  proceeding  of  the 
algorithm.  The  types  of  errors  checked  for  each  iteration  are  as  follows: 

(a)  The  value  of  the  objective  function  decreases  (maximi¬ 
zation  objective). 

(b)  The  value  ex'  the  objective  function  increases  (mini¬ 
mization  objective). 

(c)  An  infeasible  vector  enters  the  basis  -  indicates  an 
unbounded  problem. 
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FIGURE  14.  GFNERAL  FLOW  CHART  OF  MINIO  ROUTINE 
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(d)  A  chain  cannot  be  found  between  an  0/£  pair. 

(e)  An  infeasible  problem  as  stated  (minimization 

objective). 

(f)  Algorithm  has  exceeded  the  maximum  number  of 
iterations  as  stated  by  the  user  in  the  input  data,  (see  Chapter  2) 

(g)  Parity  error  on  eta  file  TAPElh  (for  CDC  6400). 

(h)  Tried  to  read  backwards  when  TAPElh  load  point 
(for  CDC  6400). 

(i)  Jhe  number  of  links  adjacent  to  constrained  nodes 
exceeds  the  maximum  number  of  700  as  set  by  the  variable  NCLIUK  in  the 
routine  RACAT. 

(j)  The  vehicle  payload  factor  (or  sortie  constraint 
factor)  is  0.  This  message  may  only  occur  when  the  V/TP  formulation 
is  used. 


d.  Inversion  Mode.  A  technique  used  in  standard  linear  programming 
codes  is  to  "invert"  the  basis  every  "n"  iterations  in  order  to  reduce 
the  size  of  the  eta  vector  file.  This  requires  the  elimination  of 
all  nonbasic  vectors,  i.e.,  those  which  entered  the  basis  at  gome 
previous  iteration  and  were  replaced  by  a  later  vector.  The  "n" 
iterations  >t  which  inversion  is  to  occur  is  specified  in  the  input 
data  file.  It  is  felt  that  experimentation  with  this  feature  will 
provide  the  analyst  knowledge  as  to  the  setting  of  this  parameter. 

Some  standard  LP  codes  set  it  at  75  or  compute  the  number  as  a  function 
of  significance.  In  E1MAM  the  analyst  has  the  ability  to  specify  it 
as  a  parameter.  The  formula  for  inversion  time  is  simply 

NNI  =  INVC  +  INN0 

where , 

INVC  =  The  inversion  iteration  parameter 

INN0  =  Iteration  at  which  the  inversion  occurred  and  if  the 
current  iteration  is  equal  to  NNI  then  an  inversion 
is  called  for. 

Figure  15  is  a  general  flow  chart  of  the  inversion  mode.  Subroutines 
S0RTAJ  and  REINV  are  required  to  control  this  processing. 

(l)  SORTAJ.  This  routine  controls  the  REINV  process.  First, 
the  basic  vectors  are  sorted  using  the  generalized  internal  "Shell 
Sort"  discussed  in  Chapter  2,  INPUT  MODULE.  The  array  in  which  the 
basis  vectors  are  stored  is  (NCHARC(u}>  1=1,  NRCW)  so  that  j 
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FIGURE  15.  GENERAL  FLOW  CHART  OF  INVERSION  MODE 


corresponds  to  a  row  of  the  basis.  The  vector  numbers  in  NCHARC  may 
be  classified  into  three  categories.  If  NCHARC(l)  =  0,  then  slack 
vector  I  is  on  row  I  (i.e.,  the  slack  vector  I  is  said  to  be  in  its 
home  position).  If  0  <  NCHARC(l)  £  NROW,  then  the  slack  vectoi  whose 
identifying  number  is  given  by  NCHARC(l)  has  been  entered  the  basis  on 
row  I.  If  NCHARC(l)  >  NRCW,  then  NCHARC (I)  gives  the  identifying 
number  of  a  structural  vector  which  has  entered  the  basis  on  row  I. 

A  count  is  made  of  the  number  of  structural  vectors  present  in  the 
basis  at  inversion  time  and  this  count  is  communicated  to  the  REINV 
routine.  Upon  the  completion  of  REINV,  the  SORTAJ  routine  resets  for 
normal  operation  and  returns  to  the  routine  from  which  it  has  been  called  — 
either  PHASE1  or  3MPLEX. 

(2)  REINV.  Figure  16  is  a  general  flow  chart  of  the  REINV 
routine.  The  vector  file  is  rewound  to  load  point  and  the  initialization 
of  the  basis  and  other  variables  is  performed  by  subroutine  SETUP. 

The  vector  file  is  then  read  to  obtain  a  structural  vector  that  is 
present  in  NCHARC,  i.e.,  a  basic  structural  vector  number  sorted  in 
SORTAJ.  Since  the  vector  file  has  been  written  sequentially  each  time 
a  vector  entered  the  basis,  the  vectors  may  be  located  one  at  a  time 
as  the  file  is  processed,  v'hen  an  active  vector  is  found  on  the  file, 
it  is  entered  into  the  tasio  via  the  SMPLEX  routine.  The  starting  basis 
for  the  invert  procedure  is  the  all  slack  basis.  The  basis  update  is 
performed  in  a  similar  manner  as  a  normal  iteration  in  the  maximization 
or  the  minimisation  operating  objective.  The  difference  is  that  the 
vector  entering  the  basis  is  allowed  to  replace  only  a  slack  vector 
which  was  not  in  the  basis  at  the  time  the  invert  procedure  was 
initiated.  A  new  eta  file  is  generated  representing  only  those  struc¬ 
tural  vectors.  Since  the  simplex  multipliers  are  not  required  to 
generate  a  new  vector  to  enter,  the  eta  file  does  not  have  to  be  read 
backward  and  the  shortest  chain  algorithm  does  not  have  to  be  used  to 
generate  a  chain  or  vector.  The  inversion  process  is,  therefore,  a 
relatively  rapid  operation. 

e.  Restart  Mode.  The  ability  to  restart  a  RACAT  operation  has 
many  advantages.  Two  important  uses  occur  when:  a  previous  run  has 
been  terminated  before  reaching  an  optimal  solution  (e.g.,  when  a  time 
limit  was  set  or  a  maximum  number  of  iterations  reached);  or  when  it  is 
desired  to  make  certain  small  changes  to  the  problem  and  obtain  a  new 
solution.  As  long  as  the  number  and  types  of  constraint  rows  have  not 
been  altered,  the  restart  mode  may  be  used  to  "crash"  in  previously 
used  vectors  without  computing  simplex  multipliers  and  generating  a 
new  vector  to  enter  the  basis  each  iteration. 

(l)  Figure  17  is  a  general  flow  chart  of  the  RESTART  routine. 

As  mentioned  in  the  discussion  on  "Input  Data  Requirements"  in  this 
chapter,  one  input  card  is  required  for  the  operation  of  the  RACAT 
Module.  A  card  with  "NORMAL"  punched  in  cols.  1-6  causes  the  RESTART 
routine  to  be  ignored  and  a  normal  solution  process  is  used.  The 
RESTART  operation  is  called  by  punching  the  word  "RESTART"  in  cols.  1-7. 
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FIGURE  17.  GENERAL  FLOW  CHART  OF  RESTART  ROUTINE 
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(2)  It  is  desirable  to  restart  a  job  as  near  to  the  cutoff 
point  as  possible.  The  save-basis  feature  of  ETNAM  saves  a  checkpoint 
at  desired  iteration  intervals.  The  frequency  of  saving  a  basis  is 
given  by  the  variable  ITSAVE  described  earlier  in  this  chapter.  The 
saved-basis  consists  of  the  numbers  of  the  structural  vectors  (or 
routes)  that  are  in  the  current  solution  and  a  bit  map,  NVTX,  which 
indicates  which  slack  vectors  are  in  the  solution.  This  is  called 
marking  the  slacks. 

(3)  Restarting  may  be  accomplished  in  two  ways:  by  marking 
the  slacks  and  without  marking  the  slacks.  Whenever  possible,  "marking 
the  slacks"  should  be  chosen  since  it  allows  the  program  to  continue 
from  the  point  where  the  last  basis  was  saved.  Method  2  is  less  effi¬ 
cient  and  should  only  be  used  when  it  is  not  possible  to  use  Method  1, 
i.e.,  when  the  save-basis  file,  TAPE9,  is  not  available  or  when  a  change 
in  the  data  has  been  made. 

(4)  To  restart  a  run  the  following  files  are  needed: 

(a)  TAFE1  -  the  input  file  from  the  run  to  be  restarted. 

If  some  solution  information  has  been  written  on  TAPE1,  it  should  be 
regenerated  by  the  Input  Module. 

(b)  TAPE2  -  the  edited  version  (from  the  VECTRS  program) 
of  TAPE2,  the  vector  file,  of  the  run  to  be  restarted. 

(c)  TAFE9  -  copied  from  TAPF  (the  save-basis  file  of  the 
run  to  be  restarted)  by  the  VECTRS  program. 

(d)  TAFE14  -  a  scratch  file  for  the  eta  file.  The  eta 
file  is  packed  into  core.  When  it  is  too  large  to  reside  completely 
in  core,  it  is  written  on  TAFE14.  For  the  CDC  6400  computer  this  must 
be  a  tape  file. 

(5)  The  card  data  to  be  input  is  as  follows:  RESTART  is 
punched  in  cc  1-7  of  the  RACAT  data  card.  The  other  parameters  on  this 
card  are  as  described  earlier  in  this  chapter.  However,  a  further  word 
of  clarification  is  needed  for  the  parameter  IREST.  When  IREST  =2, 
the  vector  numbers  and  the  bit  map,  NVTX,  must  be  input  on  cards  immed¬ 
iately  following  the  RACAT  data  card.  The  vector  numbers  are  read  in 
the  FORMAT  (I10,l4l5/(l6l5)).  The  integer  in  the  first  10  columns  of 
the  first  card  is  the  count  of  the  vector  numbers  to  be  read  in. 

(6)  The  bit  map  is  stored  in  the  first  20  locations  of  NVTX. 

For  the  CDC  6400,  NVTX  (l)  has  the  bits  numbered  60,  59>  ...»  3,  2,  1, 
counting  from  right  to  left  corresponding  to  the  first  60  slack  vectors, 
NVTX  (2)  has  the  bits  numbered  120,  119,  ...»  63,  62,  6l,  etc.  NVTX  is 
printed  in  OCTAL  format  at  invert  time,  at  the  end  of  Riase  I  and  at  the 
end  of  a  run.  The  vector  numbers  are  printed  when  NVTX  is  printed  and 
also  at  those  iterations  where  the  basis  is  saved.  NVTX  is  read  in  the 
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OCTAI.  format  (4520)  with  five  cards  being  required.  Those  words  of 
NVTX  (l)  to  NVTX  (20)  that  are  not  needed  by  the  bit  map,  have  zeros 
punched  on  the  data  cards.  For  the  I EM  System  360,  the  bits  are  numbered 
from  left  to  right  in  NVTX  and  NVTX  is  printed  in  hexa-decimal  format. 

(7)  A  separate  small  program  deck  called  VECTRS  is  used  to  edit 
the  vector  file  and  copy  the  save-basis  file,  of  a  run  which  is  to  be 
restarted.  This  procedure  is  optional  but  recommended  since  a  restart 
will  alter  the  contents  of  TAPE2  and  TAPE9.  For  a  VECTRS  computer  run, 
the  old  TAPES  (vector  file)  and  TAPE9  (save-basis  file)  are  input  as 
TAPE3  and  TAPE8  respectively.  Two  new  files  referred  to  as  TAPES  and 
TAPE9  are  generated.  These  become  the  TAPES  and  TAPE!9  needed  for  the 
restart  run.  There  is  one  parameter  data  card. 


cc  1-5,  NRES 
cc  6-10,  NCONST 
cc  11-15,  ICARD 


cc  16-20,  ISORTI 


Number  of  resources 

Number  of  node  constraints 

0,  Read  number  of  basis  vectors 
from  TAPE8,  i.e.,  old  TAPE9. 

1,  Read  numbers  of  basis  vectors 
from  cards  following  this  data 
card,  the  format  is  (110,  1415/ 
(1615)).  The  first  10  digit 
integer  is  the  count  of  the 
vectors . 

1,  if  V/TP  (sortie)  constraints  are 
used 

0,  otherwise 


(8)  A  limit  of  300  is  set  on  the  number  of  vectors  that  may  be 
crashed  into  the  basis.  The  vector  identification  numbers  are  stored 
in  the  array  NOVEC. 


(9)  The  restart  operation  is  similar  to  the  inversion  process 
in  that  the  vector  file  is  read  and  the  vector  number  checked  to  deter- 
min  whether  or  not  the  vector  is  to  be  entered  sequentially  since  the 
entire  array  NOVEC  is  searched  for  a  match  when  a  vector  record  is  read 
from  the  file.  If  the  vector  is  to  enter  it  is  structured  properly  and 
entered  via  the  SMPLEX  routine.  During  restart,  just  as  in  the  inversion 
mode,  the  eta  file  is  read  and  written  but  does  not  have  to  be  read 
backward.  This  restart  is  also  a  relatively  rapid  process. 

(10)  RESTART  checks  each  iteration  in  a  minimization  run  to 
determine  if  a  Riase  I  solution  has  been  reached.  If  this  occurs,  the 
restart  mode  is  temporarily  interrupted  and  a  Fhase  I  update  is  per¬ 
formed.  When  the  RESTART  routine  has  completed  "crashing  in"  all  the 
specified  vectors,  control  passes  to  the  appropriate  routine  for  a 
normal  continuation  toward  a  solution.  In  a  maximization  run,  MAXIO 
takes  control;  in  a  minimization  run,  either  HiASEl  or  MINIO  assumes 
operation  depending  on  the  phase  in  which  the  solution  process  is  to 
be  continued. 
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CHAPTER  4.  OUTPUT  MODULE 


1.  Purpose  and  Scope. 

a.  The  purpose  of  the  Output  Module  is  to  prepare  and  format  out¬ 
put  reports  based  on  the  results  of  the  RACAT  Module  solution.  The 
RACAT  results  are  a  series  of  mathematical  data  that  are  not  easily 
interpreted  in  their  raw  form ,  It  is  the  purpose  of  the  Output  Module 
to  prepare  and  format  these  basic  data  so  that  they  will  be  useful  and 
convenient  for  the  analyst. 

b.  The  ETNAM  Output  Module  has  been  designed  as  a  data  processing 
system  that  can  be  run  automatically  as  a  part  of  the  overall  ETNAM 
system,  or  can  be  run  separately  at  the  user's  option.  Two  files 
contain  the  output  data  from  the  RACAT  Module  and  are  used  as  input  to 
the  Output  Module.  If  these  two  files  are  saved, the  Output  Module  may 
then  be  rerun.  These  same  files  are  also  used  to  restart  the  RACAT 
Module  in  the  event  the  program  is  interrupted  or  in  the  event  that  a 
modification  to  the  original  problem  is  required. 

c.  The  purpose  of  this  chapter  is  bo  describe  the  Output  Module, 
its  inputs,  outputs,  and  the  methods  used  to  prepare  them.  The  formats 
of  the  ten  reports  that  may  be  produced  are  included  in  the  following. 

2.  Input  Data  Requirements.  The  Output  Module  has  as  its  inputs  the 
two  files  created  by  the  RACAT  Module,  the  data  file  and  the  vector 
file.  The  ETNAM  system  may  be  set  up  to  run  through  the  Input  Module, 
the  RACAT  Module,  and  the  Output  Module  in  one  computer  session,  or 
these  modules  may  be  run  separately. 

a.  Data  Files.  The  data  file  contains  much  of  the  data  prepared 
and  edited  in  the  Input  Module  as  well  as  the  solution  results  from  the 
RACAT  Module.  The  vector  file  contains  all  the  structural  vectors  that 
were  entered  into  the  "basis"  since  the  start  of  the  problem.  Each 
structural  vector  is  a  unique  mathematical  description  of  a  potential 
route  from  some  origin  to  its  related  destination.  This  is,  each  vector 
explicitly  describes  the  arcs  making  up  a  chain  from  an  origin  to  a 
related  destination,  the  resources  required  to  move  a  unit  of  the  com¬ 
modity  along  that  route,  and  identifies  the  constrained  noles.  Vehicle 
payload  factors,  if  the  arc  capacities  are  given  in  V/TP  are  also  pre¬ 
sent.  The  RACAT  algorithm  is  a  process  of  bringing  in  such  vectors  in 
an  iterative  procedure  until  uo  better  solution  is  possible.  In  the 
process  of  finding  the  solution  some  vectors,  brought  into  the  basis 

in  earlier  iterations,  may  be  replaced  by  other  vectors.  The  vectors 
which  remain  in  the  basis  are  referred  to  as  active  vectors.  Those 
which  have  been  replaced  are  inactive.  All  are  on  the  vector  file  after 
a  normal  RACAT  run  (i.e.,  no  restarts). 

b.  Data  Arrays.  The  main  program  of  the  Output  Module  reads  the 
data  file  and  loads  COMMON  with  the  data  it  contains  for  use  in  the 
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preparation  of  the  reports.  The  following  is  a  description  of  the  arrays 
in  COMMON  that  are  constructed  by  OUTGEN,  the  main  program. 

(1)  MATARC  (700,3) .  This  array  provides  basic  data  for  the 
arcs  of  the  network.  The  columns  of  the  array  contain,  respectively, 
the  FROM  node,  the  TO  node,  and  the  "mode"  of  the  arc.  It  should 

be  noted  that  arcs  are  not  necessarily  directed.  That  is,  flow  may  occur 
that  moves  from  the  TO  node  to  the  FROM  node.  (See  Chapter  2,  Inputs.) 
The  order  in  which  these  nodes  are  listed  is  the  same  as  that  used  in 
the  initial  input  form. 

(2)  RHS  (900) •  The  righ^-hand-side  array  contains  the  con¬ 
straint  values  input  to  the  Input  Module.  Each  element  of  the  array 
corresponds  to  a  row  of  the  basis.  For  example,  the  ith  element 

of  RHS  (where  i  Is  equal  to  or  less  than  the  total  number  of  arcs  in 
the  network)  will  contain  the  capacity  of  the  ith  arc. 

(3)  N0RDES  (150,4) .  This  array  contains  the  origin-destination 
pairs  and  information  about  them .  The  columns  of  the  array  contain  the 
following  information  respectively:  the  commodity  number,  the  origin 
node,  the  destination  node,  and  the  row  number  of  the  corresponding 
movement  requirement.  In  the  min- cost  and  the  min- time  modes  of  opera¬ 
tion,  movement  requirements  are  established  for  each  commodity.  These 
requirements  become  constraint  rows  of  the  basis  matrix  and  the  fourth 
column  of  the  NORDES  array  contains  the  row  number  corresponding 
movement  requirement.  In  the  min- cost  and  the  min- time  modes  of  opera¬ 
tion,  movement  requirements  are  established  for  each  commodity.  These 
requirements  become  constraint  rows  of  the  basis  matrix  and  the  fourth 
column  of  the  NORDES  array  contains  the  row  number  corresponding  to  the 
origin-destination  pair.  The  origin-destination  pairs  and  associated 
data  are  sequenced  by  commodity. 

(k)  NCD  (20) .  This  array  is  used  for  "bookkeeping"  within 
the  system.  It  indicates  the  row  of  NORDES  containing  the  last  origin- 
destination  pair  for  the  corresponding  commodity.  Up  to  twenty  commodi¬ 
ties  may  be  handled  within  the  ETNAM  system.  Thus  the  ith  element  of 
NCD  is  the  row  number  in  NORDES  of  the  last  entry  for  the  ith  commodity. 
The  origin-destination  pairs  corresponding  to  the  ith  commodity  are 
located  from  NORDES  (NCD(i-l)+l)  to  NORDES  ( NCD( i ) ) . 

(5)  SOLVEC  (900) •  The  solution  vector  array  contains  the 
values  of  the  solution  vector  at  the  completion  of  the  RACAT  algorithm. 
Thus  SOLVEC  (j)  contains  the  flow  over  the  route  represented  by  the 
activity  vector  occupying  the  jth  column  of  the  basis  at  solution  time. 
The  active  vectors  of  the  vector  file  thus  represent  routes  that  were 
included  in  the  solution.  The  flow  that  was  assigned  to  these  routes 
may  be  determined  by  finding  the  location  In  the  basis  of  these  routes 
and  referring  to  the  corresponding  element  of  SOLVEC. 


(6)  NCHARC  (3600).  This  array  contains  cross -indexing  infor¬ 
mation  that  permits  the  program  to  determine  which  vectors  on  the  vector 
file  are  active  vectors.  As  was  mentioned  above,  the  RACAT  algorithm 

is  an  iterative  procedure  that  selects  route-resource  combinations  that 
will  result  in  an  optimum  solution.  Thus  the  vector  file  contains  both 
active  and  inacti/e  vectors.  NCHARC  contains  an  ordered  set  of  active 
vector  numbers  and  is  used  to  "purge"  the  vector  file  of  inactive  vec¬ 
tors.  The  core  storage  occupied  by  NCHARC  during  the  initial  phase  of 
the  Output  Module  is  later  used  for  other  purposes,  hence  the  dimension 
of  NCHARC  is  set  to  provide  for  its  maximum  use. 

(7)  TITUS  (20).  This  arruy  contains  an  eighty  character  title 
supplied  by  the  user  and  accepted  in  the  Input  Module.  This  title  is 
printed  at  the  top  of  each  page  of  the  output  reports. 

(8)  NAMOD  (20).  This  array  contains  the  user  supplied  names 
of  the  modes. 

(9)  NACOM  (20).  This  array  contains  the  user  supplied  names 
of  the  commodities. 

(10)  NARES  (50 ).  This  array  contains  the  user  supplied  names 
of  the  resources. 

(11)  Parameters .  The  parameters  NROW,  NRES,  NARC,  NMOD,  NCOM, 
NCONST,  and  ISORTI  are  used  here  in  the  same  sense  as  described  in 
Chapters  2  and  3.  There  are  other  words  and  arrays  in  common,  some  of 
which  are  passed  from  the  RACAT  algorithm  to  the  Output  Module,  that 
are  not  indicated  in  the  above  list.  These  words  have  highly  specia¬ 
lized  purposes  having  to  do  with  the  specific  methods  of  programming 
the  output  reports. 

c.  Data  Cards.  In  addition  to  the  two  files  described  above,  one 
to  six  cards  are  read  which  indicate  the  report  options  to  be  followed 
in  the  Output  Module.  There  are  ten  reports  that  may  be  prepared  by 
the  Output  Module.  Tiree  of  these  are  always  prepared  but  the  remaining 
seven  may  be  suppressed  at  the  user's  option.  The  "Node  Workload 
Report"  summarizes  the  total  flow  by  commodity  that  passes  through 
selected  nodes  in  the  network.  The  user  must  specify  those  nodes  in 
the  network  to  be  reported. 

(l)  Card  Format  1.  The  first  field  (card  columns  (cc)  1-10 ) 
indicates  the  reports  that  the  user  does  not  want  prepared.  That  is, 
if  no  reports  are  indicated  on  this  card  all  reports  will  be  prepared. 

A  "1"  punched  in  the  ith  column  of  the  field  indicates  that  the  ith 
report  is  not  to  be  prepared.  (See  Table  1  for  a  list  of  the  reports 
and  options.)  The  second  field,  cc  11-20,  contains  the  number  of  the 
mode  (right  adjusted)  to  be  reported  in  the  Summary  of  Mode  X  Flow.  The 
total  number  of  nodes  to  be  included  in  the  Node  Workload  Report  is 
punched  (right  adjusted)  in  the  third  field,  cc  21-30. 
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TABLE  1.  OUTPUT  REPORT  OPTIONS 

1.  Run  Parameter 

2.  Summary  of  Network  Plow;  Ton-Miles  Flow  by  Mode* 

3.  Summary  of  Mode  X  Flow  -  Optional 

4.  Movement  Constraint  Report  -  Optional 

5 •  Route  Report 

6.  Origin  Workload  Report  -  Optional 

7*  Destination  Workload  Report  -  Optional 

8.  Node  Workload  Report  -  Optional*# 

9.  Summary  of  Commodity  Flow  -  Optional 

10.  Summary  of  Resources  Used  -  Optional 

♦Reported  generated  only  when  the  arc  capacities  are  expressed  in 
tons. 

**User  specifies  the  nodes  for  which  data  is  accumulated. 

(2)  Card  Format  2.  A  second  card  format  is  used  to  input  the 
node  names  of  those  nodes  to  be  reported,  as  specified  in  the  third 
field  above.  Ten  node  names  per  card  are  allowed  with  each  name  left 
adjusted  within  a  field  of  8  card  columns.  The  number  of  cards  used 
must  be  consistent  with  the  total  number  of  nodes  as  indicated  in  the 
first  card  type.  The  fields  are  cc  1-8,  cc  9-16,  ...,  cc  73-80. 

3.  Outputs .  As  indicated  previously,  ten  reports  (Tfctble  l)  may  be 
prepared  using  the  Output  Module,  Some  of  these  reports  are  so  basic 
in  their  content  that  they  are  automatically  prepared.  Other,  more 
specialized  reports  are  prepared  at  the  user's  option.  Since  all  these 
reports  are  prepared  from  the  data  and  vector  files,  the  Output  Module 
may  be  run  with  alternate  report  options  after  the  initial  running  of 
the  full  system.  A  brief  explanation  of  each  of  the  ten  reports  is 
given  below.  It  is  suggested  that  these  explanations  be  read  with  the 
sample  reports  shown  in  the  SD,  Volume  I,  User's  Jfenual,  Enclosure  1. 

a.  The  Run  Parameter  Report.  The  Run  Parameter  Report  is  a  brief 
summary  of  the  most  basic  parameters  used  in  the  run  and  the  value  of 
the  objective  function  at  solution.  The  number  of  substitutions,  the 
number  of  modes,  the  number  of  commodities,  the  number  of  resources, 
the  number  of  constrained  nodes  and  the  type  or  arc  constraints  used 
are  shown.  The  number  of  iterations  used  by  the  RACAT  Module  in  solving 
the  network  and  the  value  of  the  objective  function  are  also  shown. 

The  value  of  the  objective  function  is  the  total  flow  achieved,  the 
total  cost  required,  or  the  total  flow  unit-hours  (ton-hours  if  capa¬ 
cities  are  in  tons)  required  depending  on  the  run  objective  --  max-flow, 
min-cost,  or  min-time  respectively.  In  addition,  the  report  options 
requested  by  the  user  are  presented. 

b.  Summary  of  Network  Flow  and  Ton- Mile  Flew  by  Mode.  The  Summary 
of  Network  Flow  Report  shows  the  total  flow  over  each  of  the  arcs  in 
the  network  by  commodity.  Thus  for  any  given  arc  in  the  network  the 
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user  can  find  the  total  flow  by  commodity  over  that  arc  that  results 
from  all  chains  using  it.  The  total  flow,  the  capacity,  and  the  re¬ 
sidual  are  shown.  The  flow  by  commodity  is  also  detailed.  Only  those 
commodities  that  have  a  flow  s.5  are  shown  in  the  section  of  the  report 
labeled,  "Plow  by  Commodity...".  Those  arcs  with  zero  flow,  that  is 
those  arcs  that  were  not  used  in  the  solution,  are  listed  separately. 
Only  a  total  flow  i.01  is  shown. 

(l)  Two  potential  areas  of  misunderstanding  should  be  noted. 

As  has  been  mentioned  previously,  users  have  the  option  of  making  any 
of  the  arcs  in  the  network  "directed"  or  "undirected."  That  is,  the 
user  may  specify  whether  flow  is  to  be  permitted  in  only  the  indicated 
direction  or  in  both  directions.  If  the  arc  is  designated  as  an 
"undirected"  arc  the  RACAT  Module  may  assign  flow  to  move  from  the  TO 
node  to  the  FROM  node.  Further,  it  is  possible  that  a  single  arc  may 
have  flow  in  one  direction  for  one  commodity  and  in  the  opposite  direc¬ 
tion  for  another.  This  is  quite  normal  in  a  multi -commodity  problem. 
Thus  it  is  not  always  possible  to  give  precise  meaning  to  the  FROM  and 
TO  nodes  of  a  link.  In  this  report  the  FRCM  and  TO  nodes  are  listed 
exactly  as  they  are  designated  on  the  input  forms.  The  actual  direction 
of  any  flow  is  easily  determined  from  the  chains  shown  in  the  Route 
Report . 


(2)  A  second  area  of  possible  misunderstanding  has  to  do  with 
the  commodity  flows.  The  "Total  Flow,"  "Capacity,"  and  "Residual"  are 
in  equivalent  tons  or  v/TP,  while  the  flows  by  commodity  are  3hown  in 
their  normal  units  (people,  barrels,  etc.).  Thus  the  flows  may  net 
sum  directly  to  the  total  flow.  Additionally,  the  flows  are  shown  as 
integers  and  may  not  exactly  equal  the  total  flow  as  indicated  even 
when  all  flows  are  in  tons  because  of  rounding  error. 

(3)  While  the  total  flow  by  arcs  is  being  accumulated,  the  flow 
by  modes  multiplied  by  the  distance  of  the  arc  is  also  accumulated.  The 
ton-mile  flow  and  the  percent  of  the  total  ton-mile  flow  for  each  mode 
are  shown  in  the  Summary  of  Ton-Mile  Flow  by  Mode  Report.  This  report 
is  omitted  if  v/TP  are  used  instead  of  tons  for  the  capacities  of  the 
arcs. 


u .  Summary  of  Mode  X  Flow.  This  report  is  exactly  like  the  Summary 
o'1  Network  Flow  Report  except  that  it  shows  only  those  arcs  of  the  mode 
designated  by  the  user.  Hie  basic  idea  behind  this  >-eport  is  that  the 
user  may  be  especially  interested  in  the  flows  over  arcs  of  a  particular 
mode.  For  example,  analysts  are  often  especially  interested  in  the 
activity  on  transfer  arcs.  This  report  gives  the  user  the  ability  to 
obtain  a  special  report  detailing  a  particular  mode. 

d.  Movement  Constraint  Report.  This  report  shows  the  factors  that 
constrain  the  solution  as  found  by  the  RACAT  Module.  For  example,  in  a 
max-flow  run  certain  arcs  and  certain  resources  may  be  the  limiting 
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factors  on  the  total  flow  achieved.  The  movement  constraint  report 
will  list  the  arcs,  resources,  nodes  and  commodity  movement  requirements 
which  are  constraining  and  will  shew  their  corresponding  shadow  prices. 
The  shadow  price  of  a  constraint  row  is  defined  as  the  change  in  the 
objective  function  per  unit  increase  in  the  corresponding  right-hand- 
side  value,  assuming  that  a  change  of  basis  is  not  required  to  maintain 
feasibility.  In  other  words,  it  is  the  rate  of  change  of  the  objective 
function  with  respect  to  the  corresponding  right-hand- side  holding  the 
other  right-hand- side  values  fixed. 

e.  The  Route  Report.  The  Route  Report  is  a  detailed  listing  of 
all  the  chains  used  in  the  solution,  together  with  the  flow  assigned 
to  that  route  and  the  resources  used  to  accomplish  it.  The  node  names 
of  any  saturated  constrained  nodes  will  also  be  printed.  This  report 

is  the  most  basic  report  of  the  series  and  may  often  be  the  most  useful. 
The  report  shows  how  the  optimum  flow  meeting  the  requirements  can  be 
achieved  through  the  actual  selection  of  routes  and  the  assignment  of 
resources  to  that  route.  Only  routes  with  flow  £.01  are  shown. 

(1)  Route  Number.  The  route  number  is  assigned  sequentially 
simply  for  identification.  Routes  are  listed  in  order  by  commodity. 

The  route  flow,  cumulative  flow,  and  travel  time  for  the  route  are 
shown.  The  route  nodes  in  order  of  their  occurrence  are  then  shown. 

This,  of  course,  permits  the  user  to  trace  the  exact  route  used.  The 
resources  assigned  to  the  route  are  shown  for  each  route  and  for  each 
commodity  in  total. 

(2)  Route  label.  Each  route  is  assigned  a  label  which  in  most 
cases  will  be  unique.  If  two  or  more  runs  are  to  be  compared  the  label 
will  tell  which  routes  are  duplicates  even  if  the  route  numbers  are 
different . 

f .  Origin  Workload  Report.  The  Origin  Workload  Report  shows  the 
flow  by  commodity  out  of  each  origin.  Each  node  designated  as  an 
origin  by  the  user  in  the  input  data  is  printed  with  its  corresponding 
flow.  If  a  node  is  both  an  origin  and  an  intermediate  node  (on  separate 
chains)  only  the  originating  flow  will  be  reported.  Flows  are  shown  in 
thert  normal  units  (PAX,  barrels,  short  tons,  °tc.). 

g.  Destination  Workload  Report.  The  Destination  Workload  Report 
is  identical  to  the  origin  workload  report  except  that  the  terminating 
nodes  of  chains  are  reported. 

h.  Node  Workload  Report.  The  Node  Workload  Report  is  identical  in 
format  to  the  Origin  Workload  Report  and  the  Destination  Workload  Report. 
The  user  may  request  that  up  to  fifty  nodes  be  included  in  the  list.  In 
certain  studies,  particularly  vulnerability  studies,  it  may  be  useful  to 
know  the  flow  by  commodity  passing  through  a  particular  node. 


i.  Summary  of  Comnodlty  Flow.  The  Summary  of  Commodity  Flow  Report 
summarizes  the  flow  over  the  network  by  commodity.  In  min- cost  and 
min-time  runs  these  flews  will  correspond  to  the  movement  requirements. 
In  max-flow  runs  these  flows  will  be  the  maximum  flow  achieved,  given 
the  constraints  of  the  network  and  the  available  resources. 

j.  Summary  of  Resources  Used.  The  Summary  of  Resources  Used  Report 
shows  the  total  resources  assigned  to  the  network  in  order  to  accomplish 
the  flow,  the  residual  (unassigned)  resources,  and  the  routes  to  which 
these  resources  are  assigned. 

4.  Methods . 

a.  The  Output  Module  actually  consists  of  a  main  program  called 
OUTGEN  and  eight  subroutines.  OUTGEN  reads  the  data  file,  TAPE1,  and 
loads  COMMON  as  was  described  earlier.  OUTGEN  then  prepares  the  Run 
Parameter  Report  and  calls  each  of  the  tabulating  and  report  writing 
reports  in  their  proper  sequence.  Figure  18  is  a  schematic  diagram  of 
the  relationship  between  the  main  program  and  the  subroutines  and  their 
corresponding  reports. 

b.  The  JOTS  subroutine  reads  the  vector  file,  TAFE2,  purges  that 
tape  of  inactive  vectors,  and  tabulates  arc  flows  by  commodity  for  the 
use  of  the  SUMRY  and  XTRANX  subroutines.  The  tabulation  of  flows  by 
arc  by  commodity  requires  a  large  data  array  in  memory.  The  actual 
space  required  is  equal  to  the  number  of  commodities  times  the  number 
of  arcs  in  the  network.  Thus  if  the  maximum  values  of  these  variables 
are  used,  a  table  of  14,000  cells  would  be  required  for  the  700  arc 
mode.  JOTS  prepares  the  table  and  SUMRY  and  XTRANX  use  the  data  to 
prepare  the  Summary  of  Network  Flows  and  the  Summary  of  Mode  X  Flows 
Reports  respectively.  Flows  by  commodity  are  shown  in  the  same  units 
of  measure  as  inpu  (e.g.,  PAX,  barrels,  tons). 

c.  JOTS,  SUMRY,  and  XTRANX  are  programmed  to  permit  a  variable 
data  array.  If  enough  core  storage  is  available  to  accumulate  all  the 
data  in  one  pass,  these  subroutines  handle  the  problem  in  this  manner. 

If  the  network  size  is  such  that  multiple  passes  are  required,  each  sub¬ 
routine  is  set  up  to  handle  the  data  in  this  fashion.  A  flow  chart  of 
the  JOTS  subroutine  Is  shown  in  Figure  19 . 

d.  The  SUMRY  and  XTRANX  subroutines  are  straight-forward  data  pro¬ 
cessing  routines  that  simply  format  the  data  and  produce  the  reports. 
SUMRY  prints  the  data  as  accumulated  by  JOTS.  XTRANX  scans  the  data 
prepared  by  the  JOTS  routine  printing  only  those  arce  of  the  designated 
mode. 


e.  MOVCON  is  a  straight-forward  data  processing  routine  that  formats 
and  prints  the  movement  constraints.  The  shadow  prices  we  simply  the 
simplex  multipliers  from  the  RACAT  Module  and  are  contained  in  the 
COMMON  array  ALPHA. 
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FIGURE  18.  OUTPUT  MODULE  SCHEMATIC 
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FIGURE  19.  (continued) 
(Sheet  2  of  2) 
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f.  ROUTES  is  a  somewhat  more  complex  routine.  Its  basic  information 
is  obtained  from  the  purged  vector  file  which  contains,  in  a  raw  state, 
all  the  data  required.  Routes  must,  however,  determine  the  nodes  con¬ 
tained  in  the  chain  and  place  them  in  their  proper  sequential  order. 

The  vector  file  record  contains  arc  numbers  and  their  order  is  not 
necessarily  correct. 

g.  LODOUT  and  NODLOD  prepare  the  Origin  Workload  Report,  the  Des¬ 
tination  Workload  Report,  and  the  Node  Workload  Report.  In  preparing 
these  reports  it  is  essential  that  the  data  be  carefully  handled  in 
order  to  avoid  possible  double  counting.  For  example,  an  origin  node 
may  have  traffic  that  goes  through  the  node  (coming  from  other  origins) 
as  well  as  traffic  originating  at  the  node.  LODOUT  and  NODLOD  handle 
this  problem  by  examining  and  tabulating  data  from  chains  rather  tlian 
from  arcs.  Thus  both  LODOUT  and  NODLOD  subroutines  are  called  by  the 
ROUTES  subroutine  when  a  completed  chain  has  been  constructed. 

h.  SUMM  is  a  relatively  straight-forward  data  processing  subroutine 
that  prints  out  the  Summary  of  Commodity  Flow  and  the  Summary  of  Resources 
Used  Reports. 
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CHAPTER  5.  POST  OPTIMAL  PACKAGE  (POP)  MODULE 

1.  Purpose  and  Scope,  The  POP  Module  provides  the  analyst  with  a  tool 
for  sensitivity  analysis  of  the  right-hand-side  values.  This  is  known 
as  right-hand-side  (RHS)  ranging.  RHS  ranging  is  the  only  post  optimal 
procedure  possible  for  the  current  structure  of  the  ETNAM  system.  RHS 
ranging  gives  the  range  of  values  a  RHS  element  may  assume,  other  RHS 
elements  being  held  fixed  at  their  given  value,  for  which  the  current 
solution  remains  optimal.  If  the  RHS  element  were  to  move  just  beyond 
the  range,  one  of  the  vectors  in  the  basis  would  have  to  be  replaced  by 
a  vector  currently  not  in  the  basis  for  the  solution  to  remain  optimal. 

If  such  a  vector  to  enter  is  not  found,  then  the  solution  would  be 
infeasible  for  the  RHS  element  beyond  that  end  point  of  its  range. 
Enclosure  4  gives  a  mathematical  treatment  of  the  applicability  or 
inapplicability  of  the  various  post  optimal  procedures  to  ETNAM. 

2.  Input  Data  Requirements .  Two  types  of  data  are  required  to  operate 
the  POP  Module:  the  eta  file  and  some  parameters  and  data  arrays  writ¬ 
ten  on  TAPE14  by  RACAT,  and  card  input. 

a.  TAPE14  Input.  This  is  a  binary  file  and  must  be  a  magnetic  tape 
in  the  Control  Data  6400  version  since  it  is  accessed  by  LTRIO — a  Control 
Data  system  subroutine — which  can  only  access  magnetic  tape  drives.  The 
contents  of  TAPE14  are  given  below.  On  the  IBM  System/ 360,  this  may  be 

a  disk  file. 

(1)  The  Eta  File.  In  the  RACAT  Module,  the  eta  buffers  and 
file  records,  if  any,  are  all  written  on  TAPE14.  This  is  the  complete 
eta  file  from  which  the  basis  inverse  is  computed. 

(2)  Parameters .  These  are:  the  array  TITLE,  NARC,  MRES, 

NCONST,  KETA  (the  number  of  eta  vectors),  ISQRTI,  NCQM  and  the  array 
CONFX.  The  definitions  for  them  are  given  in  Chapters  2  and  3. 

(3)  Right -Hand -Side.  The  RHS  elements  sire  given. 

(4)  HCHARC.  This  is  the  array  identifying  the  vector  loca¬ 
tions  in  the  basis. 

(5)  SOLVEC.  This  is  the  solution  vector. 

(6)  NGRDES.  This  is  the  data  array  relating  to  the  origin- 
destination  pairs. 

b.  Data  Cards.  These  data  cards  are  used  to  specify  an  opera¬ 
tional  parameter  NRCW  and  the  selection  of  RHSs  for  which  a  range  is 
desired.  The  RHS  selections  may  be  made  by  key  words  or  by  specifying 
the  lover  and  upper  limits  of  a  range  of  row  numbers.  There  may  be  as 
many  of  these  pairs  of  limits  as  desired.  The  data  cards  are  followed 
by  an  END  card.  Figure  20  gives  a  schematic  of  the  data  card  structure. 
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(1)  Operational  Parameter.  This  is  the  parameter  NROW.  Its 
value  may  be  obtained  from  the  first  Input  Module  report.  It  is  the 
number  of  rows  in  the  problem.  This  value  is  punched  as  an  Integer, 
right  adjusted,  in  card  columns  (cc)  1-5  on  the  first  card.  It  is  the 
only  value  on  the  first  card  and  must  be  present  for  the  eta  file  to  be 
read  correctly. 

(2)  Key  Words.  The  key  words  are  punched  left  adjusted  in 
cc  1-4.  Bach  key  word  is  alone  on  a  card  unless  it  is  the  key  word 
SOME,  which  is  explained  below.  Each  of  the  key  word  cards  is  optional. 
They  may  be  in  any  order. 


(a) 

ALL. 

This  ranges  each  of  the  RHSs. 

(b) 

ARCS. 

This  ranges  each  of  the  RHSs 

associated  with 

the  arcs. 

(c) 

NODE. 

This  ranges  each  of  the  RHSs 

associated  with 

the  constrained  1 

nodes . 

(d)  RES.  This  ranges  each  of  the  RHSs  associated  with  the 
resource  inventories. 

(e)  REQ.  This  ranges  each  of  the  RHSs  associated  with 
the  movement  requirements. 

(f)  SOME.  Ch  each  SOME  card  at  least  one  set  of  limits 
must  be  specified.  There  may  be  up  to  seven  pairs  of  limits  on  each 
card  and  as  many  cards  as  desired.  If  the  lower  limit  of  a  given  pair 
is  n  and  the  upper  limit  m  then  the  RHSs  corresponding  to  rows  n  through 
m  are  ranged,  n  s  m.  A  SOME  card,  as  has  been  noted,  must  have  the 
word  SOME  In  cc  1-4.  Card  Columns  11-20,  21-30,  31-40,  4l-50,  51-60, 
61-70,  71-80  are  the  fields  for  pairs  of  limits.  A  field  may  be  leit 
empty.  The  lower  limit  is  right  adjusted  in  the  first  five  columns  of 

a  field  (e.g.,  cc  11-15)  and  the  upper  limit  is  right  adjusted  In  the 
second  half  of  the  field  (e.g.,  cc  16-20).  Other  key  word  cards  may 
follow  a  SOME  card. 

3.  Output  Reports.  POP  generates  two  reports. 

a.  The  first  report  gives  the  title  and  identifies  the  row  numbers 
associated  with  each  type  of  constraint,  i.e.,  arc,  resource,  node,  and 
movement  requirements. 

b.  The  second  report  is  a  sunmary  by  row  showing  the  RHS  value; 
the  total  flow  on  arcs  or  constrained  nodes,  or  the  amount  of  resource 
used,  or  the  commodity  requirements  moved  (depending  on  the  type  of 
row);  and  the  lower  and  upper  value  of  that  RHS  range. 
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k.  Methods .  The  Post  Optimal  Package  consists  of  five  routines,  the 
main  program  POP  and  four  subroutines,  ETAS,  RANGE,  WRITE  and  ERROR. 

Each  routine  is  discussed  belov  and  a  block  diagram  is  given  in 
Figure  21. 

a.  POP.  This  is  the  main  program  of  the  Post  Optimal  Package.  It 
controls  the  operational  flow  of  the  whole  program.  It  calls  EEAS  to 
initialize  the  data  arrays  and  the  eta  file.  It  readB  the  input  data 
and  checks  for  errors  and  calls  ERROR  if  necessary.  It  calls  RANGE  to 
find  the  RHS  range  and  it  calls  WRITE  to  generate  the  reports. 

b.  ETAS.  ETAS  is  initially  called  by  POP  to  read  TAPE14  for  the 
parameters,  the  eta  file,  the  RHS  vector  and  other  necessary  data  as 
described  in  paragraph  2. 

(1)  H^b  is  computed  where  b  is  the  RHS  vector  and  the  basis 
inverse  B-1  is  a  product  of  elementary  column  matrices  */.  Hie  column 
of  xf  vhich  differs  from  the  corresponding  identity  matrix  is  known  as 
the  eta  (xt )  vector,  (kily  the  non-zero  elements  of  Xj  with  proper 
identification  are  stored  in  the  eta  file.  The  identification  consists 
of  an  index  for  each  xt  denoting  the  column  of  x/,  and  a  bit  map  indi¬ 
cating  the  position  in  xj  of  the  zero  and  non-zero  elements  for  each  eta 
vector  i . 

(2)  b  =  col.  (bj ,  bg ,  ...  ,  b^. ,  ...  ,  bB ) .  If  the  range  of  b^ 
is  to  be  fouid,  the  coefficient  of  1^  in  the  product  B_1b  needs  to  be 
found.  This  is  done  by  separating  b  as  b  =  col.  (bt ,  h,,  ...  ,  b^, 

^t+i>  >  \)  +  co1*  (0,  0,  . . .  ,  \  ,  0,  . . .  ,  0)  and  then  replac¬ 

ing  \  by  1.  Then  the  following  product  is  computed  x«h,,'_x  •••  Xb'ki 
[col.  (b,  ,  bg ,  ...  ,  ,  0,  ...  ,  b^  )  +  col.  (u,  0,  ...  ,  1,  0, 

...  ,  0) J  where  q  is  the  number  of  eta  vectors,  m  the  number  of  rows 
and  the  1  is  in  the  ktft  position  of  the  last  vector.  ThiB  product  must 
be  a  0. 


c.  RANGE.  If  there  are  m  rows  in  the  basis,  the  inequality  Flb  i  0 
yields  m  separate  inequalities.  The  range  of  \  is  essentially  the  set 
of  numbers  which  satisfy  the  m  inequalixies .  The  upper  limit  of  the 
range  is  the  least  upper  bound  of  the  limits  given  by  the  m  inequalities 
and  the  lower  limit  of  the  range  is  the  greatest  lower  bound  of  those 
m  limitt-. 

WRITE.  This  subroutine  generates  the  reports  discussed  in 
paragraph  3.” 

e*  «  This  subroutine  gives  error  messages  for  errors  gener¬ 

ated  by  card  input  and  by  file  input.  The  computer  run  is  terminated 
where  appropriate. 


FIGURE  21.  POP  MODULE  FUNCTIONAL  FLOW  CHART 

(Shoot  1  of  2) 


FIGURE  21. 

(Sn«et 


(l)  File  Input  Error.  When  reading  TAPE14,  there  are  three 
error  checks,  each  of  which  will  terminate  the  run. 

(a)  If  the  number  of  rows  of  the  problem  (NROW)  is  not 
equal  to  the  value  specified  on  the  card  input,  the  run  is  terminated. 
NROW  must  be  input  correctly  by  card  to  read  IAPE14  properly.  But 
TAPEL4  itself  contains  the  parameter  NROW  that  was  used  in  the  RACAT 
run.  These  two  are  compared  and  must  be  equal. 

(b)  If  the  number  of  eta  records  differ  from  those  used 
in  the  RACAT  run,  the  POP  run  is  terminated.  As  TAPEL4  is  read  for  the 
eta  file,  the  number  of  eta  records  is  tabulated  and  compared  with  the 
parameter  on  TAPED.4  which  gives  the  number  of  eta  records  with  which 
RACAT  ended. 


(c)  Tho  run  is  terminated  if  there  is  a  parity  error  on 
TAPEL4.  If  the  parity  error  persists,  the  RACAT  run  may  be  restarted 
from  the  optimal  basis  and  TAPEL4  will  be  regenerated. 

(2)  Card  Input  Error.  The  computer  run  is  terminated  if  a 
key  word  is  misspelled.  If  on  a  SOME  card  a  lower  limit  or  the  corre¬ 
sponding  upper  limit  ie  zero,  an  error  message  is  generated  and  that 
limit  pair  is  not  used.  In  this  case  the  computer  run  is  not  terminated . 


CHAPTER  6.  ACCURACY 


1.  Purpose  and  Scope.  Whenever  a  computer  program — like  tne  ETNAM 
computer  system— is  written  which  involves  data  inputs  and  many  numer¬ 
ical  computations,  the  analyst  must  always  reckon  with  the  possible 
inaccuracy  of  the  results.  This  depends  on  the  accuracy  of  the  input 
data,  the  number  of  bits  per  computer  word  (which  determines  the 
significant  digits  in  the  computation),  and  the  tolerance  used  in  the 
program  to  decide  when  two  numbers  are  considered  equal.  These  aspects 
will  be  discussed  in  this  chapter  and  an  attempt  made  to  draw  the 
analysts'  attention  to  areas  most  susceptible  to  numerical  problems. 

2.  Inputs.  The  preparation  of  the  input  data  consists  of  three  dis¬ 
tinct  activities.  Errors  are  easily  made  in  all  three  of  these  activ¬ 
ities  . 

a.  Data  Collection.  Collection  of  the  numerical  data  (like  the 
resource  prices,  productivities,  and  inventories;  the  network  link 
capacities  and  distances;  and  vehicle  speeds )  is  a  potential  source 

?  errors.  The  resource  productivities  are  usually  numbers  in  the 
xange  10"3  to  1CT8  .  It  is  difficult  to  recognize  errors  in  such  numbers 
since  they  are  not  intuitively  related  to  operating  characteristics. 
Unfortunately  even  a  small  error  in  productivity  data  may  result  in  a 
solution  quite  different  from  the  optimum.  For  example,  if  a  POL  route 
using  waterway  in  the  optimum  solution  is  closely  competitive  with  a 
POL  route  using  rail,  even  a  small  error  in  favor  of  rail  will  result 
in  a  solution  different  from  the  optimum.  It  is  difficult  to  know  in 
advance  whether  a  particular  resource  will  be  sensitive.  It  is  often 
possible  to  determine  from  the  Movement  Constraint  Report  which  items 
of  the  input  data  are  the  most  sensitive.  Thus  it  is  worthwhile  to 
recheck  these  sensitive  data  after  a  solution. 

b.  Network  Construction.  A  transportation  and  activity  network 
consists  of  various  kinds  of  links  (e.g.,  a  road  between  two  cities 
could  be  a  link  in  the  network  as  could  the  transfer  activity  of  off¬ 
loading  a  ship  onto  trucks). 

(1)  Care  must  be  taken  in  constructing  the  network  to  include 
all  the  necessary  links.  For  a  transportation  network,  separate  net¬ 
works  are  usually  constructed  for  each  transportation  mode  and  are  then 
"hooked"  together  with  links  representing  transfer  operations.  The 
network  for  each  transportation  mode  may  be  drawn  on  paper  and  checked 
quite  carefully.  The  transfer  links  ere  more  difficult  to  check  since 
the  FROM  node  is  on  the  network  of  one  transportation  mode  and  the  TO 
node  is  on  another. 

(2)  Most  links  may  be  used  in  either  direction,  but  to  model 
the  actual  activity  correctly  it  may  be  essential  that  certain  links 
be  directed.  Consider  Figure  22,  for  example. 
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FIGURE  22.  EXAMPLE  OF  NEED  FOR  DIRECTED  LINKS 

Node  1  is  the  origin,  say  a  ship,  of  a  commodity.  Another  commodity  is 
already  being  shipped  on  highway  to  Node  2  where  it  is  transferred  to 
rail  via  Link  3,  using  it  to  capacity.  Link  1  and  Link  2  are  transfer 
links  from  the  ship  to  highway  and  railroad  respectively.  More  of  the 
second  commodity  could  be  transferred  from  Node  2  to  3  than  the  capacity 
of  Link  3  by  moving  from  highway  to  ship  to  rail.  A  "nonsense"  route 
could  thus  be  created.  In  order  to  prevent  this,  Links  1  and  2  should 
be  directed  from  Node  1  to  2  and  from  Node  1  to  3  respectively. 

(3)  In  structuring  a  network  it  is  possible  to  overconstrain 
it.  That  is,  not  enough  throughput  capacity  1b  available  to  make  all 
the  delivery  requirements.  As  a  result,  an  infeasible  solution  ensues. 
This  happens  frequently  on  newly  structured  networks.  To  circumvent 
the  problem,  a  dummy  link  is  inserted  into  the  network  for  each  origin- 
destination  pair.  But  a  very  high  toll  is  assigned  to  it.  This  means 
that  the  algorithm  will  move  as  much  of  the  requirements  as  possible 
over  the  normal  network  links  and  use  the  dummy  or  escape  links  only  as 
a  last  resort.  In  this  way  the  analyst  gets  an  optimal  solution  and  is 
also  told  how  much  of  the  delivery  requirement  vaB  not  met. 

c.  Preparation  of  the  Data  in  Punched  Card  Form.  Clerical  and 
keypunch  errors  can  result  in  wasted  computer  time.  To  facilitate  the 
data  checking,  the  recommended  procedure  is  to  run  the  Input  Module 
first  and  study  the  generated  output  reports  of  the  input  data.  If  the 
network  is  a  newly  structured  one,  it  may  take  the  first  RACAT  run  to 
fully  debug  it. 

3.  Calculations.  When  doing  numerical  computations,  the  problems 
associated  with  truncation  and  tolerance — usually  referred  to  as  numer¬ 
ical  difficulties — must  be  considered.  Hie  sixty  bit  word  of  the 
Control  Data  6h00  and  double  precision  on  the  IEM  Systenv/360  reduce 
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this  problem,  but  loss  of  significance  may  still  occur.  The  ETNAM 
computer  programs  have  been  written  incorporating  a  number  of  features 
designed  to  reduce  the  effect  of  the  numerical  problems. 

a.  Re -Inversion.  An  example  of  such  a  feature  is  the  re- inversion 
technique.  The  inverse  of  the  basis  matrix  is  used  repeatedly  in  the 
algorithm  and  must  therefore  be  available  to  the  program.  Since  it  1j 
too  large  to  be  maintained  in  core  storage,  the  basis  inverse  is  kept 
in  the  "product  form"  on  the  eta  file.  Each  iteration  of  the  algorithm 
generates  an  "eta  vector"  and  the  product  of  all  the  eta  vectors  on  the 
file  is  the  basis  inverse. 

(1)  After  100  iterations  there  are  100  eta  vectors  on  the  file 
and  100  vector  multiplications  must  be  performed  to  generate  the  data 
necessary  for  the  next  iteration.  This  is  time  consuming  and  it 
increases  the  likelihood  of  numerical  difficulties  because  of  the 
repetitive  multiplication.  It  is  typical  of  these  problems  that  routes 
(vectors)  entered  into  the  basis  in  the  early  stages  may  not  remain 
through  the  later  stages.  Thus  it  is  quite  probable  that  after  100 
iterations,  a  far  lesser  number  (e.g.,  40)  of  the  routes  may  actually 
be  active.  Ry  a  process  called  re-inversion  we  can  reduce  the  length 
of  the  eta  file  to  precisely  the  number  of  active  routes  in  the  current 
basis.  Thus  periodic  re-inversicn  is  usually  worthwhile. 

(2)  There  is  a  trade-off  involved.  A  shorter  eta  file  will, 
of  course,  reduce  the  time  required  per  iteration,  but  re- inversion 
itself  requires  time.  In  the  ETNAM  system,  the  user  may  specify  how 
many  iterations  are  to  be  executed  before  re-inversion  takes  place. 

b.  "Tie "  Routes .  Routes  consist  of  a  sequence  of  connected  links 
of  the  network  and  the  resource  required  to  sustain  a  flow  over  them. 

The  desirability  of  including  any  route  in  the  basis  is  determined  by 
the  sum  of  the  costs  of  the  links  of  the  route.  The  cost  of  each  link 
is  in  terms  of  its  toll,  the  simplex  multipliers  (shadow  prices)  of  the 
link  and  of  the  adjacent  constrained  nodes,  and  of  the  resources  used 
on  the  link.  Sometimes  two  or  more  routes  may  have  the  same  value, 
which  means  that  the  routes  are  "tied"  for  selection  as  the  "shortest" 
route.  Ties  are  broken  by  choosing  the  route  with  the  least  distance. 

If  a  secondary  test  results  in  a  tie  too,  then  the  first  route  found 
is  used.  "Tie"  routes  are  not  uncommon.  In  the  maximum  flow  problem 
the  objective  Is  to  maximize  throughput  without  regard  to  time  or  cost 
and  hence  maximum  capacity  routes  rather  than  minimum  cost  or  time 
routes  are  chosen.  These  maximum  capacity  routes  are  normally  driven 
out  of  the  solution  when  resources  are  constrained.  The  shorter  routes 
are  chosen  in  order  to  effectively  use  the  resources. 

c.  Tolerance.  In  the  RACAT  algorithm  there  are  calculations  which 
would  result  in  a  zero  value  if  numerical  accuracy  were  absolute.  Since 
there  are  problems  in  numerical  significance— -and  there  would  be  no 
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matter  how  large  the  computer  word— these  calculations  can  reeult  in 
sane  very  small  number  rather  than  zero.  For  example,  in  the  compu¬ 
tation  of  simplex  multipliers  the  program  treats  any  simplex  multiplier 
less  than  .00000001  as  if  it  were  zero.  This  value  is  known  as  the 
"tolerance"  value.  It  eliminates  any  iteration  that  might  result  in 
only  a  tiny  improvement  or  which  might  be  caused  by  small  errors  due 
to  numerical  problems. 

4.  Output  Processes.  Each  module  of  the  ETNAM  computer  system  has 
some  output  processes.  The  output  consists  of  files  and  printed  output. 
All  files  are  written  in  binary  rather  than  BCD]  this  precludes  the 
need  for  converting  the  data  to  BCD  and  introducing  round-off  errors. 

In  the  printed  output,  ,ne  FORTRAN  F  format  has  been  used  whenever 
possible  to  print  decimal  numbers  (e.g.,  126. 732)  for  ease  in  reading. 

If  range  of  a  number  is  large,  to  avoid  overflow  or  lack  of  significance, 
the  FORTRAN  E  format — the  scientific  notation — has  been  used  (e.g., 
the  number  .000032  would  be  printed  as  3*2  E-05,  the  E-05  standing  for 
10'6). 
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ABSTRACT 


The  problem  of  determining  multicommodity  flows  over  a 
capacitated  network  subject  to  resource  constraints  may  be  solved 
by  linear  programming.  However,  the  number  of  arce  and  re¬ 
sources  in  many  applications  is  each  that  the  standard  formula¬ 
tion  becomes  very  costly.  This  paper  describes  an  approach— an 
extension  of  the  column  generation  technique  used  in  the  multi- 
commodity  network  flow  problem— that  simultaneously  considers 
r  twork  chain  selection  and  resource  allocation,  thus  making  the 
problem  both  manageable  and  optimal  in  the  sense  that  flow  at¬ 
tainable  is  constrained  by  resource  availability  and  network  ca¬ 
pacity.  Extension  to  the  minimum-cost  formulation  is  proposed 
end  computational  experience  is  discussed. 
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1.  INTRODUCTION 

The  problem  of  multicommodity  flows  in  capacitated  networks  has  re¬ 
ceived  considerable  attention.  Ford  and  Fulkerson1  suggested  a  computational 
procedure  to  solve  the  general  maximum-flow  case.  Tomlin2  has  extended  the 
procedure  to  include  the  minimum-cost  case.  Jewell2  has  pointed  out  the  strong 
historical  and  logical  connection  between  this  solution  procedure  for  the  multi- 
commodity  problem  and  the  decomposition  algorithm  of  Dantzig  and  Wolfe.4 

A  related  problem,  which  has  not  been  directly  addressed,  is  the  deter¬ 
mination  of  multicommodity  flows  in  a  system  constrained  by  resource  avail¬ 
ability.  For  example,  flows  in  transportation  networks  are  constrained  by 
available  resources  that  must  be  shared  by  two  or  more  arcs  in  the  network. 
The  determination  of  the  set  of  routes  and  the  allocation  of  resources  to  these 
routes  to  maximize  multicommodity  flows  or  to  minimize  system  cost  in  meet¬ 
ing  fixed  flow  requirements  can  be  applied  to  many  problems  in  logistics  and 
other  areas.  This  paper  discusses  solution  procedures  for  multicommodity 
network  flows  with  resource  constraints  in  the  maximum-flow  and  minimum- 
cost  cases. 


2.  THE  MAXIMUM  MULTICOMMODITY  NETWORK  FLOW  PROBLEM 

Consider  the  multimode,  multicommodity  network  G(  N,fl).  N  is  the  set 
of  all  the  nodes  of  the  network.  S  is  the  subset  of  all  ordered  pairs  (x,y)  of 
the  elements,  of  N  that  are  arcs  of  the  network.  2lf .  .  .  ,2m  is  an  enumeration 
of  the  arcs.  Each  arc  has  an  associated  capacity  b(x,y)  2  0  and  an  associated 

cost  (or  distance)  y)  *  0. 

* 

For  each  commodity  k(  k  =  1, .  .  .  ,q)  there  is  a  source  sk  and  a  sink  tfc. 
The  flow  of  commodity  k  along  a  directed  arc  ( r,y)  is  y**  y> ,  (k  =  1, .  .  .  ,q) 
and  these  y(x  y) ,  (k  =  1, .  .  .  ,q)  must  satisfy  the  capacity  constraints 

i  vk  <  b  [( x .  y )  ,  a] 

k  =  l  U.y)  <».»*-  J 
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The  Kirchhoff  conservation  requirements  must  be  satisfied  for  each  com 
modity  fc(fc  =  1 . q)  and  for  all  x  f  N  as 


I 

n  <  N 


it  *  -  Sfc 

0  if  x  4  sk,  ik 

if  *  *•  'k 


where  vk  is  the  total  flow  of  commodity  k  from  sk  to  tk. 

The  multicommodity,  maximum-flow  problem  as  formulated  by  Ford  and 
Fulkerson1  is  as  follows:  Define  the  set  Pk  =  jp^l  Pj^  is  a  chain  connecting  sk 
and  ikj.  Nov  let  P  be  the  union  of  the  sets  Pk(k  =  1, .  .  .  ,q).  Further,  let  Pj(l), 
Pj(U, .  .  .  ,  Pj(W, .  .  .  ,  P^  be  the  enumeration  of  the  chains  p|W  c  P  such  that 
the  subscript  j  is  sufficient  to  id  itify  the  chain,  its  origin-destination  pair, 
and  the  commodity  with  which  it  is  associated. 

Thus  the  kth  set  is  defined  as 

]  k  =  1 1  |  P.lk)  is  a  chain  from  sk  to  ,  k  -  1 

The  arc-chain  incidence  matrix  is 

A  -  [«,,] 


where 


1  if  a,  f  p(k) 

0  otherwise 


for  i  =  1 . m;  j  *  1,  .  .  . 

Consider  the  network  used  as  an  example  in  Ref  1,  augmented  by  sfc  and 
tfc(k  =  1,2),  with  source  Si  and  sink  f!  for  commodity  1  and  source  s*  and 
sink  1 2  for  commodity  2.  Figure  1  illustrates  the  network  and  Fig.  2  shows 
the  arc-chain  incidence  matrix  A. 

Letting  ijk)  (  /  =  1, .  . .  ,n)  be  the  flow  of  commodity  k  in  chain  Pj-k>  (j  =  1, 
.  .  .  ,n;  k  implicit)  and  the  flow  capacity  of  Sit  the  multicommodity,  maxi- 


mum-flow  linear  program  is: 

Maximize 

rt 

l  *<k) 
j-i  i 

subject  to  capacity  constraints 

n 

2,  a  x 

i  - 1  •» 

4 


2 


Fig-  1— Network  A 
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Fig.  2— Arc-Chain  Incidence  Matrix  A 
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for  i  =  1, .  .  .  ,  m. 

Thus  the  objective  is  to  maximize  flow  over  all  possible  chains  from 
origins  to  their  respective  destinations  subject  to  the  capacity  constraints  of 
the  arcs. 

The  number  of  variables  in  the  aforementioned  linear  program  is  very 
large  since  the  number  of  possible  chains  is  very  large  in  most  applications. 
The  procedure  proposed  by  Ford  and  Fulkerson1  is  to  treat  the  nonbasic  vari¬ 
ables  implicitly;  i.e.,  nonbasic  chains  are  not  enumerated.  The  column  vector 
to  enter  the  basis  is  generated  by  applying  the  simplex  multipliers  to  the  arcs 
as  pseudo  costs  and  selecting  the  candidate  chain  using  the  shortest  chain 
algorithm.* 

The  remainder  of  this  paper  is  concerned  with  an  extension  of  this  pro¬ 
cedure  to  allow  the  addition  of  resource  constraints  on  the  flow  over  the  net¬ 
work.  The  maximum-flow  problem  will  be  discussed  first,  followed  by  a  mini¬ 
mum-cost  formulation.  A  brief  description  of  experience  with  a  computer 
program  using  this  procedure  concludes  the  paper. 

Extension  To  Include  Resource  Constraints 


In  the  linear  programming  problem  stated  previously,  flow  is  to  be  rnax- 
imii’.ed  subject  to  the  constraints  imposed  by  the  capacities  of  the  Individual 
arcs  of  the  network.  In  some  applications  additional  constraints  on  flow  are 
imposed  by  the  limited  availability  of  resources  used  jointly  by  two  or  mr'1 
arcs  of  the  network.  An  example  of  this  type  of  network,  which  will  be  useu 
throughout  the  remainder  of  this  paper,  is  a  transportation  network. 

It  is  clear  that  the  simultaneous  consideration  of  both  types  of  constraints 
is  an  important  problem  in  transportation  networks.  Roadways,  rail  lines,  etc, 
have  capacity  limitations  that  may  limit  the  maximum  movement  of  men  and 
materials,  particularly  in  less-developed  areas.  The  vehicles  and  resources 


*Professor  Mandell  Bellmore  of  The  JohnB  Hopkins  University  and  Mr.  Donald 
Boyer  of  the  Logistics  Research  Project,  The  George  Washington  University,  have  de¬ 
veloped  computer  programs  to  solve  the  problem  using  thiB  procedure. 


. ! 
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available  to  use  the  network  can  actually  impose  a  greater  constraint  on  total 
movement  than  the  arc  capacities.  In  a  highly  developed  transportation  system 
the  capacity  of  the  network  may  greatly  exceed  that  required;  the  effective 
limitations  of  movement  result  from  too  few  vehicles  or  other  resources. 

To  effect  this  simultaneous  consideration  of  resources  as  related  to  the 
network,  the  equipment  required  to  sustain  a  flow  of  1  unit  across  an  arc  (x,y) 
could  be  considered.  For  example,  consider  some  of  the  equipment  that  might 
be  necessary  over  the  network: 

(1)  General-purpose  trucks 

(2)  Tank  trucks 

(3)  Trains  (locomotives  and  cars) 

(4)  Barges 

(5)  Mechanics 

(6)  Drivers,  laborers,  etc 

•  •  •  • 

(p)  Other  resources 

Such  a  list  may  be  as  detailed  or  as  aggregated  as  desired.  For  example, 
standard  notional  trains  consisting  of  a  locomotive  and  its  complement  of  rail 
cars  may  be  used,  or  the  resources  may  be  considered  separately.  Some  no¬ 
tional  truck  combining  all  the  various  trucks  available  may  be  considered,  or 
general-cargo  trucks  and  tank  trucks  may  be  listed  separately. 

Let  the  resource  matrix  for  commodity  k  be 

Rk  «■  [r^]  ( i  •  1,  s  -  1,  .  .  •  ,  p) 

where  is  the  quantity  of  resource  s  required  to  sustain  a  unit  flow  of  com¬ 
modity  k  over  arc  i;  r*s  2  0.  Note  that  for  some  arc  commodity  combinations 

rk  -  w  (s  =  1,  .  . .  ,  p) 

1$ 

e.g.,  if  the  arc  represents  a  pipeline  and  the  commodity  is  passengers. 

Letting  pg  be  the  quantity  of  resource  $  available  (e.g.,  in  inventory)  for 
assignment  to  the  network  (s  *  1, .  .  .  ,p),  the  linear  programming  problem 
with  resource  constraints  is: 
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maximize 


i  i'*' 
i-i  i 


subject  to  (a)  capacity  constraints 


for  i  =  1, .  . .  and  (b)  resource  constraints 

ill  fl„*<k>  rk  <  a. 

i-i  k » i  j<;k  ''  i  “ 

for  s  ■  1, .  . .  ,p. 

In  matrix  form  the  flow- maximizing  linear  program  with  arc  capacity 
and  resource  constraints  is: 


maximize 


I  xtk> 
f  =  i  > 


subject  to  AX  s  b  where  A  is  a  matrix  ( »  +  p  x  n)  formed  of  two  submatrices 


A  and  E. 


A-  ± 
*  E 


A  ■  [a,,]  is  the  arc-chain  incidence  matrix.  E  =>  [e  ( ]  is  the  resource- 


chain  requirement  matrix  where 


esj  -  2  rk 

sl  IS  'i 


for  k  where  P*k)  connects  sk  and  tk,  s*l . p. 


X  -  col.(x<k>,x(k> . x,k)  \ 

\  1  2  m*fJ 


b  -  col.  (  bj ,  b2 . p2 pp) 


In  words,  the  A  matrix  is  composed  of  column  vectors  that  represent 
chains  from  sk  to  tk  in  a  natural  way.  The  chain  is  identified  by  the  arcs  that 
form  it  and  by  the  resources  required  to  sustain  a  unit  flow  from  sk  to  tk .  A 
typical  column  vector  of  A  is 

A,  *  c°l'  («ira2, . %•  . 
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Solution  Using  the  Chain  Generation  Technique 

The  flow-maximizing  linear  program  with  arc  capacity  and  resource 
constraints  will  be  quite  large  for  most  applications  and  will,  in  addition,  re¬ 
quire  considerable  preliminary  computation  to  obtain  the  coefficients  of  the  A 
matrix.  (The  authors  have  solved  several  small  problems  using  a  standard 
linear  programming  code.)  It  will  now  be  show:,  uii:  the  column  generation 
procedure  suggested  by  Ford  and  Fulkerson1  can  be  modified  to  apply  to  the 
problem  extended  to  include  resource  constraints  so  that  it  is  never  necessary 
to  form  the  A  matrix  explicitly. 

The  linear  program  described  previously  in  matrix  form  is  now  modified 
to  a  standard  linear  programming  form.  (An  attempt  has  been  made  throughout 
this  discussion  to  follow  the  notation  used  by  Hadley.8)  The  inequalities  are 
changed  to  equalities  and  *+  p  slack  variables  are  added.  A  cost  or  profit 
vector  is  included  in  the  calculation  of  the  objective  function.  The  problem 
then  becomes: 
maximize 

CX  -  z 

subject  to 

ax  -  b 

where  A  =  [A;Im+p] 

C  =  (Cj, . . . ,  cn+m+p) 

Cj  =  the  profit  obtained  for  each  unit  of  flow  achieved  on  P-k)  * 

*  =  Ml-  («i“ . 


An  initial  basic  solution  is 

B0  *b  -  * 

with  objective  function 

CbXb  -  *  =  0 


*In  the  previous  formulation  of  this  problem  the  value  of  the  c,’b  ( j  =*  1,  .  .  .  ,  it) 
was  assumed  to  be  1.  That  is,  the  value  of  a  unit  of  flow  achieved  on  any  chain  is  equal 
to  that  achieved  on  any  other  chain.  This  is  the  usual  case.  If  any  Cjji  1  then  one  must 
require  that  cr  =  c,for  all  r,  s  c  J* ,  for  k  =1,  .  .  .  ,q,  in  order  for  the  proposed  algorithm 
to  be  valid.  For  c,  /  c,  (r,  sc  Jk)  implies  that  of  the  chaine  rom  sk  to  ik  some  are 
preferred  over  others. 
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where  B0  =  lm+p 

^b  =  ^cji«  cjj . Cj»t+P* 

y  3  coi  /-(W  T(k)  JM  » 

Ab  co*‘  '*U  >  . *jiti+p' 

such  that  j, ,  j, , .  .  .  ,  jm  +p  are  selected  so  that  Cb  and  Xb  conform  properly  to 
the  Bo  matrix,  which  is  the  basis. 

Thus  one  begins  with  a  feasible  solution,  and  the  simplex  procedure  is 
fundamentally  a  method  for  improving  this  solution  by  selection  vectors  from 
the  A  matrix  to  enter  the  basis.  Given  a  feasible  solution  to  the  linear  pro¬ 
gramming  problem  as  expressed  previously,  the  value  of  the  new  objective 
function  (following  the  introduction  of  a  new  vector)  is:* 

x  «  z  +  x B r / y, j  (Cj  -  Xj)  «  X  +  e  [cj  -x() 

Thus,  for  a  suitably  selected  8  the  value  of  the  new  objective  function  will  be 
greater  than  the  old  if,  and  only  if, 

CJ  ~  x,  >  0  (1) 

Now,  since  B  is  a  basis,  every  A^  (a  candidate  to  enter  the  basis)  can  be  ex¬ 
pressed  as 

■♦p 

A(  *  y(Bi  +  vjbi  +  ■  •  • +  -  (Vi8! 

the  Bj ’s  being  the  *+  p  column  vectors  of  B  and  the  y^s  the  unique  scalars 
such  that  A(  is  a  linear  combination  of  the  Bj .  Now  let 


so  that 


A(  -  BY, 


Y,  .  B-‘A, 
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Now  the  present  objective  function  (i.e.,  before  the  introduction  of  a  new  Aj ) 
is  in  terms  of  the  present  basis: 

"ip  <»> 

1  “  CB*B  “  j5j£BtxB( 

But  since  Aj  is  expressible  in  terms  of  this  basis,  there  is  for  each  Aj ,  a  can¬ 
didate  to  enter  the  basis,  a  new  £j  as  follows: 

“  Vlj  CBI  *  >21  CB2  *  •  •  •  >M  +  p,f  ‘  B 

PP  “  P 


or 


*1  -  £avi 


(3) 


Thus,  substituting  Eq  3  in  Eq  1, 

C/  -cay,  >0  (4) 

and  substituting  Eq  2  in  Eq  4, 

c(  -  Cg  B“l  Aj  >  0 

Now,  bringing  Cj  to  the  right-hand  side  and  multiplying  through  bv  - 1 , 

c  B  B  1  Aj  <  £j  (5) 

which  is  the  simplex  rule  and  a  criterion  for  choosing  a  new  Aj  to  enter  the 
basis.  If  CgB-'Aj  <  ^  then  Aj  can  be  brought  into  the  basis  and  the  new  value 
of  the  objective  function  will  be  equal  to  or  greater  than  the  preceding  value. 

Applying  the  simplex  criterion  in  the  stated  maximization  problem,  one 
may  choose 

min  (z,  -  ct) 

I  1  1 

or 

i»in(cB  B-1  Aj  -  Cj) 

and  when  all 

cBB-‘Aj  -  Cj  >  0 


an  optimum  solution  has  been  found. 

It  is  the  hypothesis  of  the  Ford  and  Fulkerson  suggested  solution  to  the 
multicommodity  problem  (and  of  the  present  resource-constrained  algorithm) 
that  the  shortest  chain  algorithm  can  be  used  to  develop  the  Aj  that  will  satisfy 
the  simplex  rule.  Further,  if  the  shortest  chain  algorithm  can  find  no  chain 
satisfying  the  requirement,  an  optimum  has  been  reached. 
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Let  cbB-‘  =  (ou,  aa , .  .  .  ,  am,  Vt,  *>••••  >  ffp).  The  shortest  chain  in 
the  multicommodity  problem  without  resource  constraints  may  be  found  by 
assigning  the  to  the  corresponding  arc  i  so  that  d{][  y)  =  a(  ,*  where  3(  =  (x,y); 
the  shortest  chain  algorithm  can  then  be  used  to  find  the  chain  connecting  sk 
to  tk  (k  »  1, .  .  .  ,4)  such  that  cB B_1Aj  —  Cj  is  minimized.  Further,  if  the  shortest 
chain  c8 B“l  Aj  2  c(  then  no  chain  exists  for  which  cbB-1Aj  <  Cj  and  the  optimum 
has  been  reached. 

Now  in  the  resource-constrained  problem  this  procedure  must  be  modified 
since  all  the  elements  of  cBB~k  cannot  now  be  assigned  directly  to  arcs.  This 
is  because  ir,  to  wp  correspond  to  resources,  not  arcs.  Therefore  a  "pseudo 
cost"  must  be  devised  that  can  be  assigned  directly  to  arcs  such  that  the  short¬ 
est  chain  will  be  the  chain  with  minimum  (cBB-lAj  -  ).  Now 


cb®  (  a 1 •  aa  >  •  • 


• ,  v 


where 


"fl»‘ 

- 

*aM 

a-l 

a  _ 

a 

"J 

mi 

m 

ij 

1  rk  fl,, 

C7) 

i-i  •*  " 

■  . 

t  m 

2  r‘  «„ 

Pi 

l-l  T  11 

L-  J 

•If,  however,  any  a,  <  0  or  W,  <  l  the  shortest  chain  algorithm  may  (ail.  Ths  slack 
vecicr  corresponding  to  the  negative  or,  or  should  therefore  be  brought  in.  This  is  a 
unit  vector  (  u,)  with  a  1  in  the  1  th  position  (a  1  in  the  *  *  ith  position  in  the  case  of  f,  <0) 
and  zeros  elsewhere.  For  if  a,  <  0,  then  cuB^U.-c  c, ,  since  t,  *  0  (|  » 1 . *  + 
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for  k  where  pjk)  connects  $k  k .  Thus 

CB  B-1  A|  -  ^  ^oc,  atj  +  Ji( 

and  tho  contribution  of  each  arc  i ,  member  of  the  chain  f](k> ,  is 

01,  atj  +  it,  (for  k,  where  P^1*1  connect*  $k  end  1„). 


Define  a  real  numbered  function  d*  as 

dk  «  a,  a,,  +  £  ir,  r* 

i  '  */  ,.i  *  i<  i) 


and  assign  d*  to  each  arc  i.  The  shortest  chain  algorithm*  may  then  be  used 
to  find 


min  I 

I 


[a  *  d*"l  *  *«  „  S  |  a,  e„  +  £  ir,  rk  «  I 

M  'J  i  \  M  *  “  *7 


-  min  c8  B-,A^  (for  k  -1 . 4) 

Thus  the  candidate  Ay 'is  found  to  enter  the  basis.  The  column  vector  to  leave 
the  basis  may  be  determined  in  the  normal  simplex  manner.  The  procedure 
for  determining  the  maximum  multicommodity  network  flows  with  resource 
constraints  is  as  follows: 

(1)  Form  an  initial  basis  B„  -  I„+p. 

(2)  Calculate  B"1, cbB-1  «(a, . . tp),andi»  cBB-‘k. 

(3)  Evaluate  (at . a*,  . . itp  )  as  follows: 

(a)  If  any  a,  <  0,  enter  into  the  basis  the  slack  vector  for  that  arc, 
which  is  a  unit  vector  with  a  1  in  the  ith  position,  zeros  elsewhere,  and  go  to 
step  2. 

(b)  If  any  ir,  <  0,  enter  into  the  basis  the  slack  vector  for  that  re¬ 
source,  which  is  a  unit  vector  with  a  1  in  the  (is  +  t  )th  position,  zeros  else¬ 
where,  and  go  to  step  2. 

(c)  If  no  a, ,  t ,  <  0,  go  to  step  4. 

(4)  Determine  the  min  [c8B_1Af  -  c,  ]  as  follows: 

(a)  For  k  -  1, .  .  .  ,  4,  calculate  and  assign  d*  to  arcs  <,...,■  and 
find  the  shortest  chain  from  sk  to  (k . 
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(b)  From  the  q  shortest  chains,  select  that  chain  with  the  minimum 
[cBB-lAj  -  Cj  ].* 

(5)  If  the  minimum  chain  results  in  cbB-1Aj  a  Cj  ,  terminate  the  algorithm 
since  there  exists  no  chain  that  will  increase  the  value  of  z ;  otherwise  go  to 
step  6. 

(6)  Enter  the  minimum  Aj  of  step  5  into  the  basis.  Determine  the  vector 
to  leave  the  basis  in  the  normal  simplex  manner.  Go  to  step  2. 


3.  THE  COST -MINIMIZATION  PROBLEM 


In  many  applications  the  objective  is  to  minimize  the  cost  of  resource 
utilization  to  meet  fixed  delivery  requirements  at  given  destinations.  Tomlin 
has  proposed  such  an  extension  to  the  arc -chain  formulation  cf  the  maximum - 
flow  linear  program  2  In  this  section  of  the  papei  omlin’s  formulation  is 
modified  to  include  resource  constraints  and  a  revised  column  generation 
scheme  as  a  solution  procedure  is  proposed. 

The  minimum-cost,  multicommodity  flow  problem  with  resource  con¬ 
straints  may  be  formulated  as  follows :t 
minimize 


subject  to  (a)  capacity  constraints 


X  o  .  i<l<)  <  b,  for  i  -  1 . m 

j=  ')  ) 


(b;  i  t  source  constraints 


XXX  a..  x,,">  r*  <  p„  for  s  =1,  . 
.  =  1  =  l  jdk  ‘1  J  is  5 


and  (c)  delivery  requirements 


X  X*'  =  Ak  for  k  =  1 .....  q 


where  Afe  is  the  delivery  requirement  at  tfe  (  fe  =  1 . q),  (\k  *  0). 


*The  algorithm  may  be  modified  to  accept  the  first  chain  where  cbB-1  Aj  <.  Cj. 
tThe  notation  used  in  this  section  parallels,  where  applicable,  that  used  in  the  pre¬ 
ceding  maximum-flow  section. 
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vamps 


The  cost  coefficient  c ^  may  be  defined  as: 

*  £ 

Cj  -  ^  r,  at)  +  4sesi  for  j  -1 .....  n 

-  £  r, +  £  £  i.rk  0,,  (for  j  «1,  .  .  .  ,  »;  k  where  Plk)  connects  sk  and  tk) 

j.l  1  'I  s-l  i.i  *  i*  'J  )  *  * 


where  r<  is  the  cost  (or  toll)  for  a  unit  flow  over  arc  i ,  and  0S  is  the  cost  of 
using  a  unit  of  resource  s. 

Define  the  matrices  G  and  A  as  follows:  G  is  a  commodity  delivery  in¬ 
cidence  matrix  ( q  x  it) 

G  - 


where 


8kj 


1  if  jrjk* 
0  otherwise 


A  is  a  matrix  ( m  +  p  +  q  x  n)  formed  of  the  submatrices  A ,  E,  and  G  as  'ollows: 


A  = 


Now  the  typical  column  of  A  is 

A.  m  col.  («jj,  •  •  •  .  <lj . «pj.  8lj . 


In  the  matrix  notation  the  cost-minimization  problem  is  then  to 
minimize 

CX 

subject  to 

AX  =  b 

where 

A  = 
and 

C  =  (Cj.  •  •  •  .  cn  +  m  ■*  p  +  <1  ^ 

The  Cj  are  the  costs  associated  with  the  chains  O'  =1,  .  .  .  ,n),  the  slack 
variables  (f=rs  +  l,...,n+m+p),  and  the  artificial  variables  (/  =  n  +  m  +  p+l, 

.  .  .  ,»  +  iiitp+(|). 

'♦Note  that  all  colimn s  of  g  are  unit  vectors.  No  chain  meeting  the  Kirchhoff  con¬ 
servation  requirements  may  be  a  member  of  more  than  one  of  the  sets  Jk ,  h 
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Adopting  the  standard  two-phased  solution  procedure  for  the  minimum- 
cost  formulation,  the  Phase  I  procedure  minimizes  to  zero  the  value  of 

|  =  n<  n*p*  1  * 

to  obtain  n  initial  basic  feasible  solution.  This  effectively  assigns  a  cost  of 
1  to  the  artificial  variables  and  a  cost  of  zero  to  the  other  variables  in  Phase  I. 
Phase  II  begins  with  the  basic  feasible  solution  determined  in  Phase  I  and  pro¬ 
ceeds  to  minimize 


Phase  1 

In  Phase  I  Iw+p^maybe  used  as  the  initial  basis  and  the  simplex  rule  is 
to  enter  a  chain  in  the  basis  if,  and  only  if, 

c,  -‘bb"Ia;  <  0 

where 

c B  0-1  =  (ai - a*.  . V  . V 


so  that  the  simplex  multipliers  are  associated  with  the  arcs,  the  irs  are 
associated  with  the  resources,  and  the  ok  are  associated  with  the  artificial 
variables.  Thus  the  vector  Af  is  entered  if 


c(  -  cB  B  1  A( 


m 

Z  a 

i  =  i  1  11 


9ltj  <  0 


P 

1 

S  =  I 


”s  “i,  r~s  - 


<  0 


Since  each  chain  is  associated  with  one  and  only  one  of  the  tfe  (k  =  1, . . .  ,q), 
and 

then  A.  is  to  enter  the  basis  if 
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Thus  the  contribution  of  each  arc  to  Cj  -  cBB~lAj  in  Phase  I  is 


ik  =  - 


p 

2 


S  =1 


” s 


One  may  then  use  the  shortest  chain  algorithm  as  in  the  maximum-flow  case 
to  find,  for  each  it , 


Thus,  one  may  find  the  candidate  A(  to  enter  the  basis.  The  column  vector  to 
leave  the  basis  may  be  determined  in  the  usual  manner. 

Should  any  a{  or  rr,  be  positive  one  can,  in  a  manner  analogous  to  the 
maximum-flow  case,  bring  in  the  corresponding  slack  variable.  Should  the 
minimum  cbB~*  Aj  be  greater  than  or  equal  to  zero,  there  is  no  feasible  solu¬ 
tion  and  the  procedure  would  be  terminated.  If  there  is  a  feasible  solution, 
Phase  I  will  terminate  when  all  artificial  variables  have  vanished  and  2=0. 


Phase  II 

Phase  II  begins  with  a  basic  feasible  solution  with  all  artificial  variables 
equal  to  zero.  Its  simplex  criterion  states  that  a  chain,  represented  by  the 
vector  A. ,  may  enter  the  basis  if,  and  only  if, 


‘,-cB  B-'A,  <0 


In  Phase  n, 


:i  ~CB°  lAj  " 


Skj 


P  m 


X  ajf  ( r.  -  a ,)  +  2  2  au  r* ( <t> .  -  n J  -  2  ok  gki 

j  =  l  1  '  S  =  1  i  =  l  •)  IS  S  S  It  - 1  * 


m  f  P  "j  1 

2  flg.  I  (r(  -Kj)  +  2  rk  -  n  )  -  2  ( 

i  =  l  '•  L  *  =  1  «  5  s  J  k  =  1 


k  8kj 
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**?*:!: 


Thus 


<*■'  -  ri  -<** + 


is  the  analogue  of  d*  in  the  maximum-flow  case  and  may  be  assigned  as  a 
pseudo  cost  to  each  arc.  The  shortest  chain,  i.e.,  the  chain  with  the  least, 

1  ik'  -  a,  <  0 

i»l  >  * 

for  it  =  1, .  .  .  ,  q,  may  then  be  entered  in  the  basis. 

The  Phase  II  simplex  rule  is,  then,  enter  the  vector  Af.  if 


or 

As  in  Phase  I,  a  positive  at  or  trs  will  indicate  that  the  corresponding 
slack  variable  should  enter  the  basis.  Phase  II  is  terminated  and  the  value  of 
z  is  minimized  when,  for  the  minimum  j , 

* 

UjfPj  1 


4.  COMPUTATIONAL  EXPERIENCE 

A  computer  program  for  the  IBM  7044  has  been  developed  to  solve  both 
the  maximum-flow  and  minimum-cost  problems.  This  test  pi  ogram  is  an 
all-in-core  routine  that  updates  the  basis  inverse  each  iteratior  A  series  of 
test  cases  has  been  solved  successfully.  The  largest  of  the  series  had  108  arcs, 
9  types  of  resources,  and  3  commodities.  Several  variations  of  this  case  have 
been  tested  including  both  network-constrained  and  resources-constrained  prob¬ 
lems.  Two  typical  problems  were  also  solved  by  "numerating  the  possible 
chains  for  the  LP  3  linear  programming  code.  A  comparison  of  the  results  of 
the  two  methods  indicates  that  the  proposed  algorithm  is  accurate  and  signifi¬ 
cantly  faster  than  the  standard  code. 
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ABSTRACT 


In  a  prior  paper  (RAC -P -44,  "A  Model  for  Optimal 
Multicommodity  Network  Flows  with  Resource  Alloca¬ 
tion,”1)  a  methodology  for  determining  optimal  multicom¬ 
modity  flows  overs  capacitated  network  with  simultaneous 
consideration  of  network  and  resource  constraints  was 
developed.  This  paper  describes  an  extension  to  this  basic 
work  to  Incorporate  the  ability  to  handle  resource  sub¬ 
stitutions.  Substitution  of  resources  is  the  capability  to 
employ  alternative  resource  combinations  to  accomplish 
flown  in  the  network.  The  extended  model  permits  the 
selection  of  an  optimal  set  of  resources  and  chains  to  ac¬ 
complish  required  flows  in  a  network. 
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INTRODUCTION 

In  an  earlier  paper  (RAC-P-441)  a  model  was  presented  for  the  solution 
of  multicommodity  network  flow  optimization  problems  with  resource  con¬ 
straints.  Included  in  the  paper  was  a  description  of  computational  experience 
to  d.ce.  The  model  formulation  was  based  on  an  extension  to  the  shortest- 
chain  procedure  for  solving  multicommodity  flow  problems  suggested  originally 
by  Ford  and  Fulkerson1  in  1958  and  further  developed  by  Tomlin,9  Bellmore,4 
Boyer,1  and  others. 

The  extended  formulation  considers  the  problem  in  which  resources  must 
be  assigned  to  the  arcs  of  a  chain  in  order  to  sustain  a  flow  over  that  chain. 

For  example,  consider  a  transportation- network  problem  in  which  trucks, 
drivers,  gasoline,  trains,  barges,  etc,  must  be  allocated  to  the  various  routes 
if  any  movement  is  to  be  accomplished.  The  extended  formulation  thus  Includes 
(in  addition  to  arc-capacity  constraints)  resource-utilization  constraints  that 
prevent  the  assignment  of  more  of  any  type  of  resource  than  is  available.  There¬ 
fore,  in  the  maximum -flow  mode,  routes  are  selected  and  resources  allocated 
to  attain  maximum  network  flow.  In  the  minimum-cost  mode,  routes  ara  se¬ 
lected  and  resources  allocated  so  as  to  minimize  the  cost  of  meeting  flow  re¬ 
quirements  subject  to  the  arc  capacity  and  resource-utilization  constraints. 

This  paper  presents  a  modification  of  the  initial  extended  formulation  to 
allow  for  the  substitution  of  resources.  In  the  previous  paper  only  one  com¬ 
bination  of  resources  was  permitted  to  be  applied  to  an  arc  in  orde.  to  move 
one  uni:  of  commodity  fc  over  arc  i.  In  economic  terms  the  arc-commodity 
pair  is  similar  tc  a  production  function  with  constant  returns  to  scale  and  fixed 
technical  coefficients  (see  R.  G.  D.  Allen,  Macro-Economic  Theory,*  p  36). 

In  some  applications  this  may  be  a  significant  limitation.  Consider  again  a 
transportation  network.  A  highway  arc  might  be  considered  for  the  transport 
of  manufactured  products.  Closed  vans  wi‘h  a  driver  and  an  alternate  driver 


3 


might  be  the  most  efficient  combination  of  resources.  A  combination  of  a  van 
and  one  driver  would  be  less  efficient,  perhaps,  since  more  rest  periods  would 
be  required,  but  it  is  nevertheless  a  feasible  combination.  Similarly  a  third 
alternative  would  be  the  utilization  of  stake  and  platform  trucks  with  containers, 
possibly  more  expensive  than  the  first  two  alternatives  but  still  feasible. 

Specific  inventories  of  trucks  and  drivers  ma,  exist  and  the  objective 
might  be  to  assign  these  sets  of  resources  in  the  most  efficient  way  over  all 
arcs  even  if  some  arcs  or  commodities  are  assigned  a  less  than  most  efficient 
set  of  resources. 

An  extension  to  the  multicommodity  network  flow  model  with  resource 
constraints  is  proposed  that  permits  alternative  resource  combinations  to  be 
employed  to  accomplish  flows  over  arcs  ani  chains. 

NOTATION 

To  familiarize  the  reader  with  the  notation  used  in  the  development  of  the 
extended  multicommodity  network  flow  model,  this  section  of  the  report  is  de¬ 
voted  to  a  brief  review  of  the  mathematical  structure  of  the  resource-constrained 
multicommodity  network  flow  model.  This  section,  as  wei.'  as  the  following  sec¬ 
tion  on  substitu'.  in  of  resources,  emphasizes  the  minimum-cost  mode  formula¬ 
tion  rather  than  the  maximum-flow  mode  (minimum  cost  being  the  more  detailed 
of  the  two  similar  formulations). 

Consider  the  multimode  multicommodity  network  G  (N,  3).  N  is  the  set  of 
all  the  nodes  of  the  network.  3  is  the  subset  of  all  ordered  pah«  (x  ,y)  of  the 
elements  of  N  that  are  arcs  of  the  network.  3k , .  .  .  ,flm  is  an  enumeration  of 
the  ares.  Each  arc  has  an  associated  capacity  b<i  y)  2  0  and  an  associated  cost 
(or  distance)  d(xy)2  0. 

For  each  commodity  k  (k  -  1 . q)  there  is  a  source  sk  and  a  sink  tk. 

The  flow  of  commodity  k  along  a  directed  arc  (x,y )  is  y),(k  ■  1 . q),and 

these  y^  ,  (k  =  1, ...  ,q)  must  satisfy  the  capacity  constraints 

£(i! 

Define  the  set  Pk  ■  ( P(lk>  |  ?,<k)  is  a  chain  connecting  sk  and  tk }.  Now  let  P 
be  the  union  of  the  sets  Pk  { k  -  1,  .  .  ,q).  Further,  let  Pj1’ ,  P,1’ ,  .  .  .  ,  °fk) , 
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.  .  .  ,  pW  be  the  enumeration  of  the  chains  Pfk)  €  P  such  that  the  subscript  j 
is  sufficient  to  identify  the  chain,  its  origin -destination  pair,  and  the  commodity 
with  which  it  is  associated. 

Thus  the  kth  set  is  defined  as 

] k  -  I  >•  8  chain  from  sk  to  t ,  k  - 1,  .  .  .  ,4 

The  arc -chain  incidence  matrix  is 

A  -  [.„] 

where 

( 1  if  C,  £  pJ<fc) 

a„  -{ 

\  0  otherwise 

for  1  *  1  f  •  •  •  1  wtj  j  2  1^  •  •  «  jNi 

Let  the  resource  matrix  for  commodity  k  be 

Rk  -  [rk,  ]  ( i  -  1 . .  s  »  1,  .  .  ,p) 

where  r{“s  is  the  quantity  of  resource  s  required  to  sustain  a  unit  flow  of  com¬ 
modity  k  over  arc  i ;  *  0.  Note  tlu'.t  for  some  arc-commodity  combinations 

r,  -  -  (s  -  1 . p) 

e.g.,  if  the  arc  represents  a  pipeline  and  the  commodity  is  passengers. 

Let  x(k)  (/  =»  l . n)  be  the  flow  of  commodity  k  in  chain  Pfk>  (j  *  1, 

.  .  .  , n ;  k  implicit),  b,  the  flow  capacity  of  3j,andpsthe  quantity  of  resource 

s  available  (e.g.,  In  inventory)  for  assignment  to  the  network  (s  =  l . p). 

The  minimum-cost  multicommodity  flow  problem  with  resource  constraints 

may  be  formula*  xi  as  follows: 

minimize 


subject  to  (a)  capacity  constraints 


£  a,,  x!kl  £  6,  for  i  -  1,  ....  m 
j-l  'I  >  ‘ 


(b)  resource  constraints 
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Now  the  typical  column  of  A  is 

A.  -  col.  (atj . amj,  tl-.  .  .  .  ,tpj,  §lr  ■  ■  ■  . 

In  the  matrix  notation  the  cost-minimization  problem  is  then  to 
minimize 

ex 

subject  to 

AX  .  b 

where 

[a 

and 

C  -  l  cl . (n  *m  *p  ♦ 
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The  Cj’s  are  the  costs  associated  with  the  chains  (j  =  1, the  slack 
variables  (j--n  +  l,...,it+m+p,  and  the  artificial  variables  (j  =  .1  +  m  + 
p  +  1, .  .  .  ,n  +  m  +  p  +  q  ). 

Adopting  the  standard  two-phast  tolution  procedure  for  the  minimum- 
cost  formulation,  the  Phase  I  procedure  minimizes  to  zero  f  he  value  of 

f  =  <fe) 

jon+m+p+1 


to  obtain  an  initial  basic  feasible  solution.  This  effectively  assigns  a  cost  of 
1  to  the  artificial  variables  and  a  cost  of  zero  to  the  other  variables  in  Phase  I. 
Phase  II  begins  with  tee  basic  feasible  solution  determined  in  Phase  I  and  pro¬ 
ceeds  to  minimize 


DISCUSSION 

Substitution  of  Resources 

In  the  previous  formulation,  for  each  arc-commodity  pair  a  single  com- 

*  ^  y 

bination  of  resources  is  required  and  represented  by  the  vector,  R"  =  (ra  ,  rj2  , 

.  .  ,  ,  rjp ) ,  of  the  matrix  R1* . 

Define  a  new  resource  matrix:  T  =  [tjk]  (i  »  1,  k  ■  1,  ...  ,4) 

A  L  A  L 

where  tjfc  =  (Rj  |  R"  is  any  feasible  resource  vector  for  arc  i ,  commodity  It); 
Rj  =  (ra  1  •  •  •  )• 

In  words,  each  element  of  T  is  the  set  of  alternative  resource  vectors 
for  a  movement  of  one  unit  of  commodity  k  over  arc  i . 

The  contribution  of  each  arc  to  Cj  -  cB  B"1  A.  in  Phase  II  of  the  minimum- 
cost  procedure  is 

d,k  -  rt  -  a,  +  Us  -  ns) 

The  possiuiiity  of  employing  alternative  methods,  i.e>,  alternative  combinations 
of  resources,  affects  this  by  allowing  for  a  number  of  vectors  Rf .  Thus  to  find 
the  minimum  c,  -  cB  B-1  A(  one  must  find  the  minimum 
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over  the  permissible  Rk  as  "ell  as  over  all  feasible  combinations  of  arcs.  The 

elements  0S  (s  =  1 . p)  are  fixed  for  any  problem  and  the  elements  irs  (s  = 

1, .  . .  ,p>  arc*  fixed  for  any  iteration.  One  may  therefore  find  the  vector 

-k  -k min  r  v  .  i 

R/  ‘  R. £  ‘(it  ~n*)j  for  ‘  -1*  •  •«; k  -1 . « 

The  kth  matrix  of  these  minima  may  then  be  defined  as: 

R*  =  [rf,]  u  -  1 . .  s  n  1,  .  .  . ,  p) 

Each  column  vector  Rk  *  (F,\  , .  .  .  ,  rfk  )  is  the  alternative  combination  of  re¬ 
sources  such  that 

is  minimized  for  arc  i  and  commodity  k.  Now  Rk  may  be  substituted  in  the 
minimum-cost  procedure  previously  discussed  and  the  appropriate  A,  selected 
for  entry  into  the  basis.  Thus  new  Rfc  is  constructed  for  each  commodkj 
iteration. 


Summary  oi  the  Procedure 


To  summarize,  the  proposed  procedure  is: 

( 1)  Calculate  CgB  *  t  oe  i  ,  .  .  .  ,  o  ^ ,  .  .  .  ,  a  $ ,  Oi ,  .  .  .  ,  ) 

(2)  For  each  arc -commodity  pair  find  tht-  least-cost  applicable  resource 
vector,  “cost*  meaning  cost  in  terms  of  the  simplex  multipliers  and  resource 
prices. 


min  i-  p  -i 

[,*/» (*. -».)J 


(3)  For  commodity  k  =»  1, ...  ,q  calculate 


dk'  - 


a  + 


(*3  ~  {°r  ‘ 


,  m 


and  assign  the  dj  to  the  arc  i  as  a  pseudo  cost. 

(4)  Using  the  shortest -chain  algorithm  find  the  chain  with  least 


di  ‘  ,§» °'J  [r<  ~  a>  *  ,*i7*»  ‘  "‘’j  ,pr 


k  -l, 
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(5)  Find 


(6)  If  the  minimum  [  d*  -  <yfc]  <  0,  the  vector  Af  is  entered  in  the  basis. 
If  [d*  -  ak]  s  0,  there  is  no  chain  that  may  improve  the  value  of  the  objective 
function,  and  the  procedure  is  terminated. 


Validity  of  the  Procedurr 


Consider  the  linear  programming  formulation  of  the  substitution  problem. 
It  is  identical  to  the  original  cost-minimization  problem  except  that  every 
column  vector  in  the  criminal  problem  will  be  replaced  by 


II  N(M,)  [wh;re  N(M,)  ia  the  number  of  alternate  reiourte 
“ i  *  vcctora  applying  to  arc  i  ] 


alternate  chains.  The  expanded  substitution  matrix  will  be  many  limes  larger 
than  the  original  matrix,  should  either  actually  be  enumerated. 

It  is  claimed  that  the  procedure  outlined  here  will  find  the  least-cost  (in 
the  sense  previously  described)  vector  to  enter  the  basis.  It  should  be  noted 
that  if  the  procedure  does  not  find  the  least-cost  vector,  but  some  ether  vector, 
sa;  the  nth  least-cost  vector,  the  algorithm  will  progress  toward  an  optimum 
solution  in  the  early  stages  but  will  terminate  early.  That  is,  any  vector  that 
satisfies  the  simplex  rule  may  be  brought  into  the  basis,  but  since  the  algo¬ 
rithm  is  terminated  when  the  “shortest”  chain  does  not  satisfy  the  simplex 
rule,  the  validity  of  the  procedure  depends  on  the  validity  of  the  shortest-chain 
procedure. 

Suppose  that  the  chain  produced  as  a  candidate  is  not  the  shortest  chain 
and  there  is  some  other  candidate  chain  j  *  for  which 

M'  s  A*' 

dJ.  -  •». <  di  -  °k 

Two  possibilities  for  this  other  chain  exist: 

(1)  The  shorter  chain  consists  of  the  same  arcs  as  our  candidate  chain 
but  has  different  (allowable)  resource  vectors  associated  with  one  or  more  of 
these  arcs. 
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(2)  The  shorter  chain  consists  of  different  arcs  altogether  with  some 
allowable  set  of  resource  vectors  assigned  to  their  respective  arcs. 

The  first  case  In  a  chain  that 

di*  <  ‘i 

or  that 

m  m  f  *k 

,5Ia«(,l“ai)  +  ,*,-U  ~ak 

Is  less  than 

.SjV'i-V  +  J,flo 

but,  since  the  first  and  last  terms  of  each  expression  are  Identical,  that  Is  a 
claim  that  for  at  least  one  arc,  common  to  both  chains, 


P  u 

but  since  r's  ts)  (j  »  1,  ...,*;  k  =>  1, .  .  ,,q)  Is  the  minimum  e.vall- 

t’jle  (step  2),  the  claim  that  d£  -  vk,  <  d*’  -  oj,  is  inconsistent,  and  hence 
cate  1  cannot  occur.  A  true  shortest  chain  must  employ  the  least-cost  allow¬ 
able  resources  on  each  arc  that  Is  a  member  of  the  chain. 

Case  2  resolves  Itself  to  a  claim  that  there  is  some  chain  ‘''.at  uses  the 
least-cost  allowable  resources  on  each  of  its  member  arcs  and  has  a  lower 
( -  okJ  than  that  of  the  candidate  chain.  Sine*,  the  proposed  procedure 
evaluates  the  pseudo  cost  of  eich  arc  incorporating  the  minimum  resource 
costs  [i.e.,  L  f*  '$,-*,)]  and  identical  arc -use  pseudo  costs  [i.e., 

M  S  “  1 

L  ( t(  -  at, )  ] ,  a  claim  that  case  2  exists  is  simply  a  claim  that  the 
shortest-chain  algorithm  does  not  find  the  shortest  chain. 

Usefulness  of  the  Procedure 

In  order  to  be  useful  in  application,  the  routes  selected  and  resources 
assigned  must  be  feasible  In  the  object  system.  Routes  through  the  network 
are  composed  ot  a  series  of  arcs  and  the  resources  assigned  to  them.  Again 
using  a  transportation  network  as  an  example,  It  is  undesirable  to  have  dif¬ 
ferent  vehicle  types  assigned  to  contiguous  arcs  of  the  same  mode  In  a  chain. 
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That  is,  one  wants  the  same  vehicle  to  carry  the  commodity  over  all  contiguous 
arcs  of  the  same  mode  in  a  chain.  Quarter-ton  and  12-ton  trucks  may  be  feasible 
subctitutes,  but  one  does  not  wish  to  transfer  from  one  to  another  at  a  node. 

This  is  an  Important  consideration  if  the  results  of  the  solution  are  to  be 
used.  It  is  simply  not  feasible  in  practice  to  use  chains  that  employ  different 
vehicles  on  various  arcs  of  the  same  chain  unless  th*  chain  is  multimode  and 
transfer  arcs  are  included.  A  procedure  that  is  computationally  simpler  than 
the  general  method  just  described  is  available,  and  it  guarantees  that  the  same 
resource  combinations  will  be  used  on  all  arcs  of  a  chain  that  are  of  a  particular 
mode. 

A  "master*  resource  vector  representing  the  resources  required  to  sus¬ 
tain  a  unit  flow  over  a  standard  arc  of  unit  length  is  provided  for  every  com¬ 
modity,  mode,  '.id  method.  Each  arc  then  has  a  mode  identifier,  a  condition 
factor,  and  a  length  factor  assigned.  The  minimum-cost  (in  terms  of  the  sim¬ 
plex  multipliers)  method  is  then  selected  for  each  iteration,  and  the  resource 
vectors  for  each  arc  are  generated  using  the  condition  and  length  scalars. 

Thus  a  single  master  vector,  representing  a  particular  method,  is  selected  as 
the  minimum-cost  method  for  all  arcs  of  that  mode  for  each  Iteration.  Con¬ 
tinuity  of  vehicle  type  is  ensured  for  all  chains  in  the  solution. 


COMPUTATIONAL  EXPERIENCE 

A  computer  program  in  FORTRAN  IV  for  the  Control  Data  6400  has  been 
developed  for  both  maximum-flow  and  minimum-cost  formulations  incorporating 
the  substitution  feature.  The  program  uses  the  product  form  of  the  Inverse  and 
will  accommodate  up  to  150  commodities,  550  arcs,  and  50  resources.  Up  to 
10  modes  are  permitted  and  each  mode  may  have  up  to  three  alternative  re¬ 
source-requirement  vectors.  Thus  each  arc  may  use  any  of  three  feasible 
combinations  of  resources  to  accomplish  the  move.  A  series  of  test  cases  has 
been  solved  successfully  and  the  results  are  encouraging  with  respect  to  ac¬ 
curacy  and  speed  of  solution.  The  use  of  the  substitution  feature  does  increase 
the  time  required  tor  solution,  but  this  Increase  has  been  small  in  the  ccses 
tested  to  date. 
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ABSTRACT 


The  minimum-cost  formulation  of  the  problem  of  determining  multi- 
commodity  flows  over  a  capacitated  network  subject  to  resource  con¬ 
straints  has  been  treated  in  previous  papers.  In  those  treatments 
only  capacitated  arcs  were  assumed  and  a  uniform  unit  of  measure  like 
short  tons  was  used  for  all  commodities.  This  paper  treats  the  effect 
of  constraints  on  the  nodes  of  the  network,  allows  the  commodities  to 
be  measured  in  their  "natural"  units  and  allows  the  network  capacities 
to  be  expressed  in  vehicles  per  time  peidod  —  in  some  cases  giving  a 
more  accurate  representation  of  the  capacities  of  the  network.  The 
solution  procedure  using  the  column  generation  technique  is  described 
and  computational  experience  is  discussed. 
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INTRODUCTION 


In  previous  papers1  the  maximum  flow  multicommodity  network  pro- 

2 

blem  as  formulated  by  Ford  and  Fulkerson  and  the  corresponding  minimum 

3 

cost  problem  as  stated  by  Tomlin  were  extended  to  include  resource 
constraints  and  resource  substitution.  The  purpose  of  this  paper  is 
to  describe  other  modifications  to  the  multicommodity  network  flow 
model  that  have  been  found  to  be  useful  in  studies  of  transportation 
networks. 

These  modifications  will  be  described  in  two  steps.  First,  the 
problem  will  be  revised  to  permit  the  expression  of  cornmo  .ity  flows  in 
terms  of  the  natural  unit  of  the  commodity  and  the  expression  of  arc 
capacities  in  terms  of  vehicles  per  time  period.  In  the  standard 
formulation,  flows  and  requirements  for  commodity  movement  a’re  con¬ 
verted  in  into  some  common  unit  such  as  3hori  tons  for  processing  by 
the  algorithm.  It  is  desirable,  both  for  convenience  and  accuracy, 
to  represent  all  commodities  in  terms  of  their  natural  units  (e.g., 
numbers  of  passengers,  barrels  of  gasoline,  tons  of  coal,  etc.).  It 
is  also  convenient  to  express  arc  capacities  in  terms  of  vehicles  per 
time  period  since  most  traffic  statistics  are  collected  in  this  fashion. 

Second,  node  flow  constraints  will  be  added  to  the  problem  formu¬ 
lation.  In  the  analysis  of  transportation  networks,  there  are  fre¬ 
quently  cases  where  the  intersection  of  three  or  more  arcs  has  itself 
a  capacity.  The  capacity  of  a  node  is  not  easily  expressed  through 
the  addition  of  dummy  arcs.  It  will  be  shown  that  node  constraints 
can  be  added  to  permit  explicit  statement  of  these  capacities. 
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MINIMUM  COST  MULTICOMMODITY  NETWORK  FLOWS  WITH  RESOURCE  CONSTRAINTS 


For  sake  of  continuity,  the  minimum  cost  problem  with  resource 
constraints  is  restated. 

Consider  the  multimode,  multicommodity  network  G(N,  G).  N  is  the 
set  o *  all  nodes  of  the  network.  G  is  the  subset  of  all  ordered  pairs 
(x,  y)  of  the  elements  of  N  that  are  arcs  of  the  network.  G^,  G^, 

. ..,  (^  is  an  enumeration  of  the  arcs.  Each  arc  has  an  associated 
cost  (or  distance)  d(x,  y)^0  and  each  arc  has  an  associated  capacity 
b(x,  y)aO. 

For  each  commodity  k(k  =1,  . ..,  q)  there  is  a  source  s^  and  a 
sink  t^.  The  flow  of  the  commodity  k  over  arc  (x,  y)  is  denoted  by 
F^x,  y),  k  =  1,  ...,  q.  The  flows  and  the  arc  capacities  must  be 
stated  in  terms  of  some  common  unit  such  as  short  tons.  These  flows 
must  satisfy  the  capacity  constraints. 


Z  F*(x,  y)  £  b(x,  y),  (x,  y)€G 

k=l 


In  keeping  with  previous  notation,  define  the  set  =  -jP 


00 


(k)  1  "  I 

P:  '  is  a  chain  connecting  s,  and  t,  k  Now  let  P  be  the  union  of  the 

j  k  kJ  /i  \  /i  \  /,.\ 

sets  r  (k  =  1,  ...,  q).  Further,  let  P^ 


(<l) 


be  the  enumeration  of  the  chains  P, '  'GP  such  that  the  subscript 


(1)  „(1) 
00, 


y  Pg 


Lj 


(k) 


n  3 

j  is  sufficient  to  identify  the  chain,  its  origin-destination  pair,  and 


the  commodity  with  which  it  is  associated. 

Thus  the  k1"*1  commodity  set  is  defined  by 
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00 


is  a  chain  from  sfc  to  t. 


k  * 


q. 


The  arc-chain  incidence  matrix  is 


A  .  [a13l 


where 


ij 


-{ 


X,  If  QjSPj 
0,  otherwise 


00 


for  1  J  —  • • • j  n • 

Define  the  resource  matrix  for  commodity  k  by 
Rk  ■«  [rjg]  (i  =  1,  ...,  a;  s  =  1,  ...,  p) 
where  r^s  is  the  quantity  of  resource  s  required  to  sustain  a  unit  of 
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It  .  » 

flow  of  commodity  k  over  arc  i;  r,  so.  And  latp  (s  =  1,  . ..,  p)  be 

».s  s 

the  quantity  of  resource  s  available  (e.g.,  in  inventory)  for  assign¬ 
ment  to  the  network. 

The  objective  is  to  minimize  the  total  cost  of  the  resources  used 
end  the  tolls  on  the  arcs  while  meeting  delivery  requirements  and 
satisfying  the  arc  capacity  and  resource  inventory  constraints.  Lett¬ 
ing  x^k^(j  =  1,  ....  n)  be  the  flow  of  commodity  k  in  chain 
(j  =  1,  . ..,  n;  k  implicit),  b^  the  flow  capacity  of  in  arc-chain 
terms,  the  minimum  cost  network  flow  linear  program  with  resource 
constraints  is: 


Minimize  n 


E  cVk^  =  z 
j=l  J  J 


subject  to 

(a)  arc  capacity  constraints 


n 


(k)  <  bt,  i  -  1, 


A  « J 

(b)  resource  constraints 

m  q  i,  00 


E  E  E 
i=l  k=l  j€J. 


,  k 
a.  4r4  x  ^ 
ij  is  d 


m 


S  —  1  j  •  •  •  y  p 


and 


(c )  delivery  requirements 


.00  . 


J 


”  V  k  -  1, 


where  X^.  is  the  delivery  requirement  at  t^  ,  X^C,  (k  =  1,  ...,  q). 

The  cost  coefficient  c.  may  be  defined  as: 
m  Jp  m  . 

Ci  =  S  Tiai1  +  1  2  Vis  V  (for  j  =  1}  n;  k' 

(k)  J  1=1  ^  6=1  1=1  J 

where  P^  connects  s,  and  t,  )  and 

J  K  K 

where  t  is  the  cost  (or  toll)  for  a  unit  flow  over  arc  i,  and  $  is 
i  s 

the  cost  of  using  a  unit  of  resource  s* 
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THE  REVISED  PROGRAMMING  PROBLEM 


Use  of  Natural  Units  and  Vehi ^les  Per  Time  Period 

The  capacity  b(x,  y)  of  arc  (a,  y)@G  will  now  be  expressed  in 
vehicles-per-time-period.  Sane  arcs  will  have  vehicle  types  naturally 
associated  with  them.  Highway,  rail,  waterway  and  airway  arcs  will  be 
naturally  associated  with  trucks  or  busses,  rail  cars,  barges  and  air¬ 
craft  respectively.  For  other  modes  such  as  pipelines  and  transfer 
the  "vehicle"  must  be  contrived.  One  alternative  is  to  use  a  dummy 
vehicle  which  contains  exactly  one  unit  of  the  commodity. 

It  should  be  stressed  at  this  point  that  all  commodities  will  be 
counted  in  their  natural  unit.  That  is  petroleum  products  will  be  ex¬ 
pressed  in  barrels,  passengers  as  numbers  of  passengers  and  bulk  cargo 
in  tons  if  these  are  the  units  best  suited  to  them.  Throughout  the 
remainder  of  this  paper  each  reference  to  flow  or  movement  of  a 
commodity  will  be  in  the  natural  unit  of  that  commodity. 

To  expresr  the  flow  Fk(x,  y)  in  the  natural  unit  of  commodity  k 
and  the  arc  constraints  in  vehiele-per-time-period,  re-write  the 
capacity  constraints  as  follows: 

E  *  b(x,  y),  (x,  y)02; 

k=l  u  (x,  y) 

where  Ax,  y)  is  the  quantity  of  the  commodity  k  transported  over  arc 
(x,  y)  in  a  single  vehicle,  u  (x,  y)>0;  and  b(x,  y)  is  now  in  vehicles 
per  time  period. 

Define  a  matrix  of  vehicle  per  commodity  unit  coefficients  as: 

V  =  [Vik] 

where  , 

Ir  _ ^  th 

v^  =  yj  (the  arc  (x,  y)  is  denoted  as  the  i  arc 

by  the  index  i). 

k 

In  words,  v^  is  the  reciprocal  of  the  quantity  of  commodity  k  trans¬ 
ported  over  arc  i  in  a  single  vehicle. 

With  the  above  change,  the  minimum  cost  linear  program  can  now 
be  stated  as  follows : 


Minimize 


z 


Node  Constraints 

Node  constraints  allow  a  direct  limitation  on  the  throughput 
capacity  of  a  node  without  the  necessity  of  additional  assumptions  in 
depicting  the  node  graphically  as  a  small  subnetwork.  Examples  may 
clarify  this  point. 

In  figure  1,  the  node  given  in  (a)  is  to  be  constrained  by  re¬ 
placing  it  with  the  subnetwork  given  in  (b).  The  capacity  limitation 


on  arc  5  is  that  of  the  node  xt  has  replaced.  Note,  however,  that  the 
path  1-4  in  (a)  contributes  to  the  flow  through  the  node  but  in  (b) 
it  does  not  contribute  to  the  flow  on  arc  5.  In  some  cases,  this  pro¬ 
blem  may  be  solved  by  directing  the  arcs,  ^.e.,  by  permitting  flow  in 
only  one  direction  as  in  figure  1-c.  This  method  is  not  always  satis¬ 
factory  because  in  a  multi-comnodity  problem  the  optimum  solution  may 
require  that  one  commodity  flow  over  path  1-4  while  another  is  required 
to  flow  over  path  2-4. 
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A  representation  such  as  that  shown  in  figure  2  is  sometimes 
attempted,  but  this  too  is  unsatisfactory.  What  capacity  should  be 
assigned  to  arcs  5>  6,  7,  and  8? 


Figure  2 

To  get  an  accurate  node  throughput  limitation,  a  constraint  row 

is  added  in  the  problem  formulation  for  each  node  for  which  a  capacity 

A 

limitation  is  desired.  Let  N  be  the  subset  of  N  of  those  nodes  which 

are  to  be  constrained.  To  consider  the  node  throughput  limitations, 

a  A  .  . 

let  N^,  ...,  be  an  enumeration  of  the  nodes  N^GJ.  Let  b  (N^)  be 
the  capacity  associated  with  the  node  N  which  is  expressed  in  vehicles 
per  time  period. 

Define  a  matrix  of  vehicle  per  commodity  unit  coefficients  for 
node  constraints  similar  to  those  for  the  arc  constraints: 

V  =  [$/] 

where 

A  k  1 

and  ^(N^)  is  the  quantity  of  the  commodity  k  transported  through  node 

in  a  single  vehicle,  U^(N^ )X).  Let  be  the  throughput  capacity 

of  node  N. . 

I 

The  node-chain  incidence  matrix  is  F  =  [f. .1 

x  j 


. r*-"- ™,->  .(kpsr. •■■■  •' 


where 


it 


-  f’  lf  HX€PJ 


00 


0,  otherwise 
f  Or  jfc  =  1 ,  .  •  •  ,  t  J  j  =  1  ,  .  »  .  ,  11  • 

In  arc-chain  terms,  the  minimum  cost  network  flow  linear  program 
with  resource  constraints,  node  constraints  and  putting  the  capacity 
constraints  in  vehicles  per  time  period  is  as  follows : 

Minimize 

l  cxto  -z 

j=l  J  J 


m 


subject  to 

(a)  arc  capacity  c  Distraints 

2  a  v  k  x(‘r^  £  b  ,  i  =  1,  . 
j=l  1  J  1 

(b)  node  capacity  constraints 

f«  V  xj(k)  s  v 1  ■ i: . * 

(c)  resource  constraints 

m  q  k  (k) 

I  2  2  a  r  x  s  p  ,  s  =  1,  ...,  p 

i=l  k=l  J€J,  18  J  8 


and 


(d)  delivery  requirements 


J€J, 


xjk)  =  xv,  k  =  1, 

j  k 


SOLUTION  PROCEDURE  USING  THE  COLUMN  GENERATION  TECHNIQUE 

We  have  defined  a  rtraight  forward  linear  programming  problem 
different  from  that  previously  defined  only  in  that  node  constraints 
have  been  added,  the  natural  units  are  used  for  the  flows  and  the 
delivery  requirements,  and  vehicles  per  time  period  are  used  as  capa¬ 
cities.  It  is  now  necessary  to  show  that  the  column  generation 
technique  can  be  adapted  to  these  changes.  This  technique  is  more 
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fully  described  in  reference  1. 

Define  the  matrices  A*,  E,  F*  and  G  as  follows : 

A*  -  '“h1 


where 


where 


a*  =  a  v^,  k  Implied  t!y  Pj  . 

E  -  [eBjl 

e  .  =  Z  r*  a.  k  implied  by  • 
sj  i=1  is  i«j  J 

F*  =  [f*J 


where 


and 


where 


fh  =  '  k  *  PJ 


00 


0  “  [skJ] 


1,  if  j€J 
k 

otherwise 


A  is  a  matrix  (m  +  p  +  t  +  q  x  n)  formed  of  the  sutmatrices  A*,  E,  ?*, 
and  G  as  follows: 


A 

A  = 


A* 
E  _ 
F*_ 

G 


%■  Hr 


The  typical  column  of  A  is 

=  col.  (afy  a*j,  e1Jf  epj,  ffy  . 

gqj’- 

A  will  be  quite  large  due  to  the  fact  that  there  is  typically  an 
enormous  number  of  chains  connecting  s^  and  t^  (k  =  1,  q).  The 

shortest  chain  algorithm  can  be  used  to  develop  the  Aj  that  will 
satisfy  the  simplex  rule.  Further  if  the  shortest  chain  algorithm 
can  find  no  chain  satisfying  the  requirement,  an  optimum  has  been 
reached. 

This  formulation  can  be  solved  by  adopting  the  standard  two-phased 


g 


8 


4 


I 


„  -.v.**.**^^  ?-  wwf-*1 '  ■  a-nss*?  Ae^tfiwr1*-  -iffa--  ■  .•< -->vf  w  -■ 


O 


o 


procedure.  Phase  I  minimizes  to  zero  the  value  of 

j=n+m  +  t+  p+  q 

Zx  (k) 

1  Xj 

j=n+m+t+p+l 

to  obtain  an  initial  basic  feasible  solution.  This  effectively  assigns 
a  cost  of  1  to  the  artificial  variables  and  a  cost  of  zero  to  the  other 
variables  in  Phase  I.  Phase  II  begins  with  the  basic  feasible  solution 
determined  in  Phase  I  and  proceeds  to  minimize 


n 

L 


n  x  GO  -  Z 

~J  J  ' 


j=l 

In  Phase  I,  I  ,  may  be  used  as  the  initial  basis  and  the  sim- 
’  m+p-tq 

plex  rule  is  to  enter  a  chain  in  the  basis  if,  and  only  if, 

YV’V0 


where 


cbb_1  =  (V 


■'  V 


■>  V  h> 


’>  V  av 


">  CTq) 


so  that  the  simplex  multipliers  a ^  are  associated  with  the  arcs,  the 

ira  are  associated  with  the  resources,  the  0^  are  associated  wi'-h  the 

nodes,  and  the  a.  are  associated  with  the  artificial  variables.  Thus 
A  * 

the  vector  is  entered  if 


m 

-I  a 
i=l 


ij 


Vi*  +  \  ns  ris 


t 

Lj 


K  Bx  <  v 


The  contribution  of  each  arc  to  c,-cJB*  A. 

J  B_  J 


1=1 

in  Phase  I  is 


=  - 


Vi 


8=1 


k 

tt  r. 
s  is 


;  noting  that  c,  =  0  for  struc¬ 


tural  vectors. 

The  contribution  of  each  node  to  c.-c^B-1^.  must  be  distributed 

J  B  j 

over  the  arcs  to  .acilitate  the  use  of  the  shortest  chain  algorithm. 
Disallowing  the  origin  and  destination  nodes  from  being  constrained, 
it  is  clear  that  if  a  constrained  node  is  contained  in  a  route,  the 
route  also  contains  precisely  two  of  the  arcs  which  emanate  from  or 
terminate  at  this  node.  !Hiis  suggests  a  method  for  the  distribution 
of  the  node  contributions.  Namely,  that  half  of  the  contribution  be 
added  to  each  arc  which  has  this  node  as  an  end  point. 
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To  do  this,  define 

if  node  N.  is  an  end  point  of  axe  G. 


6  „  =  f  ’ - 1 

0,  otherwise 


.-14 


The  effective  contribution  of  each  arc  to  c.-c_B  A.  becomes 

j  u  j 


Vi  *  _Wu  +  *s  6 


S=1 


1=1 


A  k 
v  9 
4i  4  -4 


The  shortest  chain  algorithm  may  be  used  to  find 

r 


min 

J 


L«i®J 


00  ^ 


=  min 


_  /  k  ^  k  l  ~  ^ 


over  all  k. 

Where  G.  is  the  ith  arc  and  P^)  is  the  chain  from  s.  to  t,  . 
i  j  k  k 

The  minimum  over  all  commodities  is  selected  as  the  candidate  to  enter 
the  bails.  The  column  vector  to  leave  the  basis  may  be  determined  in 
the  standard  simplex  fashion.  Should  any  a . ,  tt  or  0  be  positive  the 
corresponding  slack  variable  may  be  entered  into  the  basis. 

In  Phase  II 
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m 
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6  k„ 


C4-CJ3  A  =  E  T.a  +  E  le  .  ■  Sv.  a,a. ,  -  E  Tte  ,  -  E  f.,v.  p 


'i  V  j 

X 


i=l 


m  V  *  ***  v.  u  .  v 

=  i=l  &i^Ti"Vi  +  Ji  ^ s  "  fA3V4  PX 


i  U 


:r«  sj  1^1i  i  ij  Sa!  s  sj  “,4J'4H4 


s=l 


m 


4f 1 


* 


£ 


£  ^ 

where 


and 


Thus 


m  k 

e  =  E  rf"  a . 
sj  is  ij 


=  f ■  if 

l0,  otherwise 


,K  K.  .  ~  X.  \ 

di  ‘  T1  -  T!  “l  +  Z,  rlS  (V"s> 

8=1 
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nMnum? <  w  r 


E  j 
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and 


^  =  Ti  -  VA  +  \  ris<VV  ”  i  ?,  SH  V*i 

S— x  x~x 

may  be  assigned  to  arc  i.  The  shortest,  i.e.,  the  chain  with  the  least 

E  £  -  *  <  0 
1=1  K 

for  k  =  1,  . . . ,  q,  may  then  be  entered  in  the  basis. 

Phase  n  is  terminated  and  the  value  of  z  is  minimized  when,  for 
the  minimum  j, 

% 

S  (k)  di  *  °k* 

Gi€Y  ; 

MODIFICATION  OF  THE  SUBSTITUTION  PROCEDURE 

The  substitution  procedure  is  affected  by  the  change  in  that  the 

use  of  a  different  group  of  resources  may  require  that  a  different 

coefficient  of  vehicles  per  commodity  unit  be  us  id.  Under  the  procedure 

described  above  the  matrix  V  is  a  two  dimensional  matrix  by  arc  and 

by  commodity.  When  substitution  is  permitted  this  matrix  must  be  given 

the  additional  dimension  of  substitution  group. 

a 

Redefine  the  matrices  V  and  V,  the  matrices  of  vehicles  per 
commodity  unit,  as: 

V  =  [vfh]  and  V  =  [v^]  respectively 
k 

where  is  the  reciprocal  of  the  quantity  of  commodity  k  transported 

over  arc  i  in  a  single  vehicle  determined  by  substitution  group  h  and 

is  the  reciprocal  of  the  quantity  of  commodity  k  transported 

through  node  A  in  a  single  vehicle  determined  by  substitution  group  h. 

The  coefficients  and  are  therefore  selected  in  each  iteration 
ih  £h 

based  on  the  minimum  "cost"  (in  terms  of  simplex  multipliers)  alter¬ 
native  resource  vector  to  be  used. 

Summary  of  the  procedure 

To  summarize,  the  proposed  procedure  is: 
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(1)  Calculate  C^B  -  (q^>  •  ••,  tt^,  •••»  TTg.«  •••> 

**]_>  •  ••,  o^) 

(2)  In  the  previous  formulation,  for  each  sure- commodity  pair 
a  single  combination  of  resources  is  required  and  represented  by  the 


vector, 


P*  =  (rfp  r\2,  ...,  rjp),  of  the  matrix  Rk. 


Define  a  new  resource  matrix:  T  =  [t^]  (i  =  1,  ...,  m;  k  =  1, 
...,  q)  where  t^  =  {Rk  J  is  any  feasible  resource  vector  for  arc  i, 

commodity  kj-  Rk  =  (r^,  ••.,  £kp). 

In  words,  each  element  of  T  is  the  set  of  alternative  resource 
vectors  for  a  movement  of  one  unit  of  commodity  k  over  arc  i.  For 
each  arc- commodity  pair  find  the  least- cost  applicable  resource  vector, 
"cost"  meaning  cost  in  terms  of  the  simplex  multipliers  and  resource 
prices  and  call  it  the  vector  h. 

5k  fAk  /vk  / 1  _  \  j  _  _  j _ _ Ak  r.  1 


“ih  =  R  £  >  u  minlnnM  for  ^  ««} 

S=i 


where 


Rk  -  ;k ) 

lh  "  ^ril'  ***'  rip; 

(3)  For  commodity  k  =  1,  . ..,  q  calculate 

%  -  Ti  ’  Vto  +  \  <V"s>  -  i,\  6ii  "to  = 

S=1  x=i 

•  •  • )  in 

and  assign  ctk  to  arc  i  as  a  pseudo  cost. 

(4)  Using  the  shortest  chain  algorithm,  find  the  chain  with  the 

least 

m  %  P  — k  ^  Air 

a ,  =  I  a  T . -Of. v~  +  2  r  ($  -tt  )  -  \  2  6..  V~  B  for  k  =  1,  ..., 

J  1=1  ij  i  i  ih  fl:5l  is  s  s'  z  ^  Ji  jgh  Ki 

(?)  Find  min  [lk  -  o^] 
k 

(6)  If  the  minimum  [clk  -  o.]<0,  the  vector  is  entered  in  the 

J  k  j 

basis.  If  minimum  [clk  -  oj*0  there  is  no  chain  that  may  improve  the 

J  it 

value  of  the  objective  function,  and  the  procedure  is  terminated. 
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CONCLUSIONS 


It  seems  to  be  the  consensus  of  the  literature  that  the  basic 
Ford-Fulkerson  multicommodity  algorithm  attains  its  speed  based  in 
part  upon  the  fact  that  the  conventional  method  result s  in  a  zero- one 
matrix.  It  should  be  emphasized,  however,  that  the  literature  contains 
no  direct  test  of  this  idea.  The  assumption  of  3peed  seems  to  be  based 
largely  on  the  simplicity  of  the  zero-one  matrix  and  the  fact  that  the 
decomposition  algorithm  (that  is  directly  related  to  the  Ford -Fulkerson 
suggested  method)  typically  does  require  large  amounts  of  computer  time. 

The  Control  Data  6400  computer  program  for  the  RACA'f  algorithm 
which  had  been  written  at  the  Research  Analysis  Corporation  was  revised 
to  include  the  changes  discussed  in  this  paper.  It  was  of  interest  to 
make  a  comparison  of  the  performances  of  the  unrevised  and  the  revised 
RAC  AT  algorithms.  Two  problems  were  used  as  tests.  One,  a  small  pro¬ 
blem  with  52  arcs,  5  resources,  2  origin-destination  pairs  and  no  node 
constraints,  a  total  of  59  rows.  The  computer  time  statistics  on  the 
CDC  6400  of  two  variations  of  the  problem  formulation  are  given  in  . 
Table  1  below.  One  formulation  has  the  arc  capacities,  the  movement 
requirements  and  the  route  flows  in  terms  of  tons  per  day.  The  other 
has  the  arc  capacities  in  vehicles  per  day,  the  movement  requirement 
and  the  route  flows  in  the  commodities  *  natural  units. 

Table  1 

COMPUTER  TIME  STATISTICS  FOR  PROBLEM  1 


Central  processor 

No.  of 

Time  per 

Constraint  type 

time  (sec) 

iterations 

iteration  (sec) 

Vehicles  per  day 

13.3 

19 

.70 

Tons  per  day 

11.8 

23 

•  51 

In  this  simple  problem  one  arc  constraint  type  was  directly  con¬ 
verted  from  the  other  to  make  the  problems  as  much  the  same  as  possible. 
As  a  result,  the  objective  function  values  were  the  same  and  the  flows 
over  respective  arcs  were  the  same.  Every  effort  was  made  to  keep  the 
problems  similar  so  that  the  timing  comparisons  would  be  meaningful. 

The  central  processor  time  for  the  run  with  the  arc  constraints  in 
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vehicles  per  day  was  12. 7$  greater  than  that  for  the  run  with  the  arc 
constraints  in  tons  per  day. 

The  second  problem  which  was  used  as  a  test  did  not  have  the 
exact  conversion  from  the  tons  formulation  to  the  vehicles  formulation 
as  the  first  one  had.*  In  the  vehicles  formulation  two  arcs  of  the 
network  were  combined  into  one  with  the  capacity  being  adjusted  accord¬ 
ingly  and  several  resources  with  identical  productivities  were  combined 
and  the  resource  inventories  were  adjusted  accordingly.  The  tons  formu¬ 
lation  of  the  problem  has  290  arcs,  42  resources,  9  origin  destination 
pairs  (2  of  them  had  zero  requirements ),  a  total  of  3^1  rows.  The 
vehicles  formulation  has  the  same  statistics  except  for  one  less  arc 
and  five  less  resources.  Table  2  gives  the  run  time  comparisons  on 
the  CDC  6400.  As  a  result  of  the  inexact  conversion  of  the  tons  formu¬ 
lation  to  the  vehicles  formulation,  the  objective  function  for  the 
vehicles  formulation  was  slightly  smaller,  as  expected,  end  the  solution 
was  slightly  different  from  the  tons  formulation.  The  central  pro¬ 
cessor  time  for  the  run  with  the  arc  constraints  in  vehicles  per  day 
was  55$  greater  than  that  for  the  run  with  the  arc  constraints  in  tons 
per  day. 

Table  2 

COMPUTER  TIME  STATISTICS  FOR  PROBLEM  2 


Central  processor 

No.  of 

Time  per 

Constraint  type 

time  (sec) 

iterations 

iteration  (sec) 

Vehicles  per  day 

i860 

265 

7-0 

Tons  per  day 

1200 

156 

7.7 

It  has  been  suggested  in  the  past  that  the  multicommodity  network 
flow  problem  with  resource  constraints  might  best  be  solved  using  the 
decomposition  method  of  Dantzlg  and  Wolfe.  With  the  addition  of  the 
V  matrix  and  the  destruction  of  the  zero- one  nature  of  the  major  sub¬ 
matrix,  it  could  now  be  argued  more  Btrongly  that  decomposition  might 


*This  problem  was  solved  for  purposes  other  than  that  of  testing 
the  algorithm  and  the  time  to  solve. 
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><e  u  better  approach.  It  is  possible  that  this  is  true.  No  systematic 
evaluation  of  the  decomposition  alternative  has  been  made  in  the  course 
of  the  study  leading  to  this  paper.  It  is  argued,  however,  that  from 
a  practical  standpoint  of  problem  preparation  and  computer  running  time 
the  proposed  method  has  distinct  advantages  over  the  decomposition 
method. 

The  proposed  method  is  a  self-contained  procedure  that  will 
arrive  at  a  solution  from  the  relatively  simple  inputs  of  the  network, 
resource  productivity  data  and  movement  requirements.  The  decomposition 
approach  would  require  an  elaborate  computer  program  to  generate  all 
possible  routes  using  all  possible  combinations  of  allowable  resources. 
It  is  possible,  however,  that  a  completely  new  model  using  ohe  decom¬ 
position  principle  would  result  in  shorter  running  times.  Budget  limi¬ 
tations  prohibit  the  investigation  of  this  possibility  at  this  time. 

The  revised  RACAT  algorithm  will  increase  the  accuracy  in  the 
handling  of  capacity  limitations  if  it  is  assumed  that  a  vehicle  limi¬ 
tation  is  a  more  accurate  expression  of  capacity  than  a  limitation 
expressed  in  short  tons.  Since  most  capacities  expressed  in  short  tons 
are  simple  conversions  from  vehicle  per  time  period  figures  using  the 
average  vehicle  load,  it  seems  likely  that  this  is  true. 

The  inclusion  of  the  node  constraints  allows  for  a  direct  limi¬ 
tation  on  the  throughput  capacity  of  a  node,  rather  than  the  more 
complex  way  of  depicting  a  node  graphically  as  a  small  subnetwork. 
"Networking"  a  node  is  possible  only  in  specialized  caseB,  or  additional 
assumptions  like  directing  the  arcs  need  to  be  made.  These  additional 
assumptions  tend  to  unduly  restrict  the  network. 

The  revision  al30  permits  much  more  flexible  modeling  of  speciali¬ 
zed  processes  such  as  transfer,  depot  and  other  intra-nodal  functions. 
The  capacities  of  arcs  representing  such  processes  do  vary  depending 
on  the  commodity  involved.  The  revised  procedure  will  permit  a  more 
realistic  handling  of  these  probler- . 
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INTRODUCTION 


The  purpose  of  this  paper  is  to  study  the  use  of  post-optimal 
procedures  in  the  ETNAM  I  system  of  programs.  These  procedures  consist 
chiefly  of  two  types:  l)  parametric  linear  programming  and  2)  sensi¬ 
tivity  analysis  or  ranging. 

In  parametric  linear  programming,  a  coefficient  (usually  selected 
fran  the  objective  function  or  the  right  hand  side  constraints)  is 
treated  as  a  parameter  and  allowed  to  vary  between  specified  limits. 

As  the  parameter  changes,  the  solution  is  changed,  Whenever  necessary, 
to  keep  it  optimal. 

Sensitivity  analysis  is  done  to  find  the  range  of  a  coefficient 
for  Which  the  current  solution  remains  optimal.  This  is  a  very  useful 
tool  for  analyzing  objective  function  or  cost  coefficients  and  right 
hand  side  constraint  coefficients.  For  these  two  types  of  coefficients 
this  method  is  known  as  cost  ranging  and  right  hand  side  ranging 
respectively. 

Currently,  the  ETNAM  optimal  solution  consists  of  the  objective 
function,  the  optimal  routes  and  the  amount  shipped  an  each  route,  the 
marginal  costs  or  shadow  prices  and  various  other  specific  reports 
gleaned  from  the  basic  opt imp 1  solution.  Questions  of  the  type  "what 
happens  to  the  solution  if  I  lower  my  requirements?"  or  "what  happens 
to  the  solution  if  I  increase  a  certain  resource  inventory?"  are  usually 
handled  by  parametric  linear  programming  in  the  standard  linear  pro¬ 
gramming  codes. 

In  the  usual  linear  programming  formulation  of  the  ETNAM-type 
problem,  all  possible  routes  are  generated  and  evaluated  by  the  algorithm. 
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The  essential  difference  and  the  great  advantage  of  the  ETNAM  method 
is  that  all  possible  routes  need  never  be  generated.  Only  those  routes 
that  are  needed  are  found  by  the  shortest  chain  algorithm.  Parametric 
programming  is  a  relatively  simple  although  time  consuming  procedure 
in  the  standard  linear  program  because  the  possible  routes  can  be 
evaluated  to  determine  the  amount  of  change  in  the  parameters  required 
to  bring  about  a  change  in  the  solution.  Thus  we  can  determine  the 
limits  of  the  coefficients  by  testing  the  potential  routes.  The 
shortest  chain  algorithm  must  first  be  given  the  values  of  the  coeffi¬ 
cients  in  order  to  find  the  best  route.  What  is  a  workable  short-cut 
in  the  standard  linear  program  would  be  a  tedious  and  expensive  triad 
and  error  procedure  for  the  shortest  chain. 

In  sensitivity  analysis,  cost  ranging  also  requires  an  explicit 
knowledge  of  all  the  vectors;  but,  right  hand  side  ranging  appears 
feasible  and  may  prove  to  be  a  useful  tool  for  measuring  the  sensitivity 
of  the  right  hand  side  coefficients. 

A  mathematical  treatment  of  the  post-optimal  procedures  follows. 
This  treatment  assumes  a  rudimentary  knowledge  of  the  RACAT  algorithm. 
The  reader  should  see  the  references  at  the  end  of  this  paper  for 
background  information. 


PARAMETRIC  LINEAR  PROGRAMMING 


Consider  the  following  linear  prog  Tin  in  standard  form. 

n 

min  Z  =  V  c.x, 

u  J 

J=1 


2 


< ■*> 


subject  to  L  a,  ,x .  =  b  ,  i  =  1 ,  2 ,  ,  m 

J=1  J 

x.iO  for  j  =  1,  2,  . ..,  n. 

J 

or  in  matrix  notation 


min  Z  =  CX 


subject  to  AX  =  d 


XiO 


where 


,  C  -  ( ci>c2>  cn) 


and  denote  a  column  of  A  by  =  j 


all  S12  * ‘ *  aln 
a21  S22  ' *  *  S2n 


0*0  ••• 

ml  m2  mn 


The  criterion  for  optimality  is  that  c_B-  A.  •  C,  SO  for  all  j,  or 

B  J  J 

equivalently,  max  (c  B-  A  -  c  )  s  0.  The  shortest  chain  algorithm 

j  B  J  J 

generates  the  vector  with  maximum  oJ^A,  -  c,.  C_  is  the  cost  coeffi- 

”  O  J  <0 

cient  vector  of  the  basic  variables  and  B  is  an  raxm  matrix  consisting 

of  those  Aj  which  are  currently  in  the  basis.  Assume  for  simplicity 

that  the  first  m  variables  are  in  the  basis.  Then,  C„  =  (c, ,  c„,  ....  c  ) 

B  12  m 

'  X1  [  all  al2  ‘ ‘  alm 

X2  !  a21  a22  • • *  a2m 

XB  =  |  !  '  B  =  (A1,  A£,  ...,  Affi)  =  ,  ; 
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aml  &m2 


•  •  • 


Now 


suppose  the  first  cost  coefficient  c^  is  to  he  parametriied 


between  its  present  value  c^  and  seme  larger  value  c^.  Denote  the 
parameters  by  X.  Then 
c^Xsc^ 


To  determine  which  vector  to  enter  the  basis  when  X  is  increased, 
m 


cbb  aj  ■  Cj  ■  xhlj  + 


TJ 


I  °lh: 


ij "  V 


1=2 


where 


2J 


mj 


=  B-1A. 


m 


Increase  X  until  X  h^  +  ^ 

1=2 


h  -  c  =o.  Choose  the  A  to  enter 

^  J  O  J 


for  which 

~  III 

•>* 

V  L  cihij 

hlJ>0 

i  min 

1-2 

A  j  =  1,2,... ,n 

- 

Hence,  an  explicit  knowledge  of  all  A^  is  assumed.  As  has  been  noted  in 

u 

the  introduction,  this  is  not  the  case  in  ETNAM  I.  Rather,  each  A  is 

J 

generated  as  it  is  needed.  To  generate  all  A  would  defeat  the  purpose 

J 

of  the  RAC AT  algorithm. 

An  alternative  to  parametric  linear  programming  is  to  make  a  change 
in  the  data  without  changing  the  structure  of  the  problem  and  "crashing" 
in  the  last  optimal  basis.  However,  there  is  no  guarantee  that  this 
method  will  save  any  computer  running  time.  For  relatively  few  and 
small  changes  It  may,  for  other  changes  it  may  even  have  an  adverse 


effect. 


A  similar  analysis  of  the  right  hand  side  parametric  linear  pro¬ 


gramming  case  has  been  made  with  the  same  conclusions  being  drawn. 


SENSITIVITY  ANALYSIS 


In  cost  ranging,  the  range  of  a  cost  coefficient  c^  is  determined 

for  which  the  basis  need  not  be  changed;  that  is,  for  which  the  solution 

remains  optimal.  The  c  ,  as  used  in  ETNAM,  is  a  function  of  the  tolls 

J 

"fch 

on  the  arcs,  the  arcs  used  in  the  j  1  route  (A,),  the  resource  produc- 

J 

tivities,  and  the  resource  prices.  So  that  the  range  of  a  resource 

price  or  a  toll  is  determined  by  the  range  of  the  c  's. 

J 

Suppose  the  coefficient  to  be  ranged,  c^,  is  in  the  basis.  Then 

the  following  inequalities  determine  the  range  of  c^.  CgB_1A^  -  C^-sO 

for  all  j  such  that  A  is  not  a  basic  variable.  This  means  that  all 

J 

possible  routes,  A^,  would  have  to  be  known. 

If  c^  is  not  in  the  basis,  then  CBB_1A^  -  c^sO  determines  the 
range  of  ck>  But,  since  is  not  in  the  basis,  has  not  yet  been 
generated  and  A^  is  an  unknown  making  it  impossible  to  solve  this 
inequality  for  c^. 

Since  not  all  the  A  are  known  in  the  ETNAM  formulation,  cost 
J 

ranging  cannot  currently  be  done.  However,  right  hand  side  ranging 
can  be  done. 

To  find  the  range  of  b^,  the  followirg  set  of  inequalities  must 
be  solved. 

B  or  writing  them  out , 

*  ■■■  *  Vk  *  •"  + 
s2ibi  *•••♦  S2k\  *  *  s^V0 


fL->b,  +  •••  +  8  iA,  +  •••  +  B  b 
ml  1  mk  k  mm  m 
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where 


-1 


Bn  B 


11  p12  *•*  hm 

P21  022  p2m 


Jml  pm2  0mm 


th 


Pffll  & 

The  i  *  Inequality  can  he  re-written  as 

- (ftjl^l  +  +  Pik-1  bk-l  +  ftjk+1  bk+l  +  +  b 


V 


im  m) 


P 


ik 


for  pik>0 


and 


V 


-(Pilbl  +  •“  *  ftjk-1  Vl  +  ftjk+1  bk+l  +  •••  +  Pimbm) 


P. 


ik 


for  ftik<0 


Denote  the  upper  and  lower  bounds  of  b^,  by  b^  and  b^  respectively. 
Then 


b,  =  max 
-*  i 


b.  =  min 
k  i 


-(Pjl  bl  *  ftjk-1  bk-l  +  eik^  x+1  +  Pimbm) 


P 


ik 


^ik"0 


^ft-«i  bi  +  +  Pik-1  bk-l  +  Pik+1  bk+l  +  b- 


“fil  1 


im  m) 


'ik 


Pik<° 


The  inequlalty  for  which  0  k  =  0  contributes  nothing  to  the  upper  or 
lower  bound  of  the  range  of  b^. 

It  has  been  shown  here  that  right  hand  side  ranging  can  be  done 
in  ETNAM  I. 


CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  shown  in  the  above  analysis  that  it  is  possible  to  dc 

right  hand  side  ranging  in  the  EINAM  model.  But  within  the  present 

framework  of  ETNAM  I,  parametric  linear  programming  and  cost  ranging 

cannot  be  done.  Based  upon  this  analysis  it  is  feasible  to  implement 

right  hand  side  ranging  in  the  ETNAM  system  of  computer  programs.  This 

6 


would  allow  the  user  to  analyze  the  sensitivity  of  arc  capacities, 
resource  inventory  constraints,  requirements  and  any  other  constraint 
that  may  be  added  at  some  future  date. 

The  programming  task  does  not  appear  too  great  since  it  is  possible 
to  write  the  right  hand  side  ranging  routine  as  a  separate  program 
module  with  only  minor  changes,  if  any,  to  the  other  three  program 
modules  of  ETNAM.  Virtually  all  of  the  available  computer  core  is  being 
used  in  the  RACAT  and  (XJTRJT  modules .  Incorporating  this  post  optimal 
procedure  in  either  of  these  modules  would  necessitate  a  reallocation 
and  economization  of  computer  core.  This  would  also  make  the  debugging 
more  difficult. 

It  is  therefore  recommended  that  a  separate  ETNAM  program  module  be 
written  to  do  right  hand  side  ranging,  ics  use  would  be  optional  to  the 
user. 
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SAMPLE  LISTINGS  QF  DATA  FOR  INPUT  MODULE 


1.  Sample  Problem.  The  first  listing  is  the  input  data  of  the  actual 
reports  of  the  test  problem  shown  in  the  SD,  Volume  I,  User's  Manual. 
There  are  no  substitution,  4  modes.  2  commodities,  5  resources,  2  node 
constraints.  Link  capacities  are  expressed  in  vehicles-per-time-period. 

2.  Partial  Listing.  The  second  listing  emits  the  complete  network. 

This  sample  shows  2  substitutions,  4  modes,  2  commodities,  J  resources, 
no  node  constraints.  Link  capacities  are  expressed  as  vehicles-per- 
time  period.  The  productivity  cards  PROD  and  PALD  are  listed  shewing 
the  handling  of  substitution  of  resources.  The  payload  factor  cards, 
TRIP,  also  show  the  handling  of  substitution  of  resources. 

3*  Partial  Listing.  The  last  listing  also  omits  the  complete  network. 
This  sample  shows  no  substitution  of  resources,  4  modes,  2  commodities, 

5  resources,  no  node  constraints.  Link  capacities  are  expressed  as 
tons-per-time-period.  The  conversion  factor  card  CONF  is  shown.  The 
origin-destination  nodes  ODPR  have  open  delivery  requirement  so  the 
requirements  by  commodity  cards  DELI  are  shown.  There  is  also  a  sample 
TOLL  card. 


Enclosure  5 


E5-1 


. 


...ETNAM  I  SAMPLE  PROBLEM 

...  SEPTEMBER  1970 

1 

15  tO 

16 

90 

95 

to 

is 

50 

} 

+ 

i 

+ 

1 

* 

MIN  COST 

1  L 

2 

5 

2 

1 

75 

150 

1 

10 

16 

it 

t8 

n 

\ 

♦ 

NAMO 

HINT 

RAIL 

TRAN 

PORT 

NACO 

DRY 

PAX 

NAPE 

S/PL 

BUS 

FRTC 

PAXC 

TSC 

APCS 

1 

7 

18 

n 

91 

w 

50  51 

f 

f 

t 

t 

t 

H 

OlH 

OSH 

1 

i 

129 

700 

in 

OlH 

11H 

1 

i 

90 

700 

in 

02H 

03H 

1 

i 

120 

900 

in 

02H 

06H 

.1 

i 

190 

200 

m 

G3H 

OLH 

1 

i 

90 

800 

in 

0?H 

OSH 

1 

i 

140 

700 

in 

04H 

OSH 

1 

l 

190 

700 

in 

95H 

OSH 

1 

l 

290 

300 

in 

05H 

0  2H 

1 

l 

160 

400 

in 

05H 

11M 

1 

l 

105 

400 

in 

07H 

OSH 

1 

l 

SO 

700 

in 

07H 

11H 

1 

l 

85 

700 

in 

OSH 

09H 

1 

l 

165 

400 

hi 

OSH 

02H 

1 

i 

190 

*•00 

hi 

09H 

OLH 

t 

i 

220 

400 

in 

10H 

OLH 

1 

l 

379 

50 

hi 

11H 

02H 

1 

l 

125 

400 

in 

11H 

09H 

1 

l 

150 

500 

in 

OIR 

HR 

2 

l 

95 

700 

in 

OIP 

12R 

2 

l 

220 

700 

in 

02R 

0  3R 

2 

l 

135 

700 

in 

02P 

OSR 

2 

l 

215 

30 

in 

03R 

OLR 

2 

l 

95 

700 

hi 

OLR 

OSR 

2 

l 

215 

700 

in 

OTP 

13R 

2 

l 

70 

700 

hi 

OOP 

02R 

2 

l 

200 

700 

in 

OSP 

10R 

2 

l 

110 

700 

in 

09R 

OLR 

2 

l 

250 

700 

in 

10P 

OLR 

2 

l 

410 

50 

hi 

IIP 

02R 

2 

l 

140 

700 

in 

12P 

Q3R 

2 

l 

170 

*03 

in 

hr 

02R 

2 

l 

140 

iOO 

in 

13R 

09R 

2 

i 

115 

500 

in 

07R 

OSR 

2 

l 

65 

7M0 

in 

OlH 

OIR 

3 

l 

1 

1 

in 

02H 

02R 

3 

l 

1 

* 

in 

03H 

03R 

3 

l 

l 

1 

in 

OLH 

OLR 

3 

l 

1 

1 

in 

06H 

06R 

3 

l 

1 

1 

in 

07H 

07R 

3 

l 

1 

1 

in 

OSH 

OSR 

3 

l 

t 

1 

hi 

09H 

09R 

3 

l 

1 

1 

in 

10H 

10R 

3 

1 

1 

1 

111 

04H 

550 

3 

1 

1 

1 

1101 

04H 

658 

3 

1 

1 

1 

1011 

06H 

55D 

3 

1 

1 

1 

1101 

06H 

650 

3 

1 

1 

1 

1011 

50S 

01H 

1 

4 

t 

1 

1 

1101 

50S 

07H 

1 

4 

1 

1 

1 

1101 

60S 

07H 

4 

1 

1 

1 

1011 

60S 

08H 

1 

4 

1 

1 

1 

1011 

60S 

10H 

1 

4 

1 

1 

1 

1011 

(BLANK  CARO) 

1 

10 

SO 

f 

t 

♦ 

03R 

2 

500 

HR 

2 

180 

/ 

10 

15 

SO 

30 

40 

45 

50 

60 

\ 

t 

i 

f 

{ 

1 

INVE 

1 

100000 

2 

100000 

INVE 

3 

100000 

4 

100000 

INVE 

5 

100000 

PHIS 

1 

7 

2 

7 

PHIS 

3 

1 

4 

1 

PHIS 

5 

2 

RATE 

1 

15 

2 

12 

RATE 

3 

.1 

4 

.2 

PROO 

* 

1 

1 

IE-3 

1 

2 

3 

1P4E-5 

PROD 

i 

3 

5 

3E-5 

1 

4 

5 

139E-5 

PROD 

2 

1 

2 

25E-5 

2 

2 

4 

8286E-8 

PROO 

2 

4 

5 

143E-8 

TRIP 

1 

1 

1 

10 

1 

1 

2 

10 

TRIP 

1 

1 

3 

7000 

1 

1 

4 

10000 

TRIP 

1 

2 

3 

14000 

1 

2 

4 

7000 

TRIP 

1 

2 

1 

40 

1 

2 

2 

50 

1 

10 

n 

SO 

35 

40 

44 

50 

65 

t 

nnDD 

f 

t 

t 

t 

1 

1 

t 

1 

UUPK 

OOPR 

1 

50S 

55D 

9080 

2 

60S 

65D 

13475 

OOPR 

999 

ENO 

THERE  IS  ONLY  ONE  BLANK  CARD  IN  THE  DATA,  OTHERS  ARE  FOR  CLARITY 
CARO  COLUMNS  ARE  INOICATEO  BY  ARROWS  WITH  COLUMN  NUMBERS  ABOVE  THEM 


E5-3 


• ETNAM  I  SAMPLE  LISTING  WITH  2  SUPST 


1 

15  20 

25 

30 

35 

1 

1  1 

1 

; 

♦ 

MIN  COST 

2  4 

2 

7 

0 

1 

10 

16 

•22 

28 

34, 

1 

1 

1 

I 

♦ 

1 

NAMO 

LIWY 

RAIL 

TRAN 

PORT 

ttACO 

DRY 

PAX 

NARE 

ARCS 

S/PL 

BUS 

FRTC 

PAXC 

TSC 

NO  NOOE  CONSTRAINTS,  VEMICLES/TIME  PERIO 

40  45  50 

t  A  tA 

1*0  46 

I  I 

SCAR  CCAR 


/ 

7 

18 

n 

32 

40 

50  51 

1 

1 

1 

1 

1 

1 

II 

01H 

OSH 

1 

1  125 

700 

111 

01H 

11H 

1 

1 

90 

700 

111 

02H 

0  3H 

1 

l  120 

900 

111 

02K 

06H 

1 

1  IPO 

200 

111 

03H 

04H 

• 

1 

1 

90 

800 

111 

60S 

• 

07H 

4 

1 

1 

1 

101 

60S 

OSH 

4 

1 

1 

1 

101 

60S 

10H 

4 

1 

1 

1 

101 

(BLANK  CARO) 

1 

10 

15 

20 

30 

35 

40 

45 

50 

60 

65 

1 

1 

1 

1 

f 

1 

1 

♦ 

1 

1 

t 

INVE 

1 

inoooo 

2 

100000 

INVE 

3 

2000 

4 

50 

INVE 

5 

inoooo 

6 

100000 

INVE 

7 

1C0000 

PHIS 

1 

7 

2 

7 

PHIS 

3 

1 

4 

1 

PHIS 

5 

2 

6 

1 

PHIS 

7 

1 

RATE 

1 

15 

2 

12 

RATE 

3 

.1 

I- 

.2 

PROD 

1 

1 

1 

IE-3 

1 

2 

3 

104E-5 

PROD 

1 

1 

51 

IE-3 

1 

2 

57 

104F-4 

PROD 

2 

1 

2 

25F-5 

2 

7 

6 

8286E-7 

PRC  J 

2 

1 

52 

25E-5 

2 

2 

54 

8286E-8 

PROD 

2 

4 

5 

143E-8 

2 

4 

55 

143E-8 

PALO 

1 

3 

5 

33333 

1 

1 

4 

5 

719 

l 

PALD 

1 

3 

55 

33333 

1 

1 

4 

55 

^19 

1 

TRIP 

1 

1 

1 

10 

1 

1 

2 

10 

TRIP 

1 

1 

3 

7000 

l 

1 

4 

10000 

TRIP 

1 

2 

3 

14000 

1 

7 

4 

7000 

TRIP 

1 

2 

t 

43 

1 

2 

2 

30 

TRIP 

2 

1 

2 

• 

2 

** 

2 

50 

i 

10 

11, 

20 

35 

40 

ff 

•50 

65 

1 

1 

1 

1 

1 

1 

1 

| 

| 

ODPP 

OOPP 

1 

SOS 

550 

9050 

2 

f.OS 

658 

13475 

09PP  999 
END 


THERE  IS  ONLY  ONE  BLANK  CARO  IN  THE  DATA,  OTHERS  APE  FOR  CLARITY 
CAPO  COLUMNS  ARE  INDICATED  BY  ARROWS  WITH  -.OLUMN  NUMBERS  ABOVE  THEM 


.ETNAM  I  SAMPLE  LISTING  WITH  1  SUBST,  MO  NODE  CONSTRAINTS,  TONS 


1 

15 

20 

25  30 

35 

50 

55 

50 

t 

♦ 

t  t 

1 

+ 

i 

5 

MIN  1 

DOST 

1 

4 

2  5 

0 

0 

75 

150 

1 

10 

16 

22 

28 

35 

1 

t 

1 

NAMO 

HIWY  RAIL  TRAN  PORT 

NACO 

ORY 

PAX 

NARE 

S/PL  BUS 

FRTC  PAXC 

TSC 

ARCS 

1 

7 

18 

2£  32 

50 

50  51 

1 

1 

t 

\  \ 

t 

H 

01H 

OSH 

1 

1  125 

7000 

in 

01H 

11H 

1 

1  90 

7000 

in 

02H 

03H 

1 

1  120 

9000 

in 

02H 

06H 

1 

1  190 

2000 

in 

03H 

04H 

• 

1 

1  90 

6000 

in 

60S 

• 

07H* 

4 

1  1 

10000 

101 

60S 

08H 

4 

1  1 

10000 

101 

60S 

10H 

4 

1  l 

10000 

101 

(BLANK  CARO) 

1 

10 

15 

20 

30 

50 

55 

50 

60 

\ 

1 

t 

t 

5 

t 

t 

\ 

t 

IN  VE 

1 

100000 

2 

100000 

INVE 

3 

100000 

4 

100000 

INVE 

5 

100000 

PHIS 

1 

7 

2 

7 

PHIS 

3 

1 

4 

1 

PHIS 

5 

2 

RATE 

1 

19 

2 

12 

RATE 

3 

.1 

4 

.2 

PROO 

1 

1 

1 

IE-3 

1 

2 

3 

104E-5 

PROD 

1 

3 

5 

3E-5 

1 

4 

5 

139E-5 

PROO 

2 

1 

2 

25E-5 

2 

2 

4  l 

B286E-8 

PROD 

2 

4 

5 

143E-8 

CONF 

2 

5 

DELI 

1 

9080 

2 

13475 

1 

10 

n 

20 

35 

50 

55 

50 

5 

1 

5 

\ 

5 

t 

\ 

5 

OOPR 

OOPR 

1 

SOS 

03H 

1 

50S 

12R 

OOPR 

2 

60S 

06H 

2 

60S 

03H 

OOPR 

999 

TOLL 

TOLL 

1 

03H 

04H 

100000 

TOLL 

999 

ENO 

THERE  IS  ONLY  ONE  BLANK  CARO  IN  THE  OATA,  OTHERS  ARE  FOR  CLARITY 
CARO  COLUMNS  ARE  INDICATED  BY  ARROWS  WITH  COLUMN  NUMBERS  ABOVE  THEM 


E5-5 


SELECTED  BIBLIOGRAPHY 


General  Background  Information 

Bug acker,  R.  G.  and  Saaty ,  T.  L. ,  Finite  Graphs  and  Networks ,  McGraw- 
Hill  Book  Company,  New  York,  N.  Y.  196*5. 

Cremeans,  J.  E.,  "The  Bridges  of  Konigsturg, "  ASTME  Vectors,  May-June, 

1969,  p.  4. 

Hadley,  G.,  Linear  Programming,  Addison-Wesley  Publishing  Company,  Inc., 
Reading  Mass.,  1962. 

Ford,  L.  R.  Jr  and  Fulkerson,  D.  R.,  Flows  in  Networks ,  Princeton 
University  Press,  Princeton,  N.  J.  1962. 

Orchard-Hays  W. ,  Advanced  Linear  Programming  Computing  Techniques, 
McGraw-Hill  Book  Co.,  Inc.,  New  York,  19^8. 

Theoretical  Papers 

Cremeans,  J.  E.  and  Tyndall,  G.  R.,  "A  Model  for  Optimal  Multi- Commodity 
Network  Flows  with  Resource  Allocation,"  Research  Analysis  Corp., 
RAC-P-44,  October  1968. 

Cremeans,  J.  E.,  Smith,  R.  A.,  and  Tyndall,  G.  R.,  "An  Extension  to  the 
Multi “Commodity  Network  Flow  Model  to  Permit  Substitution  of 
Resources,"  Research  Analysis  Corp.,  RAC-P-52,  April  1969* 

Cremeans,  J.  E.,  Smith,  R.  A.,  and  tyndall,  G.  R. ,  "Optimal  Multicommodity 
Network  Flows  with  Resource  Allocation,"  Naval  Research  Logistics 
Quart.,  17(3):  Sept.  1970. 

Dreyfus,  S.  F.,  "An  Appraisal  of  Some  Shortest  Path  Algorithms,"  The 
RAND  Corp.,  RM-5433-1-PR,  Sept.  68. 

Ford,  L.  R.  Jr  and  Fulkerson,  D.  R.,  "A  Suggested  Computation  for 
Maximal  Multi-Commodity  Network  Flows,"  Mgt.  Sci.,  Oct  58. 

Tomlin,  J.  A.,  "Minimum-Cost  Multi -Commodity  Network  Flows,"  Opns.  Res., 
Dec  66. 

Additional  References 


Boyer,  Donald  D. ,  "A  Modified  Simplex  Algorithm  for  Solving  the  Multi- 
Commodity  ftoimum  Flow  Problem,"  TM-14930*  The  George  Washington 
University  Logistics  Research  Project,  Washington,  D.  C.,  Mar  68. 

Busacker,  R.  G.,  et  al,  "Three  General  Network  Flow  Problems  and  Their 
Solutions,"  RAC-SP-I83,  Research  Analysis  Corporation,  Nov  62. 

Charnes,  A.,  "Optimality  and  Degeneracy  in  Linear  Programming," 
Econometrica,  Vol  20,  April,  1952. 

Fitzpatrick,  G.  R.,  et  ad .  "Programming  the  Porcurement  of  Airlift  and 
Sealift  Forces:  A  Linear  Programming  Model  for  Analysis  of  the 
Least  Cost  Mix  of  Strategic  Deployment  Systems,"  Naval  Research 
Logistics  Quart.,  14  (2):  (1967). 


Enclosure  6 


E6-1 


Preceding  page  blank 


Additional  References 


Koopmans,  T.  C.,  ed.,  Activity  Analysis  of  Production  and  Allocation, 
John  Wiley  and  Son,  Inc.,  New  York,  N.  Y.  1951* 

Rao,  M.  R.  and  S.  Zionts,  "Allocation  of  Transportation  Units  to 

Alternative  Trips  --  A  Column  Generation  Scheme  with  Out-of-Kilter 
Subproblems,"  Qpns .  Res . .  Jan -Feb  68. 

Sakarovitch,  M. ,  "The  Multi-Commodity  Maximum  Flow  Problem, "  Operations 
Research  Center,  University  of  California,  Berkeley,  Dec  66. 


